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ABSTRACT 


This  report  is  a  compilation  of  coordinate  system,  equations,  and  general  relations  need  in  aircrafts 
motion  analysis.  The  information  presented  is  useful  for  comparative  evaluation  and  for  preliminary-design 
work.  Simplified  and  approximate  solutions  are  given  for  special  flight  conditions.  Ite  material  is  defined 
and  presented  m  a  form  suitable  for  direct  application.  Derivation  and  theoretical  development  are  not 
emphasised,  but  sources  thereof  are  named. 

This  technical  documentary  report  has  been  reviewed  and  is  approved. 
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SECTION  1.  INTRODUCTION 


Data  and  information  are  presented  in  this  report  for  use  in  the  analysis  of  aircraft  motion.  TTiis  report,  which 
was  originally  intended  to  be  a  part  of  the  11SAF  Stability  and  Control  Handbook  published  ;u  October  1960,  is  a 
compilation  and  condensation  of  the  coordinate  systems,  equations,  and  general  information  related  to  aircraft- 
motion  analysis.  The  original  form  remains  essentially  unchanged. 

The  purpose  of  the  Handbook  was  to  provide  the  data,  equations,  and  relations  necessary  to  analyze  the  mo¬ 
tions,  the  stability  and  control  characteristics,  and  flying  qualities  of  aircraft  in  concise,  consistent,  and  readily 
usable  form.  In  keeping  with  this  purpose,  emphasis  is  placed  on  description,  definition,  and  application  rather 
than  on  derivation  and  theoretical  development.  Problems  of  unusual  nature  and  unconventional  configurations 
may  require  special  analysis  and  development  of  particular  equations  from  the  fundamental  theory  cited  in  the 
references. 

The  basic  kinematic  and  dynamic  relations  for  particle  and  rigid-body  motion  are  included.  Several  convenient 
coordinate  systems  ar<.  Jefleed,  and  coordinate  transformation  relations  are  given.  Force  aud  moment  components 
are  developed,  and  a  compilation  of  conventional  stability  derivatives  is  presented.  The  rigid-body  equations  of 
motion  are  simplified  for  special  flight  conditions,  and  some  approximate  solutions  are  given.  Some  material  is 
presented  pertaining  to  instrument  readings  and  fuel  slosh. 

Symbols  and  nomenclature  are  listed  and  defined  in  the  sections  to  which  they  apply-  Consistency  in  symbols 
and  notation  is  maintained,  except  in  cases  where  established  usage  dictates  otherwise.  A  complete  list  of  sym¬ 
bols  is  not  considered  to  be  necessary  and  is  not  given. 
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SECTION  2.  COORDINATE  SYSTEMS  AND  EQUATIONS  OF  MOTION 


In  order  to  describe  the  motion  of  a  dynamic  system  it  is  necessary  to  define  a  suitable  coordinate  system  and  formulate 
equations  for  the  motion  in  accordance  with  the  physical  laws  governing  the  system. 

The  diagrams  and  discussion  that  follow  consider  the  motion  of  a  particle  (point  mass)  and  the  more  complicated 
motion  of  a  rigid  body. 


PARTICLE  MOTION 

Coordinate  systems  and  equations  that  conveniently  describe  the  motion  of  a  point  mass  are  presented  in  the  following 
pages.  Rectangular,  spherical,  and  cylindrical  coordinate  systems  are  presented.  Preferred  axis  orientation  and  notation 
indicated  and  used  are  consistent,  insofar  as  possible,  with  the  reference  literature. 


RECTANGULAR- COORDINATE  SYSTEM  (FLAT  NONROTATING  EARTH) 


The  familiar  Cartesian  or  rectangular  coordinate  system  has  many  applications  in  the  analysis  of  vehicle  motion.  For 
instance,  it  may  be  used  to  describe  the  Sight  path  (trajectory)  of  a  body  with  respect  to  a  given  starting  point  on  the 
earth’s  surface.  A  typical  case  is  suggested  in  the  description  of  the  coordinate  system  below.  Generalization  to  any 
specific  problem  is  self-evident  and  requires  no  further  discussion. 

Description  of  Coordinate  System 

Origin  of  rectangular  coordinates  x,  y,  z:  arbitrary,  often  a  point  on  the  surface  of  the  earth. 

Fundamental  plane:  usually  the  XY-plane;  tangent  to  the  surface  of  the  earth  at  the  origin. 

Positive  X-axis:  arbitrary,  often  selected  along  initial  heading  or  direction  of  motion. 

Positive  Z-axis:  arbitrary,  often  oriented  in  sense  to  denote  altitude  above  the  surface  of  earth  or  the  XY-plane. 

Positive  rotation  in  fundamental  plane:  from  X-axis  to  Y-axis;  i.e.,  right-hand  system. 
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NOTATION 


U,k 

OP 

x.y,* 

v 

v„  VT,  V. 

m 

(*) 

Vector  form: 


orthooormsl  base  (unit)  rector*  shag  X-,  Y-,  and  Z-asrt,  rwpectWdy 
position  vector  of  point  P  (rectangular  coordinator  x,  y, «) 

position  coordinate*  of  P;  also  components  of  OP  aleqg  coordinate  cnee.  La,  OP=  d  +  jJ  +  * 
velocity  vector  of  point  P 

components  of  velocity  V  along  the  coordinate-axis  directions,  ia,  V  =s  Vj  +  V»|  +  VA 

mass  of  particle  at  point  P 

denotes  differentiation  with  respect  to  time 

Equations  of  Motion 


F  =  m  —  =  nV 
dt 


(1) 


Component  form: 


F,  =  mVt  =  iM 
F,  =  m  Vr  =  miy 
F«  =  m  V,  =s  na 


(2) 
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SPHERICAL-COORDINATE  SYSTEM 

The  analysis  of  motions  within  the  inertial  frame  fixed  to  the  center  of  die  earth  is  moat  conveniently  treated  in 
spherical  coordinates.  This  section  considers  both  rotating  and  nonrotating  spherical  coordinates.  In  order  to  distinguish 
between  these  two  systems,  primed  quantities  refer  to  nonrotating  coordinates  and  imprinted  quantities  refer  to 
rotating  coordinates.  Since  it  is  customary  to  refer  our  position  and  velocity  to  the  earth,  the  rotating  coordinates 
are  generally  used. 

Flight-path  coordinates  are  introduced  because  aerodynamic  forces  are  frequently  considered  in  the  analysis  of  a 
vehicle  flight  path.  Aerodynamic  forces  are  most  conveniently  related  to  the  velocity  of  the  vehicle  through  the  air, 
which  rotates  with  the  earth.  Thus  the  rotating -earth  flight-path  coordinates  may  be  used  in  the  analysis  of  missile 
and  supersonic-  or  hypersonic-vehicle  flight  paths  whenever  aerodynamic  forces  are  included. 

The  basic  development  of  the  equations  of  motion  in  this  Section  is  given  in  reference  1. 

Description  of  Coordinate  System  (reference  2) 

Origin  of  spherical  coordinates  r,  +,  $:  canter  of  the  earth. 

Fundamental  plane:  equatorial  plane. 

Reference  direction  in  fundamental  plane:  arbitrary,  e^,  Greenwich  Meridian  ased  for  longitude  reference. 
Polar-axis  positive  direction:  toward  the  North  Pole. 

Positive  rotation  in  fundamental  plane:  eastward.  La,  a  right-hand  system. 

NOTATION 

0  origin  of  routing  spherical  coordinate  system,  center  of  the  earth 

P  particle  under  coaaideretiea 
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r 


4 


OP 

r 

* 

$ 

0 

Ir,  1*  1. 

Y 

V 


V„  V*  V. 
*•»  £•«  ®* 


< 

7 

7 

F 


F„F*F. 
F.,  F.,  F, 


m 

( > 

Notes: 


position  vector  of  P  (spherics!  coordinate*  r,  #> 
radial  distance  from  origin  to  P 
angular  inclination  of  P  to  polar  axis 

angular  displacement  of  P  from  reference  meridian  plane  (reference  plane  rotates  with  the  earth) 

angular  velocity  of  rotating  coordinate  system  about  the  polar  axis 

orthonormal  base  (unit)  vectors  at  P  along  spherical-coordinate  directions 

velocity  vector  of  P  with  respect  to  rotating  coordinates 

magnitude  of  velocity  vector  V 

components  of  velocity  V  along  the  coordinate  directions,  Len  V  =  Vrl,  -j-  Y,l»  +  V#l# 

orthonormal  base  (unit)  vectors  at  P  along  flight-path  coordinates  (e,  k  aligned  along  the  velocity  vector  V ;  e.  and  e. 
are  oriented  normal  and  sidewise,  respectively,  to  the  flight  path) 

flight-path  heading  -  the  angle  between  the  meridian  plane  through  P  and  the  flight  plane  determined  by  the  radios 
vector  OP  and  velocity  vector  V.  The  angle  I  k  measured  clockwise  from  the  poL?  axis. 

flight-path  attitude  —  the  angle  between  the  velocity  vector  V  and  the  local  ooriaon  at  P 

roll  angle  —  the  angular  displacement  of  the  base  vector  e.  from  the  flight  plane  containing  OP  and  V.  The  angle  y 
k  positive  in  the  sense  of  a  right-hand  rotation  about  e,. 

real  force  vector  applied  at  P 

components  of  F  along  spherical-coordinate  directions,  Le^  F  =  F,l,  +  F#l«  +  Fsl, 
components  of  F  along  flight-path-coordinate  directions,  is,  F  =  F«e.  +  F»e«  -f  F,*, 
mass  of  particle  P 

denote*  differentiation  with  respect  to  time. 

1.  Notation  .for  nonrotating  spherical  coordinates  k  the  same  as  above  with  the  addition  of  a  prime. 

Z  OP  =  r  1, 

* 

3.  V  =  Ve, 


The  equations  of  motion  for  a  particle  mass  moving  in  the  inertial  frame  fixed  at  the  center  of  the  earth  are  given  below. 
These  equations  are  derived  from  the  basic  vector  equation  for  the  motion, 


£ 

m 


dV'  dV 
dt  ~  dt 


+  (flXV) 


It  is  important  to  note  that  V'  is  the  inertial  velocity  measured  with  respect  to  nonrotating  coordinates,  and  that 
V  is  the  relative  velocity  measured  with  respect  to  the  rotating -coordinate  system. 

The  equations  of  motion  in  nonrotating  coordinates  may  be  obtained  by  simply  considering  the  angular  velocity  O  to 
be  aero  in  the  equations  for  the  rotating-coordinate  systems  fin  which  caae  V  ==  V'). 

The  relation  between  spherical  coordinates  and  flight-path  coordinates  is  given  by  the  following  rotation  of  the  base 
vectors  lr,  If  1#  at  P.  (refer  to  figure  2). 

a.  Rotate  1#  about  1,  through  angle  90—1  to  the  flight  plane  determined  b j  OP  and  V. 

b.  Rotate  in  the  flight  plane  through  an  angle  y  such  that  the  base  vector  1#  ».  w  coincides  with  the  velocity  V. 
This  base  vector  is  then  noted  as  e,. 

c.  Rotate  in  roll  about  eT  through  an  angle  y  to  the  final  orientation  of  die  flight-path  base  vectors 
e»»  e«  ®v- 
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1.  ROTATING'EARTH  SPHERICAL  COORDINATES  (REFERENCE  1) 


L_  Vf+W  +  Qrna*)* 

m  “  V'“ 

=  f  —  r(^)*  —  r  (i  -j-  Q)8  sin8  4 

I±=1r  4.  V'  V*  (V<-f  rQ.in»)8 
m  *  r  r  tan  <f> 

—  -j-  2f  4  —  r(l+°)*8*n^C08^ 


*=v. 

m 


(Vr+2 


=  r  ?sin^  -}-  2r£  (0  -f  O)  ccs^  -}-  2r  (0  -{-  n)  sin <f> 

2.  ROTATING-EARTH  FLIGHT-PATH  COORDINATES 

f  r  .  v8 

— -  —  I  V  y - cos  y  —  r  Q8  sin  4  (cos  y  sin  <f>  ~f-  sin  y  cos  8  cos  <£) 

m  L  r 


2  O  V  sin 


1  r  *  V* 

in  8  sin  cos  if- 1-  |^V  8  cos  y  —  -■  ^  ^-cos3  y  sin  8  cos^ 

—  r  Q3  sin  8  sin  ^  cos  ^  —  2  O  V  (cos  y  cos  $  —  sin  y  cos  8  sin  sin  1 j 

F  f  • 

— -  =  I  V  8  cos  y  —  r  O8  sin  8  sin  4  cos  4  -|-  2  Q  V  (sin  y  cos  8  sin  ^ 

ra  L 

>] 


Vs 


—  cosy  cos  ^)  — — : — -  cos8  y  sin  8  cos  ^  I  cos  if 
r  sin  9  1 


fv3 


+Lt  cos  y  —  V  y  2  n  V  sin  8  sin  ^ 


-{-  r  n8  sin  ^  (cos  y  sin  ^  -f-  ®*n  y  cos  8  cos  <£)  J  sin  1 j 


Mv 

m 


in  y  sin  $)J 


(3) 


(4) 


V  -f-  r  D3  sin  ^  (cos  y  cos  8  cos  ^  —  sin 


Note:  The  equation  for  Ft  (side  force)  may  be  solved  for  the  roll  angle  ij  such  that  F,  =  0.  The  equations 
then  become  the  equations  for  motion  in  coordinates  similar  to  symmetric  wind  axes.  The  angle  *>  is  then 
the  bank  angle  (in  aircraft  terminology)  required  for  flight  with  aero  sideslip. 

3.  NONROT ATING-EARTH  SPHERICAL  COORDINATES 

JL=y  .(V^-hrV) 

m  '  r 

=  ?_r(J)»_r(j')W# 

it.=  v«,+yi^_-zk 

m  *  ~  r  r  tan  4 
=  rqk  -+■  2rj  —  r  (£)8  sin  ^  cos  + 

F,  ,  VVV%  ,  VVVV 

17- V#+  r  +rUn^ 

=  r0'sin*-j-2r4#,cos4-f2r0'sinf 


(5) 
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4.  NONRGTATING-EARTP  FLIGHT-PATH  COORDINATES 
Jk  — f"v' 


[\"  y  —  —cMy'JcMif'  +^f'kct»Y-  ~  —  6os*  y  sin  S’  coa  4  J  sin  1 j' 

—  =  F V'  i  '  cm  y'  —  CM*  y'  sin  S'  cm  + 1  cm  if 

in  i  r  sin  4  I 

+  r^-  cos  y'  -  V'  y'  J  sin  vf 


(6) 


CYIJNDRICAL-FOLAR  COORDINATE  SYSTEM 

Cylindrical  coordinates  have  limited  applications  in  particle-motion  analysis.  They  are  however  conveniently  used 
in  many  problems  of  planar  motion  where  perturbations  perpendicular  to  the  fundamental  plane  are  considered. 

The  general  cylindrical-coordinate  system  is  defined  and  illustrated  in  fignre  3.  Two-dimensional  polar  coordi¬ 
nates  are  a  special  case  of  cylindrical  coordinates  where  the  ^-coordinate  is  held  constant. 

Description  of  Coordinate  System 
Origin  of  cylindrical  coordinates  r,  $,  s:  arbitrary. 

Fundamental  plane:  the  reference  plane  normal  to  the  polar  (Z lexis. 

Reference  direction  in  the  fundamental  plane:  arbitrary. 

Polar  (Z)sxis  positive  direction  and  rotation  in  die  fundamental  plane:  right-hand  system. 


pol  an  axis 

s 


la 


m  PtmDAMXKTAL  PLANS 


PIOUU  S  OTURBIIOAMOOIOIIATI  imM 
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NOTATION 


0  origin  of  eoenHsme  «yit«m 

P  denotes  particle  wader  eonsldention 

OP  pooftka  vector  sf  P  (cylindrical  coordinates  r,  0, «) 

r  radial  distance  of  P  from  polar  ask 

$  aagalar  dkplatsmrat  of  P  about  polar  ask  from  reference  direction 

x  displacement  of  P  from  the  fundamental  plane 

1„  1*,  1.  orthonormal  base  vectors  at  P  along  cylindrical -coordinate  directions 
V  velocity  vector  of  P  with  respect  to  coordinate-ask  system 

V„  V#,  V,  components  of  velocity  V  along  cylindrical-coordinate  directions,  Le.,  V  =  V,l,  -f  V»l#  -f  V,l« 
F  reel  force  vector  applied  at  P 

Fr,  F»,  F.  components  of  F  along  cylindrical-coordinate  directions,  i.en  F  =  F,l,  +  F*l»  +  F,l, 
m  mass  of  particle  P 

( *)  denotes  differentiation  with  respect  to  time 

Note:  OP=rlr  +  iH. 


For  i  particle  moving  in  an  inertial  frame  the  equations  of  motion  expressed  in  cylindrical  coordinates  are  as  follows: 

Vr~  ^  =V-r<0)2 
m  -r 


v.  .  VV. 

~  =  V#-f  1*n  =  T6  +  2r6 
in  r 


m 


=  V,  =  I 


(7) 


As  was  noted  previously,  the  first  and  second  equations  above  may  be  used  for  planar  motion  (z  —  constant) . 


RIGID-BODY  MOTION 

The  coordinate  systems  and  equations  generally  used  in  the  analysis  of  the  motion  of  a  rigid  body  are  presented 
in  the  Sections  that  follow.  The  preceding  Sections  have  considered  the  motion  of  a  point  mass  in  several  co¬ 
ordinate  systems  where  position  coordinates,  as  functions  of  time,  are  sufficient  to  describe  the  motions.  For  the 
more  nearly  complete  caae  of  rigid-body  motion,  it  ia  necessary  to  consider  the  rotational  motion  of  the  body.  In 
most  cases,  it  is  convenient  to  refer  translation  and  rotation  to  the  body  center  of  gravity.  This  reference  center 
is  used  in  the  cases  that  follow,  unless  specifically  noted  otherwise.  The  notation  and  coordinates  used  are  con¬ 
sistent,  insofar  aa  possible,  with  the  reference  literature  and  current  usage. 

Thi»  Section  describes  various  coordinate  systems  and  gives  equations  of  motion  for  a  rigid  body  moving  with 
respect  to  a  flat,  nonrotating  reference  frame.  The  origin  of  each  coordinate  system  ia  located  at  the  vehicle  center 
of  gravity.  These  conditions  are  those  most  commonly  used  in  analysis  of  aircraft  motion. 

The  general  notation  and  terminology  of  established  aircraft  Mage  are  used  in  this  Section.  A  basic  notation  ia 
established  without  subscripts.  Various  specialised  axes  systems  and  corresponding  equations  are  denoted  by 
subscripts  added  to  the  basic  notation. 

Limiting  the  reference  frame  to  a  rectangular,  noe  rota  ting  system  alimlnnUe  consideration  of  Coriolis-type  forces 
in  the  equations  of  this  Section. 

*fa  this  report  Ik*  tan s  c— to*  «t  gravity  mA  rats  of  mss  ms  «**i  MWlugMHk.  tfco  MMUmUso  ip  wki*rirtn1  tor  oppMooSwu. 
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All  axis  systems  in  this  Section  are  right-hand  and  orthogonal. 

A  general  notation  for  the  force,  velocity,  and  inertia  terms  used  in  the  equation*  for  motion  of  a  rigid  body  are 
given  below.  These  items  refer  to  a  rectangular-coordinate  system  having  axes  designated  by  X,  Y,  and  Z, 
respectively.  The  origin  of  the  coordinate  system  is  at  the  center  of  gravity  Qf  the  vehicle.  The  symbols  below  are  used 
as  listed  for  vehicle  body  axes  and  with  subscripts  for  special  axis  systems. 

NOTATION 

oithonormal  b*»e  (unit*  vectors  along  X,  Y,  and  Z  coordinate  axes,  respectively 

external  force  vector  applied  at  vehicle  center  of  gravity;  includes  aerodynamic,  thrust,  and  gravity  forces 
external  force  vector  components  along  coordinate  axes,  i.e.,  F  =  Fxi  +  Ft]  +  Fik 
external  moment  vector  applied  at  vehicle  center  of  gravity 

external  moment  vector  components  along  coordinate  «.x«s,  i.e.»  G  =  Gxl  4-  Jtj  +  G*k 
total  velocity  vector  of  vehicle  center  of  gravity  (translation  of  origin  with  respect  to  a  remote  fixed  point) 
total  velocity  vector  components  along  coordinate  axes,  i.e.,  V  =  Ul  +  Vj  +  Wk 
total  angular-velocity  vector  of  vehicle  about  its  center  of  gravity 

angular-velocity  vector  components  along  coordinate  axes,  Le.,  m  ss  Pi  +  Qj  +  Rk  (Note:  P  is  the  angular  velocity  of 
rotation  about  the  X-axis  according  to  the  right-hand  rale  of  vector  representation  for  moments  and  angular  velocity.) 

gravity  or  weight  vector  of  the  vehicle 

mass  moments  of  inertia  of  the  body  about  the  X,  Y,  and  Z  coordinate  axes,  respectively 
mass  products  of  inertia  of  the  body  with  reference  to  the  X,  Y,  and  Z  coordinate  exes,  respectively 

Subscripts  used  with  the  above  symbols  denote  axis  systems  as  follows: 

Subscript  Coordinate  Axes  SyaSem 

e  earth  axes 

a  stability  axes 

p  principal  axes 

w  wind  axes 

wt  wind-tunnel  axes 

Various  axis  systems  used  frequently  in  the  analysis  of  vehicle  motion  are  described  and  sketched  in  the  pages  that 
follow.  All  of  the  coordinate  systems  presented  are  right-hand  orthogonal  systems  with  the  origin  located  at  the  vehicle 
center  of  gravity. 

Each  axia  system  is  defined  and  illustrated  in  vehicle  notation,  and  terminology  is  outlined  in  reference  3.  Special 
notation  is  defined  aa  required  and  equations  of  motion  are  listed  for  the  axis  systems  commonly  used  in  stability 
and  control  analyses. 

EARTH  AXES 

Earth  axea  ace  used  primarily  as  s  reference  system  for  the  gravity  vector,  altitude,  borisontal  distance,  and  vehicle 
orientation.  Fixed  earth  axes  provide  a  reference  for  reckoning  the  flight  path,  altitude,  and  borisontal  distance.  Earth 
axes  moving  with  the  aircraft  are  sufficient  to  define  the  gravity  vector  and  orientation  of  the  vehicle.  Both  fixed  and 
moving  axes  are  illustrated  below  with  the  preferred  sequence  of  rotation  to  define  the  orientation  angles.  (See 
reference  4). 

■> 

Description  of  Coordinate  System 

Origin  location:  arbitrary  for  fixed  earth  axes.  The  origin  of  moving  earth  axes  is  usually  placed  at  the  vehicle 
center  of  gravity. 
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Zt-axis.  along  the  gravity  vector  mg.  Thus  Z,  is  positive  downward,  i.e.,  toward  the  center  of  the  earth.  I 

Orientation  of  the  X-axis:  may  be  fixed  arbitrarily.  Fixed  earth  axes  often  have  Xe  directed  toward  the  North  | 
polar  axis.  In  moving  earth  axes  the  Xe-axis  may  be  directed  along  the  vehicle’s  initial  azimuth  heading.  f 

Ye-axis:  oriented  to  form  a  right-hand  orthogonal  axes  system.  | 

\ 


x  ,  X2 


riOUB*  4  BEL  ATJON8HIP8  BETWEEN  MOVING  AND  FIXED  BABTH  AXES.  AND  VEHICLE  BODY  AXES 
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The  axes  in  figure  4  are  defined  as  follows: 

X'„  Y'„  T,  fixed  earth  axes 

X„  Y„  4  moving  earth  axes  parallel  to  fixed  earth  axes 

X|,  Y|,  Z,  intermediate  axes  used  in  defining  orientation  of  vehicle 

X^,  Y  j,  Z_.  intermediate  axes 

X,  Y,  Z  vehicle  body  axes 


t 


The  sequence  of  rotations  defining  the  orientation  angles  of  the  body  axes  with  respect  to  moving  earth  axes  is  as  follows: 
Rotate  moving  earth  axes  Xr,  Y„  4  through  azimuth  angle  ■I'  about  4-axis  to  intermediate  axes  Xt,  Yi,  Zi. 
Rotate  axes  Xi,Y|,Zi  through  elevation  angle  6  about  Y,-axis  to  intermediate  axes  X2,  Y2,  Z-_-. 

Rotate  axes  X2,  Y2,  Z*  through  bank  angle  +  about  X2-axis  to  vehicle  body  axes  X,  Y,  and  4 


With  the  above  rotation  sequence  the  body- exes  orientation  angles  may  be  defined  in  the  following  terms: 

■*  Azimuth  or  yaw  angle  of  body  axes  from  reference  direction  of  earth  axes. 

9  elevation  or  pitch  angle  of  body  X-axis  from  the  horizontal  or  X,  Y,-plane. 

4  bank  or  roll  angle  of  the  body  Y-axis  about  the  body  X-axis  from  the  X,  Yr-plane. 

Note:  The  angles  4  and  9  are  not  necessarily  the  same  as  the  flighWpath  heading  and  the  Right-path  angle,  respectively. 


BODY  AXES 


The  body  axis  system  is  the  moot  general  kind  of  axis  system  in  which  the  axes  are  fixed  to  a  rigid  body.  The 
use  of  axes  fixed  to  the  vehicle  insures  that  the  inertia  terms  in  the  equations  of  motion  are  constant  and  thet 
aerodynamic  forces  and  momenta  depend  only  upon  the  relative-velocity  orientation  angles  «  and  p.  The  orienta¬ 
tion  of  body  axis  with  respect  to  earth  axis  is  defined  in  the  preceding  pan^aph. 

The  general  body  axis  system  is  denned  and  illustrated  below.  Special  body  axis  systems,  namely,  the  stability 
axis  system  and  the  principal  axis  system,  are  given  on  pages  13  and  14,  respectively,  of  this  Section. 


Description  of  Coordinate  System 

Origin :  vehicle  center  of  gravity. 

Reference  plane:  X4  usually  a  plane  of  symmetry. 

Positive  X-axis:  forward  along  a  reference  line  fixed  to  the  vehicle. 
Poaitrv  «•  Z-axis:  toward  bottom  of  vehicle. 

Positive  rotation:  about  Y-axis  from  Z  to  X,  i.e.,  right-hand  system. 
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TOTAL  VELOCITY  V 


FIGURE  8  VEHICLE  BODY  AXES 


The  angles  a  and  P  in  figure  5  define  the  orientation  of  the  velocity  vector  V  with  respect  to  the  body  axes  X,Y, 
and  Z.  The  angle  of  attack  o  and  the  angle  of  sideslip  fi  are  shown  in  the  preferred  yaw-pitch  rotation 
sequence.  (See  page  41.) 

Complete  equations  of  motion  referred  to  body  axes  are  given  below.  The  general  notation  defined  on  page  9  is 
used.  These  equations  are  applicable  to  any  rigid  body,  since  there  are  no  simplifying  conditions  of  symmetry 
used. 

F,  =:tn  (U-RV  +  QW) 

F,  =m(V  — PW+  RU)  -  (8) 

F,  ~m  (W  — QU-j-PV) 

c*  =PIx  --  Q Ixy  —  R I xb  —  Q R  (I*  -  It)  -PQlxr.  -  (Q2  -  R2)  1tZ  +  PRIXt 

Gt  = +QIy-RItb-P1xt  +  PR(Ix-1z) -QRIxy  -  (R2-P2)  Ibx+ QPI.t  .  (9) 

G*  ”  -f  ft  Ie—  £  Ixb  —  0  Irz  ~  PQ  (It  —  Ix)  —  R  P  Irz  —  (P2  —  Q2)  Inr  -f  R  Q  lex 

Note*:  1.  In  most  instances  a  vehicle  has  a  plane  of  symmetry,  the  XZ-plane.  The  pioduct-of-inertia  terms  Ixr  and  It*  are  aero  with  this 
symmetry,  and  the  equations  may  be  simplified  accordingly. 

2.  Gyroscopic  terms  resuiting  from  rotating  asses  in  the  vehicle  are  not  included. 
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STABILITY  AXES 


Stability  axes  are  specialized  body  axes  (see  preceding  paragraph)  in  which  the  orientation  of  the  “body  axes”  is 
determined  by  the  initial  flight  condition.  The  X„-axis  is  selected  to  be  coincident  with  the  velocity  vector  V0  at  the 
start  of  the  motion.  Consequently,  the  moment -of-inertia  and  product-of-inertia  terms  vary  for  each  initial  flight  condi¬ 
tion.  However,  they  are  then  constant  in  the  equations  of  motion. 

The  use  of  stability  axes  is  limited  to  symmetric  initial  flight  conditions  and  small-disturbance  motions 


Description  of  Coordinate  System 


Origin:  vehicle  center  of  gravity. 

Reference  plane:  X»Z.,  a  plane  of  symmetry. 

Positive  X,-axk:  coincident  with  velocity  vector  at  start  of  motion. 

Positive  Z,-axis:  toward  bottom  of  vehicle. 

Positive  rotation :  about  Y,-axis  from  Z,  to  X„  i.e^  right-hand  system. 

VKHICLB  REFERENCE  LINK,  OS 
BODY  X-JU£ 16 


FIGURE  •  STABILITY  AXES 


The  initial  angle  of  attack  a*,  is  the  angle  between  the  body  X-axis  and  the  steady  relative  velocity  vector  V0  at  the 
start  of  motion. 


Equations  of  Motion 

The  equations  of  motion  referred  to  the  stability  axes  of  a  vehicle  symmetric  about  the  XZ-plane  are  given  below. 
Symbols  are  aa  defined  on  page  9. 
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(10) 


F,<  =  m(U.  +  Q.W.-R.V.) 

F.r<  ”  m  (^»  -j-  R,  U,  —  P«  W,) 

F, .  =  m  (W.  —  Q»  U,  -f  P„  V,) 

G, ,  =  P,  lx.  -  R.  Ix*.  -  Q.  R.  (Iy,  - 1*.)  -  P.  Q>  hr.. 

Gt,  =  Q.  It.  -  R.  P„  ( I*,  llx.)  ~  (R.2  -  P.2)  hr. .  (11) 

G.,  =  +  R*  I*.  —  P,  Ixa.  —  P«  Q«  (lx.  —  Ir.)  +  Q.  R»  Ixx. 

PRINCIPAL  AXES 

A  special  set  of  body  axes  (see  preceding  paragraph)  aligned  with  the  principal  axes  of  the  vehicle  and  therefore  called 
principal  axes  is  used  for  certain  applications.  The  convenience  of  principal  axes  results  from  the  fact  that  all  of  the 
prod uct-of -inertia  terms  are  reduced  to  zero.  The  equations  of  motion  are  thus  greatly  simplified. 

Description  of  Coordinate  System 

Origin:  vehicle  center  of  gravity. 

Reference  plane:  XpZp,  a  plane  of  symmetry. 

Positive  Xp-axis:  forward  along  principal  axis  nearest  the  direction  of  motion. 

Positive  Zp*axis:  in  plane  of  symmetry,  toward  bottom  of  vehicle,  normal  to  Xp. 

Positive  rotation  about  Yp-axis:  from  Zp  to  X,„  i.e.,  right-hand  system. 

The  angle  <  denotes  the  angle  between  the  principal  axis  Xp  and  the  body  X-axis. 

6  6  r  r  I  VEHICLE  KBFRRKNCE 


FIGURE  T  PRINCIPAL  AXES 
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The  equations  of  motion  referred  to  the  principal  axes  are  listed  below.  Symbols  are  defined  on  page  9. 

FSp  =  m  (U,,  —  Rp  Vp  -j-  QP  Wp) 

F,p=m(Vp-PpWp+RpUp)  ■ 

F„p  =m(W„-QpUp  +  PpVp) 

GIp  =  PPIXp-QPRp(Irp-V 
Gyp  =QPItp—  R#Pp  (!»„  — Ixp)  ■ 

Gip  =  Rp  1%  —  Pp  Qp  (Ixp  —  Itp) 


(12) 


(13) 


GENERAL  WIND  AXES 


General  wind  axes  use  the  reticle  translational  velocity  as  the  reference  for  the  axis  system.  Wind  m  are  than  oriented 
with  respect  to  the  flight  path  of  the  vehicle,  i.e.,  with  respect  to  the  relative  wind. 

The  relation  between  general  wind  axes  and  vehicle  body  axes  defines  the  angle  of  attack  «  and  the  sideslip  angle  jg. 
These  angles  are  convenient  independent  variables  for  use  in  the  expression  of  aerodynamic  force  and  moment 
coefficients. 

Wind  axes  are  not  generally  need  in  the  analysis  of  die  motion  of  a  rigid  body,  becaase,  an  ia  the  case  of  earth 
axes,  the  moment-of-inertia  and  prodact-of-inertia  terms  in  the  three  rotational  equations  of  motioa  vary  with  tune, 
angle  of  attack,  and  sideslip  angle. 


The  general  wind-axis  system  is  defined  and  illustrated  below.  A  special  case  of  symmetric  wind  axes  follows 
on  page  16. 
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Description  of  Coordinate  System 

Origu;  vehicle  center  of  gravity. 

Reference  plan:  X«Z*flme. 

Positive  X«*axfe:  along  die  velocity  vector  V. 

Positive  Zr«ib:  in  the  vehicle  plane  of  reference  XZ  and  toward  the  bottom  of  the  vehicle. 
Positive  rotation  aboatYw-axis:  from  Z.  to  Xw,  i.e^  right-hand  system. 


The  angles  a  and  ft  ia  figare  8  me  shown  in  the  yaw-pitch  rotation  sequence.  The  orientation  relations  between 
the  body  axes  and  the  general  wind  axes  are  given  in  Section  3. 

The  equations  of  motion  of  a  rigid  body  referred  to  general  wind  axes  are  identical  in  form  to  die  equations  of  mo¬ 
tion  referred  to  body  axes.  Unis  the  equations  in  general  wind  axes  may  be  obtained  from  the  equations  given  on 
page  12.  The  mom oat-of- inertia,  and  product-of-inertia  terms  become  very  complex  in  the  general  wind  axes  system 
and  than  practically  preclude  the  use  of  these  axes  in  the  analysis  of  vehicle  motion. 


SYMMETRIC  WIND  AXES 

Hie  symmetric  case  of  the  preceding  general  wind  axes  may  be  usefully  applied  in  the  analysis  of  symmetric  vehicle 
motion,  &£.,  dive  recovery.  Symmetric  wind  axes  are  obtained  from  the  general  wind  axes  when  the  sideslip  sng)*  fi 
is  aero.  Thus  die  preceding  description  and  illustration  for  the  general  case  may  be  used  directly  with  ft  =  0. 


The  equations  of  motion  for  the  symmetric,  unbanked  flight  of  a  vehicle  with  a  plane  of  symmetry  are  given  below. 


Fxw  =mUf 

Ftw  =  0 

Fa.=—  mQ.U, 
Cx,  0 
Gtw  =  Q»  Itw 
Ga*  =  0 


(14) 


(15) 


WIND-TUNNEL  STABILITY  AXES 

Wind-tunnel  stability  axes  are  wed  aa  a  reference  system  for  measuring  and  reducing  aerodynamic  data  in  wind-tunnel 
testa.  This  set  of  axes  diflers  from  Ihe  previous  stability  axes  in  that  the  Zwt-axis  is  aligned  normal  to,  and  remains 
normal  to,  the  relative  wind,  whereas  the  general  stability  axes  are  body  axes  determined  by  the  initial  flight  condition. 

Since  it  ia  not  convenient  to  une  wind-tunnel  stability  axes  in  analysis  of  the  motion  of  a  vehicle,  the  equations  of 
motion  are  omitted  for  this  case. 

Description  of  Coordinate  System 

Origin  location:  in  the  reference  plane  of  the  vehicle  at  the  point  corresponding  to  the  vehicle  center  of  gravity. 
Reference  plane:  dm  X„tZ»t  plena. 

Positive  Zvt-exb:  in  the  reference  plane  of  the  vehicle  perpendicular  to  the  relative  wind  V. 

Positive  Xwrttris:  toward  the  forward  part  of  the  vehicle,  along  the  projection  of  the  relative  wind  V  upon  dm 

Positive  Y^vxh:  oriented  to  form  i  right-hand  orthogonal  axes  iphn 


M 


VEHICLE  BSFttBBNGE  LIMB, OB 
BOOT  X-AXIS 


Y,  Twt 


FIOUBB  •  WKD-PJ10IBL  ST  ABILITY  AXES 


la  the  above  figure  die  angle  of  attack  a  and  the  sideslip  angle  fi  give  tbs  orientation  of  the  relative-velocity 
vector  to  the  vehicle  body  axes.  Wind-tonnel  stability  axes  sse  the  yaw- pitch  rotation  sequence  (page  46). 


NONROLLING  AXES 


The  problem  of  formulating  equations  of  motion  for  a  symmetric  rolling  body  may  be  simplified  by  using  a  nonrolling 
axis  system.  Nonrolling  axes  are  a  special  set  of  body  axes  having  the  Y-axis  always  horisontal  and  the  XZ-piaoe 
always  vertical  This  axis  system  may  be  ased  even  though  the  body  rotates  aboat  the  X-axis.  It  is  necessary,  however, 
that  the  inertia  parameters  and  the  aerodynamic  forces,  moments,  and  derivatives  he  constant  with  respect  to  the 
aonreRiag  reference  frame.  Thus  the  body  must  have  rotational  symmetry  shoot  the  X-axis.  Applications  of  nonrolling 
axes  to  the  motion  analyses  of  aircraft,  projectiles,  and  missiles  are  given  in  reference  5. 
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SECTION  3.  COORDINATE -SYSTEM  TRANSFORMATIONS 

In  vehicle-motion  analysis  it  is  frequently  necessary  or  expedient  to  transform  coordinates,  vector  components,  inertia 
parameters,  a^d  stability  derivatives  from  one  coordinate  system  to  another.  The  following  section  gives  the  relations 
most  frequently  used  in  such  transformations. 

Equations  are  used  to  express  the  transforming  relations  whenever  these  relations  are  simple  and  not  often  repeated. 
Matrix  notation  is  used  in  the  more  complex  transformations,  and  a  tabular  presentation  is  given  when  the  forms  of  a 
transformation  relation  are  similar  for  several  cases. 

PARTICLE-MOTION  TRANSFORMATIONS 

The  transformations  between  the  coordinate  systems  useful  in  the  analysis  of  particle  motion  are  given  in  the  pages 
that  follow.  Notation  and  definitions  of  terms  are  consistent  with  those  used  in  the  preceding  sections. 

RECTANGULAR  COORDINATES 

Cartesian  or  rectangular  coordinates  are  perhaps  the  most  commonly  used  coordinates.  The  following  pages  give 
relations  for  translation  of  the  origin  and  the  rotation  of  rectangular-coordinate  systems  about  the  origin.  Composite 
flungM  involving  translation  and  rotation  of  the  coordinate  axes  may  be  accomplished  by  successive  application  of 
these  two  basic  transformations. 

The  equations  relating  spherical  coordinates  to  rectangular  coordinates  and  the  equations  relating  cylindrical  coordinates 
to  rectangular  coordinates  are  also  included  in  this  section. 

Symbols  and  notation  are  defined  when  first  used  or  as  required. 


Translation  of  the  origin  in  rectangular  coordinates  is  illustrated  in  the  figure  below. 


TRAN8L.ATRD 
z*  REFERENCE 


It 


NOTATION 


0  origin  of  rectangular  coordinate  system  with  axes  X,  Y,  and  Z 
x,  y,  z  position  coordinates  of  point  P 
x,  t,  x  component*  of  the  nxtor  OP 

a,b,c  position  coordinates  of  translated  origin  O'  in  initial  coordinate  reference  system.  (These  coordinates 
are  considered  as  constants.) 

(  )'  denotes  coordinates  and  quantities  referred  to  the  translated  coordinate  system 

From  figure  10  the  relation  between  coordinates  in  the  initial  and  in  the  translated  coordinate  system  is 

x'  =  x  —  a 

/  =  y-h 

s'  —  X  —  C 

Since  there  is  no  rotation  of  the  coordinate  axes  with  a  pure  translation,  the  components  of  a  vector  at  P  referred  to 
the  axes  X,  Y,  and  Z  are  identical  to  the  components  referred  to  the  axes  X',  Y',  and  Z'.  Consequently,  components  of 
vectors  such  as  force  and  moment  vectors  are  unchanged  by  a  translation  of  the  origin.  Velocity-  and  acceleration-vector 
components  are  unchanged  also,  except  when  the  translated  axis  system  becomes  a  moving  reference  system. 


Rotation  of  rectangular-coordinate  axes  about  the  origin  is  very  often  useful  and  sometimes  quite  necessary.  A  general 
rotation  of  an  orthogonal-axes  system  may  be  accomplished  by  three  successive  planar  rotations;  hence  a  simple  planar 
rotation  is  considered  first  and  then  extended  to  the  general  case.  Also  included  in  this  section  are  the  direction-cosine 
relations  for  defining  a  general  rotation  of  rectangular-coordinate  axes.  The  relations  given  in  the  following  pages  are 
developed  in  many  standard  mathematics  and  engineering  texts,  such  as  reference?  6.  7,  and  8.  The  tabular  presentation 
of  the  transformation  relations  is  adapted  from  reference  9.  J 

1.  PLANAR  ROTATION 

The  rotation  given  below  corresponds  to  a  rotation  in  the  XY -plane  about  the  Z-axis.  The  subscript  I  denotes  the 


The  components  of  the  vector  OP  in  the  preceding  figure  are  transformed  from  x,  y,  and  %  to  x,,  y,,  and  Zi, 
respectively,  by  the  following  equations: 

xt  —  x  cos  ^  -f  y  sin 

yi  =  —x  sin  ^  +  y  cos  $ 

*1  =  z 


These  equations  may  conveniently  be  expressed  in  matrix  form. 


"X|” 

cos  \f/  sin  $  0 

~  x~ 

’  X  ’ 

yi 

= 

—  sin  $  cos  tp  0 

y 

=  M 

y 

_*i  _ 

0  0  1. 

_  z  _ 

.3. 

The  components  of  any  vector  in  the  X Y -plane  may  be  transformed  by  the  above  relations. 

Since  the  transformation  matrix  is  orthogonal,  the  inverse  transformation  [*p]  is  given  by  the  transpose  of  [^]. 
Since  the  transpose  of  a  matrix  is  obtained  by  interchanging  the  rows  and  the  columns,  in  this  case  the  inverse  trans¬ 
formation  matrix  [f]_1  is  defined  as 

T  cos^  — aii'up  0 

[^]  _1  as  [^r]'  =  sin^f  cos if/  0 

LOO  1 

whence 


X 

"  Xi  " 

y 

=  w-1 

yi 

X 

_  *i  . 

Note  that  this  procedure  is  equivalent  to  replacing  the  angle  t/'  by  ^  and  interchanging  the  subscripted  and  unsub- 
scripted  components  in  the  first  form  of  the  equations. 


It  is  convenient  to  introduce  a  tabular  presentation  for  the  transformation  matrix  and  its  inverse.  Table  1  gives 
the  transformation  matrix  array  with  initial  position  coordinates  at  the  head  of  each  column  and  the  rotated  co¬ 
ordinates  in  front  of  each  row.  From  this  array  the  transformation  and  inverse-transformation  equations  are  written 
by  using  the  matrix  elements  as  coefficients  of  the  appropriate  vector  components  in  the  transforming  equations. 


TABLE  1. 

VECTOR  TRANSFORMATION  MATRIX 
INITIAL  AXES  TO  ROTATED  AXES 


COMPONENTS  IN  INITIAL  COORDINATE  SYSTEM  ] 

X 

Y 

Z 

COMPONENTS  IN 
ROTATED  COORD. 
SYSTEM 

X, 

coa  ^ 

ain  ^ 

0 

Y, 

—sin  <[> 

coa  f 

0 

z, 

0 

0 

1 

Direct  transformation  equations  are  obtained  by  summing  horizontally  along  ~«ch  turn. 

x,  =  (coa  *)  x  -|-  (sin  *)  y  +  (0)  t 

y,  =  ( —  rin x  +  (coa^f)  y  +  (0)  z 
a,  =r(0)x-M0)y+(l)i 
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Inverse  transformation  equations  are  obtained  by  summing  vertically  down  each  column. 

x  —  (cos^)  xi  +  (~sin  ifr)  yt  -f  (0)  *1 
y  =  (sin  >fi)  %i  -f  (cos  if,)  yt  +  (0)  zt 
x  =  (0)  xt  -f  (0)  yi  -j-(l)  *i 

2.  GENERAL  ROTATION 

He  general  rotation  of  a  rectangular  axes  system  may  be  accomplished  by  successive  planar  rotations  of  the  type 
described  in  the  preceding  paragraph.  In  making  a  general  rotation,  however,  the  sequence  of  rotation  is  important. 
The  basic  order  of  rotation  is  described  and  illustrated  below.  This  sequence  of  rotation  and  the  terminology  have 
been  used  extensively  in  aircraft  motion  analysis  (re Terences  4  and  7). 


cjv 


p*' 


rtaxmn  i*  cjbwbrax  botation  about  obichk—  ouextation-ancilk  n—owmow  or  botation 


The  transformation  matrices  are  given  in  tabular  form  for  the  basic  sequences  and  for  several  other  sequences  of  axis 
rotation.  These  transformations  may  be  applied  to  position,  velocity,  force,  moment,  and  acceleration  vectors  to  obtain 
their  components  in  the  rotated-axes  system.  Both  the  direct*  and  inverse-transformation  relations  are  obtained  from 
the  tabular  presentation  as  shown  in  the  sample  problem. 

In  the  preceding  sketch  the  coordinate  axes  are  designated  by  capital  letters  (X,  Y,  Z)  and  the  position  coordinates  by 
lower  case  letters  (x,  y,  a).  The  Greek  symbols  fa  8,  and  \j>  are  used  to  refer  to  angular  rotation  about  the  X-,  Y-, 
and  Z-axes,  respectively.  Subscripts  refer  to  the  various  rotated-axes  systems.  Thus  the  subscript  3  denotes  the  final 
axes  and  coordinates.  Similar  subscripts  are  used  with  the  rotation  angles  to  indicate  the  reference  axes  for  the 
particular  angle.  The  basic  order  of  rotation  is 

1.  “Yaw”  about  Z-axis  through  the  angle  fa 

2.  “Pitch”  about  Yi-axis  through  the  angle  8\. 

3.  “Roll”  about  X^-axis  through  the  angle  fa. 

General  transformations  for  rotation  of  rectangular-coordinate  systems  are  tabulated  in  table  2.  Both  direct  and  in¬ 
verse  transformations  are  given,  as  illustrated  in  the  sample  problems  on  this  page.  Also  included  in  this  table 
are  the  equations  for  the  instantaneous  angular  velocities  fa,  8$,  and  fa  about  the  final  coordinate  axes  in  terms 
of  the  orientation  angles  and  their  rates  of  change. 

The  first  esse  lifted  is  the  most  commonly  used  order  of  rotation.  Cases  1  and  2  may  be  considered  as  fundamental 
rotations.  The  remaining  cases  may  be  obtained  from  cyclic  permutations  of  the  initial  coordinate  and  angle  notation, 
it  sf  be  noted  that  changing  the  sequence  of  rotation  changes  the  definition  of  the  orientation  angles.  Consequently, 
arv?,k'f  Kf'ith  different  subscripts  are  not  interchangeable,  i.e.,  generally  fa  =/=  fa  fa.  Also,  the  orientation  angle 
rat»^  >j’  change,  i.e.,  fa  8U  fa,  are  not  orthogonal. 


Use  of  table  2  is  illustrated  by  the  sample  problems  below. 


Example  1.  Direct  Transformation 

Given:  Velocity-vector  components  Vx,  Vv,  Vr.. 

Rotation  order:  yaw,  pitch,  roll  to  0t  to  fa,  as  in  Case  1  of  table  2) . 

find:  Velocity-vector  component  along  Zj-axis  (VZa) . 

Solution:  Write  equation  for  V*s  by  summing  terms  along  the  Z8-row  of  the  vector  transformation  matrix. 

Vj-S  =  (cos  fa  sin  0i  cos  >/r  4-  sin  fa  sin  f)  Vx 
-f-  (cos  fa  sin  0i  sin  f  —  sin  fa  cos  Vr 
4-  (cos  fa  cos  0i )  Vz 

Example  2.  la  verse  Transformation 

Given:  Acceleration  vector  components  aXj,  at3,  aZj  along  final  coordinate  axes. 

Rotation  order:  yaw,  pitch,  roll  to  0t  to  fa,  as  in  Case  1  of  table  2) 

Fird:  Acceleration  vector  component  along  X-axis  (ax). 

Solution:  Write  equation  for  ax  by  summing  terms  down  the  x-column  of  the  vector  transformation  matrix. 
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TABLE  2 

RECTANGULAR  COORDINATE  TRANSFORMATIONS 
ROTATION  OF  AXES  ABpUT  ORIGIN 


CASE  4  zu  Z,  Z  INITIAL  VECTOR  COMPONENTS 


initial  vector  components 
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ax  =  (cos  61  cos  \p)  ax3 


4-  (sin  <p->  sin  6\  cos  i p  —  cos  <f>>  sin  ip)  ay3 
-f-  (cos  <f> 2  sin  $i  cos  i p  -}-  sin  fa  sin  \p)  a%3 

A  general  rotation  transformatic -  of  vector  components  from  one  coordinate-axes  system  to  another  may  be  interpreted 
in  terms  of  direction  cosines.  The  direction  cosines  are  defined  as  the  cosines  of  the  angles  between  each  of  the  final 
coordinate  axes  X3,  Y3,  and  Z3  and.  each  of  the  original  coordinate  axes  X,  Y,  and  Z.  Thus  there  are  nine  angle* 
required  to  describe  a  general  rotation  of  rectangular-coordinate  axes.  Direction  angles  that  locate  the  X3-axis  with 
respect  to  the  original  X-,  Y-,  and  Z-axes  are  illustrated  in  the  figure  below.  Similarly,  direction  angles  are  defined 
for  the  Y3-  and  Z3-axes. 


K 

FIGURl,  18  GENERAL  ROTATION  ABOUT  ORIGIN  -  DIRECTION-COSINE  DESCRIPTION  OF  ROTATION 

(DIRECTION  ANGLES  FOR  X#  AXIS) 


The  direction-angle  notation  used  is  as  follows: 


(X,  X3)  =  Angle  between  X-  and  X3-axes. 

(Y,  X;i)  —  Angle  between  Y-  and  X*-axes. 

(Z,  X3)  =  Angle  between  Z-  and  X3-axes,  etc. 


The  cosines  of  the  direction  angles  may  be  arranged  as  a  vector  transformation  matrix  and  ua*d,  exactly  as  in  the 
preceding  Section,  to  transform  vector  components.  The  vector  transformation  matrix  of  direction  cosines  is  shown 


in  table  3. 


TABLE  8 

VECTOR  TRANSFORMATION  MATRIX  OF  DIRECTION  COSINES 


_ — 

INITIAL  VECTOR  COMPONENTS 

X 

Y 

Z 

TRANSFORMED 
VECTOR  COMP. 

X, 

cos  (X,  X.) 

cos  (Y,  Xj) 

cos  <Z,  Xj) 

V, 

cos  (X,  Y>) 

cos  (Y,  Y>) 

cos  (Z.  Yj) 

z. 

cos  ( X,  Z») 

cos  (Y,  Z.i) 

cos  (Z,  Za) 
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Relations  (or  the  direction  cosines  in  terms  of  the  orientation  angles  4,  #,  sad  f  may  be  '***h»««*  by  equating 
corresponding  elements  of  the  above  matrix  and  the  appropriate  matrix  of  table  2.  For 

cos  ( Y,  Zs)  =  cos  sin  0i  sin  ^  —  sin  cos  ^ 

for  axis-orientation  angles  defined  by  the  rotation  sequence  of  Case  1  and 

cos  (Y,  Zs)  =  —  sin  fa 

when  the  orientation  angles  are  defined  as  in  Case  3. 

3.  RECTANGULAR  COORDINATES  TO  CURVILINEAR  COORDINATES 

The  transformation  of  rectangular  space  coordinates  to  a  curvilinear-coordinate  system  involves  a  nonlinear  co¬ 
ordinate  change.  The  relations  used  to  change  from  rectangular  coordinates  to  spherical  coordinates  are  given  as 
equation  16  and  those  used  for  the  transformation  to  cylindrical  coordinates  are  given  as  equation  17.  In  both 
cases  it  can  be  seen  (figures  14  mid  15)  that  the  transformation  equations  are  statements  of  simple  trigonometric 
relationships. 

I 


Rectangular  to  spherical  coordinates: 


Spherical  to  rectangular  coordinates: 
X  =  r  sin  4>  cos  t 
y  =  r  sin  4  sin  0 
m  —  rco*4 
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Bcdflophk'  to  cylindrical  coordinates: 


/ 


r  =  V**  +  Xs 
x 

X  =x 

Cylindrical  to  rectangular  coordinates: 

x  —  r  cos  6 
y  =  r  tin  0 

x  —  x 


s 


(17) 


It  is  important  to  note  that  base  vectors  at  each  point  in  the  preceding  curvilinear-coordinate  systems  are  defined  as 
orthogonal.  Consequently,  at  a  given  point,  transformation  of  vector  components  from  rectangular  to  spherical  (or 
cylindrical)  coordinates  is  s  transformation  between  rectangular  axes  and  correspond*  to  a  rotation  of  the  axes  system 
at  the  point.  Hies  appropriate  transformation  matrices  from  table  2  may  be  uwed  directly  to  transform  vector  com¬ 
ponents  from  rectangular  to  spherical  (or  cylindrical)  coordinates  at  a  point 

For  example,  at  a  point  the  vector 

V  =  Vxi  +  VTJ  +  V*k 

—  \&  +  V*4  4-  VX 

may  he  changed  from  rectangular  components  (Vx,  VT,  V*)  to  spherical  components  (V„  V*,  V«)  by  a  linear 
transformation  corresponding  to  oae  of  the  transformation  matrices  given  in  table  2.  This  method  ia  used  in  the 
following  paragraph  to  transform  from  spherical  to  flight-path  coordinates. 

SPHERICAL  COORDINATES 

The  relations  used  to  change  from  rectangular  to  spherical  coordinates  are  given  in  equation  16.  As  is  noted 
there,  this  ia  a  nonlinear  transformation.  However,  at  any  point  in  a  apace  described  by  spherical  coordinates,  a 
rotation  of  local  base  vectors  ia  accomplished  by  a  rotation  of  rectangular  axes. 

In  this  Section,  presentation  of  vector  transformations  is  limited  to  the  change  from  local  spherical-coordinate 
axes  to  flight-path  axes.  These  axes  are  defined  and  illostmted  in  Section  2.  This  transformation  serves  to  further 
illnatrate  the  use  of  rectaagalsrcoordinatc  transformations  (table  2). 
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Hie  rotation  of  the  local  base  vectors  at  point  P  from  a  spherical-coordinate  orientation  to  flight-path  axes  is  de¬ 
fined  in  Section  2  on  page  5.  It  is  noted  that  this  rotation  correaponds  to  die  roll-pitch-yaw  sequence  (Case  6  of 
table  2)  of  the  rectangular-coordinate  transformations.  Hins  by  identifying  quantities  of  the  transformation  (Case 
6,  table  2)  with  the  notation  defined  for  spherical  and  flight-path  coordinates  in  Section  2  the  desired  transforma¬ 
tion  matrix  is  obtained.  The  correspondence  between  terms  is  given  in  table  4. 

TABLE  4 


CORRESPONDENCE  BETWEEN  RECTANGULAR-COORDINATE  TRANSFORMATION 
AND  SPHERICAL  -  FLICHT-PATH-COORDINATE  TRANSFORMATION 


Item  in  Rectangular-Coordinate  Transformation 
(Case  t».  Table  2) 

Corresponding  item  in  Spherical  —  Flight-Path-Coordinale  System 
(Section  2) 

Initial  Vector 
Component* 

X 

G, 

Vector  Components 

Along  Spherical- 
Coord  mate  Directions 

Y 

G* 

Z 

G. 

Transformed 

Vector 

Components 

X, 

C. 

Vector  Components 

Along  Flight-Pat) 

Coordinate  Direction* 

Y, 

G. 

Z. 

G» 

Orientation 

Angie* 

♦ 

(904) 

Flight-Path  Orientation 

Angies 

#i 

y 

f. 

n 

A  general  vector  Gmay  be  expressed  as  follows: 

G  =  G4r  4*  Gpfy  +G4«  (spherical  coordinates) 

G  =  Ga  -f  Ga  -f  GtCt  (flight-path  coordinates) 

Substitution  of  the  above  items  in  the  vector  transformation  matrix  of  table  2,  Case  6,  results  in  the  transformation 
shown  in  table  5. 

Note:  sin  (90  —  J)  =  cos  4  and  cos  (90  —  4)  =  sin  I 


TABLE  5 

VECTOR  TRANSFORMATION  MATRIX 
SPHERICAL  COORDINATES  TO  FLIGHT-PATH  COORDINATES 


COMPONENTS  ALONG  SPHERICAL  COORDINATES 

G, 

c* 

G* 

COMP.  ALONG 

FLT.  PATH  COORD. 

C. 

CO*  7  coo  a 

*ia  1  «ia  o 

4- co*  1  aha  7  ce*  v 

ooo*  wm  o 
— «*■  1  ma  7«oo  v 

D 

—co*  7  *ui  * 

maicost 
—ooo  I  sin  7  *in  s 

00*  1  00*  * 

-f  *m  1  *in  7  sin  * 

H 

nb  y 

—COS  loot  7 

afatl  ooo  > 

27 


Base  (met)  vectors  1„  1^,  lt  are  aligned  along  spherical-coordinate  directions. 

(■Wwrton  *■>«»,  «t  *re  aligned  along  the  flight  path  (  eT  is  along  the  velocity  vector  of  the  point  P  ) . 
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FIGURE  IS  ROTATION  OF  BASS  VBCTOBfl  FROM 

SPHERICAL  TO  FLIOHT-PATH  COORDINATES 


CYLINDRICAL  COORDINATES 

Equation*  relating  rectangular  and  cyiiurbicai  apace  coordinates  are  given  on  page  26.  This  coordinate  change  is 
nonlinear,  as  is  the  change  to  spherical  coordinates.  The  discussion  concerning  rotation  of  local  base  vectorn  for 
spherical-coordinate  systems  also  applies  to  cylindrical-coordinate  systems.  A  similar  procedure  may  be  used  to 
obtaiu  the  transformation  matrix  for  a  rotation  of  the  local  base  vectors. 
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RIGID-BODY  TRANSFORMATIONS 


Transformations  useful  in  the  analysis  of  rigid-body  motion  follow  directly  from  the  rotation  of  rectangular- coordinate 
axes.  The  translation  of  a  rigid  body  may  be  considered  as  the  motion  of  a  particle  concentrated  at  the  center  of 
gravity  of  the  body,  and  the  foregoing  Sections  may  then  be  applied.  Rotation  of  a  body  about  its  center  of  mass, 
however,  introduces  additional  degrees  of  freedom  that  depend  upon  the  inertia  characteristics  of  the  body. 

The  transformation  relations  useful  in  the  analysis  of  rigid-body  motion  are  given  in  this  Section.  First  the  general 
transformation  of  vector  components  and  inertia  parameters  is  discussed.  Then  specialised  transformations  such 
as  those  used  in  conventional  aircraft-motion  analysis  are  given. 

GENERAL  AXES  TRANSFORMATIONS 

The  transformations  of  this  Section  are  limited  to  rotations  of  rectangular-coordinate  systems.  As  is  noted  previously 
for  spherical  and  cylindrical  coordinates,  the  rotation  of  a  rectangular  axis  system  about  ita  origin  may  be  utilized 
for  local  orthogonal  axis  systems  even  though  the  origin  of  this  system  may  be  defined  with  respect  to  a  curvilinear 
coordinate  system. 

Vector  components  ore  transformed  in  the  case  of  rigid-body  motion  in  the  same  manner  as  vector  components  are 
transformed  for  particle  motion.  Thus  die  transformation  matrices  and  method  given  on  page  23  are  directly  appli¬ 
cable  to  the  rigid-body  case. 


General  transformations  of  the  inertia  parameter,  which  are  important  in  rigid-body  rotational  motion,  do  not  have 
the  simple  form  of  the  vector  transformation.  This  transformation  is  presented  for  a  general  rotation  of  the 
coordinate  axes. 


The  general  transformation  for  moment-of-inertic  and  produet-of-inertia  terms  may  be  «1  from  the  vector 

transformation  matrix  and  its  inverse.  This  procedure,  given  in  reference  10,  is  outlined  below. 


Rotational  njotion  of  a  rigid  body  is  expressed  by  die  fundamental  relation 

1  ~  [!]•  '  G8) 

where  I  is  the  angular  momentum  vector  km  «  is  the  angular  velocity  vector.  The  inertia  matrix  [i]  is  defined  as 
follows: 


m= 


lx 

—Ixx 


L  — lx* 


-IXT 
It 
— In 


-Ixx  ‘ 
— Ita 
Is  . 


In  this  matrix  Ix,  It.  and  I*  are  the  mass  momenta  of  inertia  and  Its.  Ixx.  and  Ixt  an  the  mass  prod¬ 
ucts  of  inertia  with  reference  to  the  X,  Y,  and  Z  axes,  respectively. 


The  above  vector  equation  for  rigid -body  rotation  is  independent  of  the  coordinate  system  ejected  to  represent 
the  vectors.  Honor,  if  the  subscript  o  denotes  reference  to  the  original  and  3  to  the  trsnsfonne  i  coordinate  system, 
the  equation  for  rigid-body  rotation  may  be  written 

I.  =  (UV.  «l.  =  [!.]«. 
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Vector  transformations  from  table  2  nay  then  be  need  to  change  the  vectors  I„  and  «•„  from  the  original  to  the 
new  coonfioats  system.  These  transformation  relations  may  be  expressed  as  follows: 


Is  =  [r]  lo 

*1  =  [r]  Wo  orw0  =:  [  r]  _1Wg 

where  the  vector  transformation  matrix  [r]  and  its  inverse  [r]-1  denote  a  general  rotation  of  axes  as  given 
on  page  20. 

Combining  the  transformed  vectors  with  the  body-rotation  equation  results  in  the  desired  transformation  relation 
for  the  inertia  matrix.  Tims 


I*  =  fr]  [UJ«.  =  [r]  [IJ  [r]-»w,  =  [i3] 

Hence 

tW  =  lr][Uitr]-> 


(19) 


Expansion  of  this  transformation  relation  for  a  general  rotation  defined  by  three  orientation  angles  results  in  a 
complex  expression  having  a  large  number  of  terms.  This  is  a  straightforward  procedure  using  matrix  multiplication, 
but  in  moat  cases  it  is  impractical  and  unnecessary.  Body  symmetry  usually  reduces  two  of  the  product-of  nertia 
terms  to  sero,  and  in  many  cases  a  simple  planar  rotation  is  sufficient  to  define  the  axes  rotation.  These  practice 
considerations  simplify  the  expansion  of  the  inertia  transformation  relations  used  later. 

In  order  to  illustrate  the  transformation  of  a  matrix,  the  foregoing  procedure  is  expanded  below  for  a  planar  rotation. 
This  particular  case  provides  a  general  form  that  is  subsequently  useful  in  the  transformation  of  airplane  stability 
derivatives. 


If  [A]  and  [A]  represent  the  original  and  the  transformed  matrices,  respectively,  the  matrix  transformation 
relation  is 

[A]  =  [r]  [A]  [r]-» 


"  •« 

•v 

*11 

»« 

III 

[A]  = 

•r« 

•tt 

•r* 

and  [A]  = 

■y* 

*rr 

- 

•v 

a.  _ 

_  *n 

•v 

If  the  vector  transformation  matrix  fr]  corresponds  to  a  simple  rotation  about  the  Y-axis,  the  above  equation 
may  be  written  as  follows: 


cos  8  0 

0  1 

_  sin  9  0 


—  sin  9 

0 

cos  6 


<» 

•»r 

cos 

9 

0 

sin 

9  “ 

■r* 

•rr 

0 

1 

0 

^*1 

**J 

>u  . 

—  sin 

9 

0 

cos 

9  . 

Note  that  the  matrix  may  be  obtained  from  table  2.  For  example,  using  the  vector  transformation  matrix  for 
case  1  of  this  table,  aubstitotioa  of  9  =  0,  9\  =  9  nsd  —  Q  results  in  the  above  matrix. 


Expansion  of  the  matrix  multiplication  above  gives  the  transformed  matrix.  Thus 
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Transformations  of  vectors,  inertia  parameters,  and  stability  derivatives  used  in  vehicle  notion  analysis  are  sum¬ 
marized  in  this  Section.  These  transformations  refer  to  the  vehicle  axes  systems  defined  in  Sectioa  2  and  nse 
established  aircraft  notation  and  terminology. 

It  is  convenient  to  consider  the  vehicle  transformations  in  two  groups.  The  first  group  involves  single  rotations 
about  the  lateral  (Y )  axis.  A  change  from  body  axes  to  stability  axes  is  a  single  rotation  of  this  type.  The  second 
group  comprises  cases  of  general  rotation  such  as  a  change  from  earth  axes  to  body  axes. 


The  inertia-parameter  and  stability-derivative  transformations  are  not  given  for  the  second  group.  The  general  ro- 
tation  cases  are  used  principally  to  define  orientation  of  axes  fixed  on  the  vehicle,  with  respect  to  the  *arth  or 
the  relative  wind.  As  is  noted  previously  (page  16),  the  i-.r.  ^.  parameters  become  functions  of  the  orientation 
angles  and  the  analysis  of  the  motion  is  then  unnecessarily  complicated.  The  transformation  of  the  inertia  matrix 
for  a  multiple  rotation  may  be  developed  horn  the  general  relationship  given  by  equation  19. 

I.  SINCLE  ROTATION 

The  pitch  rotations  used  to  change  from  stability  or  principal  axe s  to  body  axes,  and  vice  versa,  are  illustrated  in 
figure  17.  These  axes  and  the  notation  are  defined  in  Section  2.  Both  stability  axes  and  principal  axes  are  fixed 
to  the  vehicle  and  are  therefore  simply  special  ‘body*  axes.  The  angles  between  these  various  axes  systems  are 
measured  as  rotations  about  the  Y-axis. 
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VEHICLE  BODY  AXES 
STABILITY  AXES 

PRINCIPAL  AXES 


FIGURE  IT  BODY  -  STABILITY  -  PRINCIPAL  AXES  ROTATION  —  PITCH 


TABLE  6 

VECTOR-TRANSFOBMATION  MATRICES 
BODY-STABILITY-PRINCIPAL  AXES 


STABILITY  AXES 

TO  BODY  AXES 

VECTOR  COMPONENTS  STABILITY  AXES  | 

X. 

V. 

*■ 

VECTOR 
COMPONENTS 
BODY  AXES 

X 

C06  am 

0 

—sin  a*  ? 

-  ..  .  ...  ...  * 

Y 

0 

1 

0 

Z 

tin  <*. 

[  0 

cos  a  o 

PRINCIPAL  AXES 

TO  BODY  AXES 

VECTOR  COMPONENTS  -  PRINCIPAL  AXES  j 

X, 

Y, 

Zp 

VECTOR 
COMPONENTS 
BODY  AXES 

X 

cot  t 

0 

—sin  • 

Y 

0 

1 

0 

Z 

tin  « 

i 

CP*  t 

STABILITY  AXES 

TO  PRINCIPAL  AXES 

VECTOR  COMPONENTS -STAB 

I  CITY  AXES 

X. 

Y. 

Z. 

VECTOR 
COMPONENTS 
PRINCIPAL  AXES 

X, 

cot  * 

0 

—sin  t| 

Y. 

0 

1 

r\ 

Z, 

tin  t 

0 

CO*  1) 

DIRECT  TRANSf ORMATION  —  Sum  horiaootally  ilong  each  row. 
INDIRECT  TRANSFORMATION  —  Sum  vertically  down  each  column. 
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(Vector  transformation  matrices  may  be  obtained  directly  from  those  in  table  2  by  letting  ^-=0,  ^  =  0,  and 
0 s  =  a«i,  e,  or  rj.  These  vector  transformation  matrices  are  given  in  table  6.  Examples  of  the  use  of  *hese  matrices 
may  be  found  on  page  20. 

Transformation  of  the  inertia  matrix  for  a  general  rotation  is  discussed  on  pages  29—31.  However,  for  the  vehicle 
axis  systems  used  in  this  Section  the  XZ-plane  is  a  plane  of  symmetry*  and  only  rotations  about  the  Y-axis  are 
considered.  These  conditions  greatly  simplify  the  inertia  matrix  and  its  transformation. 

t 

j  The  inertia  matrices  of  an  aircraft  referred  to  body,  stability,  and  principal  axes  are  given  below. 


\ 


i 

" 

Iv 

0 

-Ixz  “ 

| 

Body  Axes 

[I]  = 

0 

It 

0 

(21) 

j 

- 

— l\Z 

0 

Is  . 

j 

A 

< 

- 

Ixh 

—  I.XXg 

(22) 

| 

Stability  Axes 

[Is]  = 

0 

Iv8 

0 

3 

l 

- 

—  IX7„ 

0 

I*K 

i 

■ 

lx,. 

0 

0  " 

i 

Principal  Axes** 

[Ip]- 

0 

«v. 

0 

(23) 

_ 

0 

0 

Izr 

The  transformation  relations  for  tl  e  elements  of  the  r.bove  matncec  are  given  in  table  7.  These  equations  are  ob¬ 
tained  by  identifying  elements  ot  the  inertia  matrices  with  the  corresponding  olements  of  the  general  matrix  trans- 
fonnr*ioK  on  pages  30  and  31.  The  appropriate  angle  substitution  may  be  determined  from  figure  17. 

»  Some  additional  relations  pertinent  to  inertia-parameter  transformations  are  listed  c  *low. 

1.  — 

2.un2,  =  ^lH- 

J  *Z  —  «x 


3.  tan  2«j  ~ 


lx8  -- 


Stability  derivatives  are  used  extensively  n  the  analysis  of  aircraft  motion.  They  are  introduced  with  the  linearization 
of  the  aerodynamic  foroe  a. ad  moment  relations.  These  derivatives  may  be  conveniently  arranged  in  matrix  form; 
hence  it  is  frequently  necessary  to  tr an.. form  them  from  one  axes  system  to  another.  Transformations  relating  stability 
derivatives  in  terms  of  body  and  stability  axes  are  given  in  this  Section.  This  involve*  a  simple  rotation  of  the  axes 
system  about  tire  Y-axis  and  follows  directly  r'rom  the  matrix  transformation  on  page  30. 

The  notation  for  stability  derivative*  is  confused  in  existing  literature.  It  is  therefore  appropriate  to  reiterate  here 
the  statement  on  page  129  of  reference  4  reminding  the  reader  to  exercise  extreme  care  in  using  the  literature  tha; 
involves  stability  derivatives.  This  is  necessary  to  insure  that  the  definitions  of  the  symbols  used  are  fully  understood 
and  that  comparisons  and  results  will  be  correctly  interpreted.  (Also  see  reference  2.) 

Notation  in  this  report  is  consistent  with  that  used  in  references  3  and  4.  This  notation,  defined  in  Section  4,  im¬ 
plies  diffem»fiati~."  of  direct  forces  and  moments  with  respect  to  perturbation  quantities. 

•  By  symmetry  about  the  XZ-plane,  Ixt  and  IV7  are  eero. 

**  Principal  axes  are  defined  by  the  condition  IXTr  —  lrzr  =  IxxP  =  0 


TABLE  7 

INERTIA  MATRIX  ELEMENT  TRANSFORMATIONS 

BODY  -  STABILITY  -  PRINCIPAL  AXES 
Angles  are  positive  as  shown  on  Page  3  2 


INERTIA 

PARAMETERS 

BODY 

AXES 

STABILITY  TO  BODY  AXES 

PRINCIPAL  TO  BODY  AXES 

Coefficient  of  Element 

Transformation  Equation 

Coefficient  of  Element 

Transformation  Equation 

sin’  ao 

COS*  Co 

sin  ao  cos  a0 

1 

sin’e 

cos’* 

sin  e  cos  e 

I 

lx 

K 

Ixg 

2IxZg 

0 

Ixp 

0 

0 

Ir 

0 

0 

0 

Irg 

n 

0 

0 

It;> 

Is 

K 

K 

— 2Ixz8 

0 

Ixp 

\ 

0 

0 

Ixz 

— IxZg 

Ixzg 

K  ^X8 

0 

0 

0 

Izp— lip 

0 

INERTIA 

PARAMETERS 

STABILITY 

AXES 

BODY  TO  STABILITY  AXES 

PRINCIPAL  TO  STABILITY  AXES  j 

Coefficient  of  Element 

Transformation  Equation 

Coefficient  of  Element  § 

Transformation  Equation 

sin*'a0 

COS*  «j 

siii  a0  cos  a0 

1 

sin’ij 

C08*1J 

sin  17  cos  17 

i 

Ixg 

Iz 

lx 

— 2  Ixz 

0 

Ixp 

0 

0 

Ir. 

0 

0 

0 

It 

0 

0 

0 

Iy 

K 

Iz 

21  xz 

0 

LP 

K 

0 

m 

Ixz, 

— Ixz 

Ixz 

lx — Iz 

0 

0 

0 

I^p  I^p 

0 

INERTIA 

PARAMETERS 

PRINCIPAL 

AXES 

BODY  TO  PRINCIPAL  AXES 

STABILITY  TO  PRINCIPAL  AXES 

Coefficient  of  Element 

Transformation  Equation 

Coefficient  of  Element 

Transformation  Equation 

sin*  t 

cos’  e 

sin  e  cos  « 

1 

sin’  p  cos*  ij 

sin  17  cos  17 

1 

Ixp 

Iz 

lx 

-  2Ixz 

0 

Izg  Ixg 

2Ixz, 

0 

Itp 

0 

0 

0 

It 

0  0 

0 

It, 

\ 

lx 

Iz 

2Ixz 

0 

lx,  Iz, 

— 2IxZg 

a 

0* 

— Ixz 

Ixz 

Ix-Iz 

0 

— Ixz,  Ixz, 

a 

0 

Write  equations  for  inertia  parameters  hy  summing  across  the  row. 

Example:  l;:g  =  Ix  sin*  a o-f-Iz  cos’  »«+  21  xz  sir  a»  cos  a<>+  (0)  1 
I;.-.,—  Ix|f  sin*  i)  +  Izp  cos’  (0)  sin  a»  cos  «,+ (0)  1 
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The  basic  stability  derivatives  may  be  arranged  in  the  six  matrices  listed  below.  With  the  assumption  of  vehicle  sym¬ 
metry  about  the  XZ-plane,  these  matrices  have  been  simplified  to  the  form  shown.  Derivatives  of  symmetric  forces 
with  respect  to  asymmetric  variable  are  neglected  and  certain  negligible  derivatives  are  taken  to  be  zero.  These 
considerations  are  trussed  further  .in  Section  4. 


Matrix  Type 

Matrix  Symbol 

- 

xB 

0 

X,  ■ 

Force-Velocity 

[Fv]  = 

0 

yt 

0 

0 

Zw  „ 

- 

0 

X,, 

0  ' 

Force  Rotary 

[F«J  = 

Y„ 

0 

Yr 

0 

t, 

0 

- 

0 

4 

0  * 

Force-Control 

[Fa]  = 

4 

0 

Y«r 

0 

4 

0 

— 

0 

4 

0  ■ 

Moment- V  elocity 

[Gv]  = 

M„ 

0 

Mw 

0 

Nr 

0 

- 

4 

0 

4  ‘ 

Moment-Rotary 

[Gj  = 

0 

M„ 

0 

NP 

0 

K  . 

- 

4, 

0 

4  " 

Moment-Control 

[Ga]  = 

0 

M«e 

0 

4 

0 

4. 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 


As  was  the  case  with  elements  of  the  inertia  matrix,  elements  of  the  velocity  and  rotary  matrices  may  be  trans¬ 
formed  according  to  the  relations  on  page  31.  The  control-derivative  matrices,  however,  require  specialized  treat¬ 
ment  because  the  control  displacements  are  independent  of  the  stability  and  body  axis  systems. 

A  procedure  similar  to  that  used  for  transformation  of  the  inertia  matrix  (page  29)  may  be  used  to  find  the  trans¬ 
formation  relations  for  the  control-derivative  matrices.  Consider  a  vector  relation  of  the  form 

f  —  fB  j  5 


in  which  the  components  of  the  displacement  8  are  independent  of  the  coordinate  system  used  to  define  the 
components  of  f.  Components  of  the  vector  f  may  be  transformed  from  one  coordinate  system  (subscript  o) 
to  another  (subscript  3),  which  has  been  rotated  about  the  origin.  The  vector  relation  above  may  be  expressed  as 


fo  =  [Bo]8  or  f3  =  [B,j  S 


One  of  the  vector  transformation  matrices  [r]  from  table  2  may  be  used  for  a  general  rotation.  Thus 

fa  — [r]f0  — [r]  [B0]  8  =  [B*]  8 

The  transformation  relation  for  the  matrix  [B]  is  then 

[B»]  =  [r]  [Bc] 

Exprnsion  of  this  equation  for  a  planar  rotation  about  the  Y-avii  results  in  the  general  form  used  in  transforming 
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the  stability  derivatives  that  involve  control  variables.  If  the  initial  and  transformed  matrices  are  denoted  by  [B] 
and  [B],  respectively,  and  the  transformation  matrix  by  [0],  the  transformation  relation  becomes 

[B]=[tf]  [B] 


or 


r  ~ 
b** 

b*y 

b*. 

"  cos 

9 

0 

— sin 

9  ~ 

‘  b,* 

b*y 

b**  ' 

byx 

by* 

= 

0 

1 

0 

by* 

byy 

by* 

**** 

_  b„ 

b*y 

r 

i _ 

_  sin 

9 

0 

cos 

9  _ 

_  bzx 

b*y 

[B]  = 


Expanding  the  right  side  of  this  equation  resuits  in  the  transformed  matrix  [B]. 

"  (bxx  cos  9  —  b„  sin  5)  (b*y  cos  9  ~  b.y  sin  9)  (b„  cos  9  —  b**  sin  9) 

[B]  =  b^  uJT  by* 

_  (bx*  sin  9  -f-  b„  cos  9)  {b*y  sin  8  +  b*r  cos  0)  (b**  sin  0  -f-  b„  cos  0) 


Finally,  identifying  elements  of  [B]  with  the  above  expansion  of  [0]  [B]  gives  the  equations  transforming 
the  elements  of  the  matrix  [B]  to  the  elements  of  [BJ.  These  equations  are 


b„  =  b„  cos  0  —  b**  sin  0 
b*y  =  b*y  cos  9  —  b*y  sin  0 
b„  =  b*,  cos  0  —  b„  sin  0 

by*  —  by*  ;  byy  =  byy  ;  fay,  =  by, 

b,*  =  b„  sin  0  +  b**  cos  0 
b,y  =  b*y  sin  9  -|-  b*y  cos  9 
bB  =  b„  sin  9  +  b«  cos  9 


(30) 


The  foregoing  equations  are  used  to  transform  the  force-control  and  moment-control  matrices.  Special  trans¬ 
formations  may  be  devised  as  required  by  using  the  method  of  the  preceding  development  or  that  used  for  the 
inertia  matrix  on  page  30. 


Stability  derivative  transformations  between  stability  and  body  axes  are  tabulated  in  tables  8  and  9.  The  angle  a„ 
is  defined  in  figure  17.  A  prime  ( ' )  is  used  in  these  tables  to  designate  the  derivatives  along  body  axes.  In  sub¬ 
sequent  Sections  the  prime  notation  is  deleted  and  the  reference  axes  are  as  noted  (see  table  17). 


The  following  example  illustrates  the  use  of  the  stability  derivative  transformation  relations. 

Given:  Stability  derivatives  with  respect  to  stability  axes. 

Find:  The  derivative  of  rolling  moment  with  respect  to  rolling  velocity  L'p  referred  ta  body 

Solution:  Write  equation  from  table  8  by  summing  terms  along  the  row  of  L',,. 

L'p  =  Nr  sin2  -f-  L,,  cos2  —  (L,  +  Np)  sin  cos  a« 
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axes. 


TABLE  8 


STABILITY-DERIVATIVE  TRANSFORMATIONS 
STABILITY  AXES  TO  BODY  AXES 


MATRIX 

MATRIX 

BODY- 

AXES 

COEFFICIENT  OF  ELEMENT-TRANSFORMATION  EQUATION 

TYPE 

SYMBOL 

tin  Up 

1 

FORM 

sinao 

cos  a® 

sin  ae  cos  a© 

COSO* 

X'. 

z. 

x„ 

-Xw-Z. 

0 

0 

0 

X'w 

-z. 

X. 

X.-Zw 

0 

0 

0 

FORCE  - 
VELOCITY 

<FV) 

0 

0  - 

0 

0 

0 

Y, 

-Xw 

z. 

x.-z. 

0 

0 

■mm 

x. 

z. 

XwfZ. 

0 

0 

B 

X', 

0 

0 

0 

-z, 

x. 

FORCE  - 

(F«) 

Y'p 

Bisk  /SS’S 

0 

0 

-Yr 

Yp 

ROTARY 

0 

0 

0 

Yp 

Y, 

0 

0 

0 

X, 

z. 

mi 

x4 

0 

0 

0 

-Z*. 

4. 

0 

FORCE- 

<F.) 

Y\ 

0 

0 

0 

0 

4. 

CONTROL 

y' 

x«p 

0 

0 

0 

0 

0 

Y*r 

Z'*e 

0 

0 

0 

X*. 

4. 

0 

L't 

0 

0 

0 

-Nr 

L. 

0 

MOMENT- 

<GVI 

M'. 

0 

~M. 

M. 

0 

VELOCITY 

M'w 

0 

0 

0 

M. 

M. 

0 

N’» 

0 

0 

0 

L. 

N. 

0 

■ag 

N, 

Lp 

— Lr— Np 

0 

0 

0 

MOMENT - 
ROTARY 

-Np 

Lr 

Lp-Nr 

0 

0 

0 

M\ 

0 

0 

0 

0 

0 

M, 

N'p 

— Lr 

Np 

Lp-Nr 

0 

0 

0 

N'r 

L 

Nr 

Lr+Np 

0 

0 

0 

L'*. 

0 

0 

0 

-N.. 

4. 

0 

L'«r 

0 

0 

—N‘r 

4, 

0 

MOMENT - 
CONTROL 

(G.) 

M'*, 

0 

0 

0 

0 

0 

M., 

N\ 

0 

0 

4. 

N* 

0 

0 

0 

0 

4, 

0 

j  a„  is  positive  «s  shown  on  page  3  2 

i 

Wind-tunnel  stability  axes  are  used  as  reference  axes  for  most  wind-tunnel  data.  It  is  therefore  necessary  to  trans¬ 
form  these  data  to  vehicle  stability  axes  for  body  axes )  before  using  them  in  analysis  of  the  motion  of  an  aircraft. 
These  axes  are  defined  in  Section  2.  The  transformation  is  a  simple  pitch  rotation  about  the  lateral  (Y)  axis.  The 
rotation  angle  {  is  the  angle  between  the  wind-tunnel-axis  (or  body-axis)  angle  of  attack  a  and  the  stability 
axis  reference  angle  (reference  12). 


Wind-tunnel  axes  may  also  be  considered  as  general  wind  axes.  In  this  case  the  transformation  from  wind-tunnel  . 
axes  to  wind-tunnel  stability  axes  is  a  simple  yaw  rotation  about  the  Z-axis  through  an  angle 


TABLE  9 

STABILITY-DERIVATIVE  TRANSFORMATIONS 
BODY  AXES  TO  STABILITY  AXES 


MATRIX 

MATRIX 

ma 

COEFFICIENT  OF  ELEMENT-TRANSFORMATION  EQUATION  | 

TYPE 

SYMBOL 

sin1  as 

cos’  a» 

sin  ffecoec* 

sin  ao 

cos  a* 

1 

X. 

Z'w 

X'u 

X'w  +  Z'. 

0 

0 

m 

X. 

-Z'. 

X’. 

Z'w-X'. 

0 

0 

i 

FORCE  - 
VELOCITY 

<FV) 

Yr 

0 

0 

0 

0 

0 

Y', 

z. 

-X'w 

Z'. 

Z'w-X'. 

0 

0 

mm 

zw 

X'. 

Z'w 

-Xw-Z, 

0 

0 

n 

X, 

0 

0 

0 

Z', 

X', 

FORCE- 

<FJ 

Y, 

0 

0 

0 

Y'r 

Y'p 

ROTARY 

Y, 

0 

0 

0 

— Y', 

z. 

0 

0 

0 

-X's 

Xv. 

0 

0 

0 

Z'«. 

X.e 

0 

FORCE#- 

(Fp) 

Y.m 

0 

0 

0 

0 

0 

Y'.. 

CONTROL 

Y*, 

0 

0 

0 

0 

0 

Y\ 

0 

0 

0 

-X*. 

Z*. 

0 

L, 

0 

a 

0 

N'r 

L'r 

MOMENT- 

(Gy) 

Mu 

0 

0 

0 

M'w 

M'. 

VELOCITY 

M„ 

0 

0 

— M', 

M'w 

mm 

Nt 

0 

0 

0 

-L'r 

N'r 

m 

u 

N'r 

L'p 

L'r  +  N'p 

0 

0 

0 

L, 

— N', 

L'r 

N'r -L'p 

0 

0 

0 

MOMENT - 
ROTARY 

(GJf 

M, 

0 

0 

0 

0 

0 

M', 

NP 

-L'r 

N'p 

N'r -L'p 

0 

Nr 

L'p 

N'r 

-L'r -N'p 

0 

0 

Hfl 

c.. 

0 

0 

0 

N's. 

L'p. 

0 

L.r 

0 

0 

0 

N's, 

L'.. 

0 

MOMENT- 

(G.) 

Ms. 

0 

0 

0 

0 

0 

M'p, 

CONTROL 

0 

-L's. 

N'p. 

Ns, 

0 

0 

0 

Nsr 

0 

0 

0 

— L'sj. 

N'.r 

0 

a0  is  positive  as  shown  on  page  a  3 


TABLE  iO 

VECTOR  TRANSFORMATION  MATRICES 

WIND- TUNNEI,  AXES  TO  WIND-TUNNEL  STABILITY  AXES  TO  VEHICLE  STABILITY  AXES 


WIND-TUNNEL  AXES 

TO 

WIND-TUNNEL  STAB. 
AXES 

VECTOR  COMPONENTS -WIND-TUNNEL  AXES 

Xw 

Yw 

Zw 

VECTOR 

COMPONENTS 

WIND-TUNNEL 

STABILITY 

AXES 

XwT 

cos  f 

sin  f 

0 

Ywt 

—sin  f 

CM 

0 

ZwT 

0 

0 

1 

WIND-TUNNEL  STAB. 
AXES 

TO 

VEHICLE  STABILITY  AXES 

VECTOR  COMPONENTS -WIND-TUNNEL  STAB.  AXES 

XwT 

Ywt 

ZwT 

VECTOR 

COMPONENTS 

VEHICLE 

STABILITY 

AXES 

X. 

COS  { 

0 

—sin  { 

Y. 

0 

1 

0 

a 

•taf 

• 

CMf 

The  relations  between  wind-tunnel  axes  and  vehicle  axes  ( stability  and  body)  are  illustrated  in  following  figure 
Table  10  gives  the  vector  transformation  matrices  in  tabular  form. 


(SAME  AS  GENERAL  WIND  AXE8) 


FIGURE  18  ROTATION  FROM  WIND-TUNNEL  TO  VEHICLE  AXES— YAW-PITCH  SEQUENCE 

2.  MULTIPLE  ROTATION 

There  are  several  transformations  used  in  the  analysis  of  vehicle  .motion  that  involve  multiple  rotations  of  an  axis 
system.  Transformations  are  given  in  this  Section  that  facilitate  the  changing  of  vector  components  between  earth 
and  body  axes  and  between  wind  end  body  axes.  The  relations  between  these  axis  systems  define  the  orientation 
of  a  vehicle  with  respect  to  the  earth  and  the  relative  velocity  or  flight  path. 

The  axes  and  notation  used  are  defined  in  Section  2,  and  the  specific  vector  transformation  matrices  are  obtained 
from  the  genevul  cases  in  table  2. 

Only  transformations  for  vector  components  are  included  in  this  Section.  As  is  noted  previously,  earth  and  wind 
axes  are  not  convenient  reference  axes  for  rigid-body  motion  analysis.  This  results  from  the  fact  that  the  inertia 
parameters  become  unnecessarily  complex  functions  of  time  and  the  orientation  angles  (aee  page  16). 

Earth  axes  are  used  primarily  as  a  reference  for  the  gwvity  force  and  the  description  of  vehicle  motion  over  a 
long  period  of  time.  Ihe  orientation  angles  between  moving  earth  axes  and  body  axes  are  defined  on  page  11  and 
are  shown  in  figure  19.  The  rotation  sequence  corresponds  to  Case  1  of  table  2.  Hence  the  vector  transformation 
matrix  is  obtained  by  substitution  of  ^  =  ♦,  —  «,  and  =  *  in  the  matrix  for  Case  1.  b  the  angular-velocity 

relations  yhe  body-axis  notation,  P,  Q,  and  R,  is  aaed  instead  of  9s»  K  and  respectively. 
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REFERENCE 

HEADING 


■I 


FIGURE  19  ROTATION  FROM  EARTH  TO  BOOT  AXES  — 
YAW  •  PITCH  -  ROLL  SEQUENCE 


Z 


i 


FIXED  EARTH  AXES 


TABLE  11 

VECTOR  TRANSFORMATION  MATRIX 
EARTH  TO  BODY  AXES 


COMPONENTS  ALONG  EARTH  AXES 

X. 

Y. 

Z. 

g 

X 

cos  ©  oo*  ♦ 

cot  8  tin  ♦ 

—tin  8 

g 

1 

Y 

tilt  •  tin  8  cot  ♦ 
—cot#  tin  ♦ 

tin  #  tin  6  tin  # 

4-cot  ♦  cot  * 

tin  #eotO 

Z 

cot  #  tin  6  cot  # 
-ftm*tia* 

cot  ♦  tin  0  tin  * 

— tin#  cot# 

cot  #  oot  8 
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TABLE  12 


ANGULAR-VELOCITY  RELATIONS 
EARTH  TO  BODY  AXES 


P  =  ♦  —  ♦  sin  0 

'i'  =  Qsin*Kc0-(-Rcos*aec9 

Q  =  0cos  +  +  *8in  +  co89 

0  =  Qcos  +  —  Rsin* 

R  =  —  0»in*  +  'i,cos*cos0 

*  P  -f  Q  sin  ♦  Un  0  R  co»  +  tan  0 

General  Wind  Axes  to  Body  Axes 

General  wind  axes  are  oriented  with  respect  to  the  relative  wind.  Hie  orientation  angles  relating  general  wind 
axes  to  vehicle  body  axes  are  therefore  convenient  variables  to  use  in  expressing  the  aerodynamic  characteristics 
of  a  vehicle.  General  wind  axes  are  defined  in  figure  8.  The  yaw-pitch  definition,  the  preferred  definition,  of  angle 
of  attack  a  and  sideslip  angle  £  is  illustrated  in  figure  20.  The  pitch-yaw  sequence  is  illustrated  in  figure  21. 
Vector  transformation  matrices  corresponding  to  these  definitions  are  given  in  tables  13  and  15,  respectively. 
The  angular-velocity  relations  are  given  in  tables  14  and  16,  respectively. 


noun  M  ROTATION  FBOM  VOID  TO  BOOT  AXU  —  TAV  -  NITON  SBOOBNCB 
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TABLE  13 


VECTOR  TRANSFORMATION  MATRIX 
GENERAL  WIND  TO  BODY  AXES 


COMPONENTS  ALONG  WIND  AXES  j 

X, 

Y* 

Zw 

gg 

D 

cos  a  cos  £ 

—cos  a  sin  ft 

—sin  a 

a 

sin  ft 

COS0 

0 

o 

z 

sin  a  cos  0 

—sin  a  sin  0 

cos  a 

TABLE  14 

ANoULAR- VELOCITY  RELATIONS 
GENERAL  WIND  TO  BODY  AXES 


P  =  ft  sin  a 

ft  =  —  R  sec  a  =  P  esc  a 

Q  =  i 

i  =  Q 

R  =  —  ft  cos  a 

0  =  P+R  tan  a 

Note:  Tse  ebove  matrix  and  equations  result  from  substitution  of  —  ft,  a, and  0  for  f ,  #i,  and 
respective]),  in  Cat-a  1  of  Table  2 


TABLE  15 

VECTOR  TRANSFORMATION  MATRIX 
GENERAL  WIND  TO  BODY  AXES 


j  COMPONENTS  ALONG  WIND  AXES  j 

Xv 

Y* 

z. 

si 

X 

cm  7  am  (l 

At 

—sin  0 

—Sm  a  cm  f) 

>•£ 

as 

Y 

cm*  tirnfi 

COS  P 

,  ^ 

-i*  m  $m  p 

*8 

Z 

0 

cm  r 
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TABfife  16 


ANGULAR-VELOCITY  RELATIONS 
GENERAL  WIND  TO  BODY  AXES 


P  =  —a  sin  p" 

"a  —  Q  sec  — P  esc  y 

Q  =  gem  JJ” 

T=-R 

R  =  -? 

0  =  P  +  0  Un  'p 

om  substitution  of  a,  — /3,  and  0, 
in  Case  4  of  table  2. 


TOTAL  VELOCITY 
^  V 


SECTION  4.  REAL  FORCES  AND  MOMENTS* 


The  preceding  Sections  contain  equations  of  motion  and  transformations  useful  in  the  analysis  of  particle  and  rigid-body 
motion.  Real  force  and  moment  components  are  indicated  in  these  equations  by  a  general  notation  for  forces  and 
moments.  In  the  following  Section  the  force  and  moment  components  are  presented  more  specifically.  The  general  com¬ 
ponent  expressions  are  expanded  to  show  contributions  of  gravity,  aerodynamic  force,  and  direct  thrust  force.  Aero¬ 
dynamic  force  and  moment  components  are  then  further  expanded  for  the  case  of  small  disturbances. 

Stability  derivatives  for  airpiane-type  vehicles  are  summarized . 


GENERAL  FORCE  AND  MOMENT  DESCRIPTION 


The  particular  flight  path  or  motion  of  a  rigid  body  is  the  result  of  the  external  forces  and  moments  that  are  applied. 
Thus  the  applied  forces  and  moments  may  be  considered  as  the  “driving  functions”  to  be  used  with  the  equations  of 
motion  of  the  vehicle.  Solution  of  these  equations  then  provides  the  motion  or  response  of  the  vehicle  to  the  applied 
forces  and  moments.  Inversely,  the  problem  may  be  formulated  to  find  the  force  and  moment  input  required  to 
accomplish  a  specified  motion. 

The  real  forces  and  moments  involved  in  the  motion  of  a  body  through  the  atmosphere,  in  the  gravitational  field  of  the 
earth,  may  be  separated  into  contributions  of  gravity,  aerodynamic  force,  and  direct  thrust.  In  the  case  of  particle 
motion,  moments  about  the  center  of  mass  are  zero,  and  only  the  force  vectors  need  to  be  considered. 

Components  of  the  external  force  and  moment  vectors  are  usually  resolved  along  vehicle  body  axes .  Relations  trans¬ 
forming  these  vector  components  to  the  body  axis  system,  or  any  other  desired  reference  axis  system,  may  be  obtained 
directly  or  derived  from  the  preceding  Section  . 

Separation  of  the  force  and  moment  vectors  into  gravity,  aerodynamic,  and  direct  thrust  contributions  is  outlined  below. 
For  illustrative  purposes  the  general  force  vector  F  and  general  moment  vector  G  are  resolved  into  components  along 
vehicle  body  axes.  Body  axes  are  usually  the  most  convenient  reference  axes. 

The  general  vectors  are  resolved  into  components  along  the  reference  axes. 

Thus 

F  =  Fxl  -4-  Fy  j  -f  F,k 

(31) 

G  =  G*i  +  Gyj  +  G,k 


Separation  of  the  several  components  into  gravity,  aerodynamic,  and  direct  thrust  contributions  results  in  the  following 
equations: 


F,  —  Xj-f-X-f  Xt 

F,  =  Y,  +  Y  +  Yt  • 

F.  =  Z,  +  Z  +  Zr  . 

G,  =  l*  4-  L  +  L* 

Gy  —  Mg  -f-  M  -f"  M«r  ■ 
G ,  =  N,  4-  N  4-  Nr 


(32) 


(33) 


where 


Xg,  Y„  Zg  are  components  of  gravity  force  along  reference  axes 

M«,  Ng  an  moment  components  about  referer^e  axes  due  to  gravity  force.  (These  are  usually  aero.) 

*  The  tera  “renT  is  need  to  deai ‘r" e  nrniniriisl  forces  and  woassnt*.  Thus  the  apparent  forces  mch  as  centrifagsi  force  or  Csrtslis  force 
are  rxdnoed.  Gravitational  sr  efectromagnstk  forces,  prapnisive  system  throat,  and  aerodynamic  force  ate  examples  W  “real*  farces. 
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X,  Y,Z 
L,M,N 
Xt,  Yt,  Zt 


are  components  of  aerodynamic  force  along  reference  axes 
are  moment  components  about  reference  axes  due  to  aerodynamic  force 
are  components  of  direct  thrust  force  aloig  reference  axes 
Lr,  Mf,  Nt  ere  moment  components  about  reference  axes  due  to  direct  thrust  force 

Note:  YT,  Lr,  vad  WT  are  usually  zero  because  of  vehicle  symmetry. 

GRAVITY- FORCE  COMPONENTS 

The  gravitational  force  upon  a  vehicle  is  most  naturally  given  in  terms  of  earth  axes.  With  respect  to  earth  axes 
the  gravity  vector  mg  is  directed  along  the  Ze-axis  (page  10).  Components  along  vehicle  body  axes  are  readily 
obtained  by  using  the  transformation  given  in  table  11.  The  gravity-force  components  are 

X,  =  —mg  sin  8 
Yt  =  mg  cos  8  sin  * 

Zg  —  mg  cos  0  cos 

along  the  vehicle  body  axes  X,  Y,  and  Z,  respectively. 

There  are  no  moments  resulting  from  the  gravity  force  uhen  the  origin  coincides  with  the  vehicle  center  of  gravity. 
However,  if  the  origin  is  displaced  from  tire  rc-zb-i  o t  gravity,  the  same  transformation  (table  11)  may  be  applied 
to  the  components  of  the  gravitatiosa!  moment  about  the  origin  to  obtain  Lg,  M#,  and  N  . 

The  components  of  gravitational  force  upt-a  the  r chicle  are  functions  of  the  vehicle  pitch  and  roll  attitudes  only. 
Heading  angle  does  not  affect  the  resolution  of  the  gravity  force  to  body  axes. 

AERODYNAMIC  FORCES  AND  MOMENTS 

The  contributions  of  aerodynamic  fc-ree  to  the  general  force  and  moment  vector  components  are  outlined  in  this 
Section.  These  components  are  referred  to  vehicle  body  axes  (Section  2,  pages  11  and  12)  and  aircraft  terminology 
and  notation  are  used.  Lift  and  (bag  forces  are  thus  introduced  and  transformed  to  die  body  axis  system.  The 
general  form  for  expansion  of  aerodynamic  terms  for  small  disturbances  is  included  in  this  Section  also 

It  is  convenient  to  use  dimensionless  coefficients  to  describe  the  behavior  of  aerodynamic  forces  and  moments. 
These  coefficients,  defined  according  to  established  usage,  are  discussed  and  analyzed  in  aerodynamic  texts  and 
in  reference  13.  The  aerodynamic  parameters  and  their  derivatives  should  be  evaluated  from  experimental  data, 
i.e.,  wind-tunnel  or  flight-test  data,  or  estimated  from  the  appropriate  data  given  in  reference  13  or  a  similar 
source.  Mach  number  and  Reynolds  number  effects  upon  aerodynamic  parameters  are  assumed  to  be  included. 

Induced  effects  of  the  propulsive  system  are  assumed  to  be  included  in  the  aerodynamic  coefficients,  since  these 
effects  are  normally  included  in  complete-model  wind-tunnel  data.  Direct  thrust  forces  and  moments  are  discussed 
later,  on  page  50. 

Aerodynamic  coeftcieeta  depend  upon  the  orientation  of  the  relative  wind  or  velocity  vector  with  rerpect  to  the  body. 
The  nogbw  of  attack  a  and  lideaiip  ft  which  define  this  orientation,  are  thus  convenient  independent  variables  fer 
expressing  the  variations  of  tha  aerodynamic  characteristics  of  a  body.  Them  angles  a**o  determine  the  velocity  com 
patents  U,  V,  and  W  along  the  rehick  reference  axes  X.  Y,  and  Z,  respectively. 

The  eagles  of  attack  sad  sideslip  are  shown  in  figure  22  with  the  velocity  components  alor,*  the  reference  uses. 
Both  the  yaw-pitch  rotation  (a,  ft)  mi  the  pitch-yaw  rotation  la,  ft)  are  given;  however,  the  former  is  the  preferred 
rotation.  The  relations  far  the  velocity  components  II,  V,  and  W  follow  directly  from  the  transformations  from 
general  wind  axes  to  vehicle  body  axes  in  tables  13.  14.  15.  and  16. 


AS 


TW  KMisiMt  of  the  total  M^rMnc  force  i«  the  vehicle  plane  of  ey— atn  is  shorn  is  figure  23.  Lift,  l.,  i«i 
drug,  D,  sre  the  familiar  aerodynamic  forces  a  ones  1  sad  parallel,  respectively,  to  V  cos  /9,  the  component  of  the 
total  velocity  is  the  vehicle  plane  of  symmetry.  Alternatively,  lift  and  drag  my  be  defined  as  the  aerodynamic 
force  component*  in  the  plane  of  symmetry  along  *  instantaneous*  liability  axes,  symmetric  wind  axes ,  or  wind- 
tunnel  stability  axes.  It  should  be  noted  that  lift  and  drag  are  defined  to  be  positive  as  illustrated.  Thus  these 
quantities  have  a  negative  sfenae  with  respect  to  the  usual  vehicle  axis  systems. 

The  relations  for  the  aerodynamic  force  components  along  body  axes  are  included  with  figure  23.  These  equations 
may  be  obtained  directly  from  this  figure  or  from  the  vector  transformation  from  stability  axes  to  body  axes 
(table  6). 


Aerodynamic  forces  and  moments  are  usually  given  in  terms  of  basic  aerodynamic  coefficients.  These  coefficients  are 


defined  by  the  following  relations: 

/ 

L*  =  CL^S 

lift  force 

L*  =  CjqSb 

rolling  moment  about  X-axis 

D  —  CDqS 

drag  force 

M  =  CmqSc 

pitching  moment  about  Y-axis 

C  =CcqS 

cross-wind  force 

N  =  CBqSb 

yawing  moment  about  Z-axis 

X  =iOxqS 

aerodynamic  force  along  X-axis 

Y  =  CynS 

aerodynamic  fore®  along  Y  axis 

Z  —  Cg]S 

l  .  v  .  - 

aerodynamic  force  along  Z-axis 

* 

where  the  quantities  in  the  above  expressions  are  defined  as  follows: 


Cl*  Cd,  Qj.Cx,  Cy,  Cz 
Cl*  Cm,  CB 


aerodynamic  force  coefficients 
moment  coefficients 

dynamic  pressure 


P 


atmospheric  density 


y 

s 

b 


total  velocity 

reference  area  (usually  wing  area) 
reference  length  (usually  wing  span) 
reference  length  (usually  wing  M.AjC.) 


The  force  coefficients  C x,  CY,  and  Cz  are  expressed  in  terms  of  lift,  drug,  and  cross-wind  coefficients  by  the 
same  equations  that  relate  the  forces  in  figure  23.  Thus 

*  Duplication  of  tile  symbol  L  for  lift  and  rolling  moment  has  persisted  throughout  the  aircraft  industry  dace  early  times.  However,  since 
coefficients  are  usually  used,  this  ambiguity  is  avoided  by  taking  L  as  ho  subscript  in  the  lift  .coefficient  sad  /  as  the  subscript  in  the 
rolling-moment  coefficient. 
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Co  esc  a 


Cj  ~  Ci  sic  a  — 


C?  —  Cc 

Cg  —  — Cl  cos  *  —  Co  sis  a 


Co  ~  — C*  esc  a  —  Co  sic  •  t 

Cc  =  Ct 

Cl  —  Cj  sin  a  —  Co  cos  a 


(35) 


The  moment  components,  as  treated  in  this  Section,  are  defined  with  respect  to  the  body  axis  system  and  as  such  may  be 
used  directly.  However,  it  should  be  noted  that  moments  and  moment  coefficients  may  be  defined  with  respect  to 
stability  or  wind-tunnel  axes  and  that  in  these  cases  the  appropriate  transformation  from  Section  3  must  be  used  to 
obtain  the  desired  aerodynamic  moment  components. 

A  word  of  caution  is  in  order  concerning  the  transformation  of  moment  coefficients.  The  reference  length?  and  are?? 
used  in  defining  the  moment  coefficients  may  be  different  and,  if  so,  this  difference  must  be  accounted  for  in  the 
transformation  of  the  moment  coefficients  from  one  axis  system  to  another.  Thus,  although  the  transformation  is 
appropriate  for  the  vector  components,  it  does  not  directly  transform  these  components  when  expressed  in  coefficient 
form.  Another  important  item  relative  to  moment  coefficients  is  the  location  of  the  moment  reference  center.  In  any 
particular  case  this  must  be  checked  to  assure  that  (he  moment-of-inertia  and  product-of-iriertia  terms  on  the  right  side 
of  the  equations  of  motion,  e.g.,  equation  9  ,  are  consistent  with  the  center  of  reference  for  the  external  moments. 

The  aerodynamic  forces  and  moments  are  involved  functions  of  many  variables.  Test  data  are  the  best  source  of 
aerodynamic  force  and  moment  characteristics;  however,  in  many  instances  a  particular  configuration  may  be  in  the 
preliminary  design  stage  and  test  data  may  not  yet  be  available.  When  it  is  required  to  estimate  aerodynamic  charac¬ 
teristics  of  a  configuration,  data  and  techniques  such  as  those  contained  in  reference  13  should  be  used. 

A  summary  of  the  major  variables  that  affect  the  aerodynamic  characteristics  of  a  rigid  body  or  a  vehicle  is  given  below. 

Velocity,  temperature,  and  altitude:  These  variables  may  be  considered  directly  or  indirectly  as  Mach  numbers, 
Reynolds  numbers,  and  dynamic  pressures.  Velocity  may  be  resolved  into  components  U,  V,  and  W  along  the 
vehicle  reference  axes. 

Angles  of  attack,  and  sideslip:  Angle  of  attack  a  and  angle  of  sideslip  /£>  may  be  used  with  the  magnitude  of  the 
total  velocity  V  to  express  the  orthogonal  velocity  components  U,  V,  and  W.  It  is  more  convenient  to  express 
variation  of  force  and  moment  characteristics  with  these  angles  as  independent  variables  rather  than  the 
velocity  components. 

Angular  velocity:  The  angular  velocity  is  usually  resolved  into  components  PT  Q,  and  R  about  the  vehicle  reference 
axes. 

Control-surface  deflection:  Control  surfaces  are  used  primarily  to  change  or  balance  aerodynamic  forces  and  moments. 


Since  the  above  variables  are  identified  with  a  steady  motion,  the  variation  of  aerodynamic  forces  and  moments 
with  time  is  asau  ned  to  be  negligible.  Asia  noted  in  references  4,  11,  and  14,  this  assumption  is  reasonable  for 
most  problems  in  analysis  of  vehicle  motion  in  the  atmosphere.  However,  aerodynamic  forces  and  moments  are  the 
result  of  the  pressure  of  the  air  exerted  on  the  body  and  this  pressure  depends  upon  the  flow  field  about  the  body. 
Because  air  has  mass,  the  flow  field  cannot  adjust  instantaneously  to  sudden  changes  in  these  variables,  and 
transient  conditions  exist.  In  some  cases,  these  transient  effects  become  significant.  Analysis  of  certain  unsteady 
motions  may  therefore  require  consideration  of  the  time  derivatives  of  the  variables  listed  above. 

Two  typical  functional -depend®  ice  relations  for  the  aerodynamic  force  component  along  the  body  X-axis  are  expressed 
below.  Similar  expressions  Y  and  2  force  components  and  the  aerodynamic  moment  component?  L,  M  and  N 
could  be  written: 


X  =  Xi  (U,  V,  W,  U,  V,  W, ....  P,  Q,  R,  P, Q,  R, 8n  in  •■•*/*» M,  R,*. . .)  ] 
—  Xs  (V,  a,  j9,  V,  a,  /9, . . . ,  P,  Q,  R,  P,  Q,  R,  • .  •  >  8i»  8*,  In  I«>  I?  In  •  •  • » p*M,  R»  •  •  ■) 


Ms MACH  NUMBER 
R-RBTNOLDS  NUMB3B 


(36) 
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It  is  apparent  that  in  the  practical  ewe  of  solving  engineering  problems  much  simplification  of  the  above  functional 
relations  is  required.  Fortunately,  it  is  possible  to  make  certain  assumptions  that  simplify  the  mathematics  considerably, 
while  still  permitting  solutions  of  practical  significance.  Expansion  of  the  aerodynamic  components  for  small  disturbances 
is  illustrated  in  subsequent  paragraphs,  and  the  simplifications  of  the  equations  of  motion  for  practical  solutions  are 
discussed  in  Section  S  . 

There  are  many  cases  of  practical  interest  in  the  analysis  of  aircraft  motion  in  which  the  disturbance  from  a  steady-flight 
condition  ia  small.  In  these  cases  it  is  permissible  and  convenient  to  express  the  aerodynamic  force  and  moment  com¬ 
ponents  in  a  Taylor  so;  ies  expansion.  This  expansion  is  formed  in  terms  of  perturbations  from  a  reference  steady-flight 
condition.  The  use  of  the  expansion  is  limited  to  problems  where  the  perturbations  are  small  and  where  the  second-  and 
higher-order  derivatives  of  the  variable  quantities  and  the  products  of  the  perturbation  quantities  are  therefore  negligible. 
Thus  they  may  be  omitted  in  the  simplified  expressions  for  the  aerodynamic  force  and  moment  components.  Of  course, 
this  type  of  expansion  requires  the  first  derivatives  of  the  aerodynamic  force  and  moment  components  with  resp-ct  to 
the  aerodynamic  variables  that  affect  these  forces  and  moments.  These  derivatives,  commonly  referred  to  as  “stability 
derivatives,”  will  be  discussed  in  detail  later. 

The  general  procedure  for  expanding  the  aerodynamic  force  and  moment  relations  for  small  disturbances  from  steady 
flight  may  be  found  in  many  places  in  the  literature,  e.g.,  references  4  and  11  .  The  notation  for  small 
disturbances  from  a  reference  flight  condition  is  as  follows.  Lower  case  symbols  are  used  to  designate  the  perturbations 
of  velocity  and  orientation  variables.  Upper  case  symbols  with  subscript  xero  are  used  to  denote  the  reference  values 
of  these  quantities.  For  example,  U  =  U0  -j-  u,  P  =  P0  +  p,  and  ♦  =  ♦<>  +  The  incremental  changes  in  the  aero¬ 
dynamic  force  and  moment  components  are  indicated  by  the  prefix  A.  Thus  X  =  X„  +  AX,  Y  =  Y,  +  aY,  M  =  M, 
-j-  aM,  etc.  The  control  deflection  angles  8„  8„  and  8,  are  used  just  as  they  are  and  are  interpreted  as  the  control 
perturbation  angles  from  the  steady-state  trimmed-flight  condition.  To  illustrate  the  general  expansion  of  the  aerodynamic 
force  component  along  the  body  X-axis  the  first  functional  relation  for  X  on  page  47  may  be  expressed  as  follows: 

X  =  X«  +  ax  =  X.  +  X.U  -I-  XTV  -1-  xww  -I-  Xii  +  X*v  +  Xiw  + . , .  + 

Xpp  -J-  +  X,r  -f-  Xjp  4*  X^q  4*  Xjr ...  4*  Xsa  8* 

Xs#  8.  -|-  Xsr  8,  -f  X«m  8,  4- . . .  -f  (higher  order  terms)  * 
where  X,  =  iX/  du,  etc. 

Variables  may  be  added  to  represent  additional  aerodynamic  effects  such  as  occur  in  higher  order  unsteady  aerodynamics. 


(37) 


l 

1  "a'1 

*  Just  the  second-order  terms  ia  the  expension  would  require  >  page  or  more  to  write  out  The  ex  pension  woeld  have  the  following  form: 
1/2  [(XoaU*  +  Xr»v*  ■+••••+  XsfU*  +  .  -  .  +  X**p’  -f-  .  .  .  +  Xjjp'  +  .  .  .  + 

Xv.*.1  +...-+•  Xv^.'  4- . . .)  4-2  (Xmiiv  -f- . . .  +  X»cev  4"  •  •  •  4* 

X-*p*  -J-  . . .  -p  X4pr  4  •  •  •  4-  X*«,l.v4- . . .  4-  X4,l.v  4  . . .)  4- 
2<X„;uw  4- ...  4-  X.,pw  4- ...  4"  -f-  •  ■  •  4"  X»«fl»w 4-  •  •  •)  + 

2(Xuvu  4- ...  4-  -f  •  •  •)  +  •  ••  4-  2  (X^  W.  +  • . .  4-)  4-  •  •  •  ] 
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TABLE  17 


AERODYNAMIC  FORCE  AND  MOMENT  INCREMENTS 

FOR  SMALL  DISTURBANCES 
| Components  Along  Vehicle  Axes) 


DISTURBANCE 

VARIABLES 


Linear 

Velocity 

Components 


Linear 

Acceleration 

Components 


Angular 

Velocity 

Components 


Angular 

Acceleration 

Components 


Control 

Deflection 


Control 

Deflection 

Rates 


Force  Components 

AY 


Moment  Components 
AL  AM  AN 


4*** 

M" 

Nj.  *** 

4* 

Mi*’* 

Ni* 

L;*** 

Mr** 

N;  *** 

K* 

K 

V 

K' 

y<;** 

Za." 

x«;** 

7,5  •  • 

Le 

Xlr” 

Y^** 

Z*r  * 

ur 

u; 

u 


Mi; 
Mi;** 
Mir'» 


Ni; 

n,; 

Ni  — 


•  AiyamWo  loro*  and  noaul  ooapoaaat  derivative,  with  raapaot  to  ayanaatrio  di  at  urban  oe  varlablea. 
**  Syaattrlo  foroa  and  moaiant  oonpoaant  derivatives  with  raapaot  to  aay  mate  trio  dlaturbaaoa  varlablea. 
(Both  aata  of  taiaa  ara  idantloally  aaro  tor  dlaturbaaoaa  from  a  atata  of  ataady  aymnatrlo  notion.  I. a., 
V,  P,  B,  and  thalr  darlvatlvaa  ara  all  aaro,  of  an  airoraft  with  an  XZ-plana  of  ayanatry.) 

•••  Taraa  tha4  ara  usually  aapllgtbla. 


THRUST  FORCES  AND  MOMENTS 

The  propulsive  system  of  a  vehicle  generally  produces  both  a  direct  thrust  force  and  indirect  or  induced  effects  upon  the 
aerodynamic  forces.  These  contributions  of  the  propulsive  system  to  the  force  and  moment  components  are  presented 
and  discussed  in  reference  13. 

Direct  thrust  force  and  moment  components  should  be  used  in  accordance  with  the  force-moment  component  resolu¬ 
tion  of  equations  32  and  33  on  page  44,  i.e.,  XT,  YT,  ZT,  LT,  Mx> and  N  r.  These  components  may  be  developed 
directly  from  the  georoetriy  relation  between  the  direct-thrust  line  of  action  and  the  moment  reference  center  of  the 
vehicle.. 


Induced -thrust  and  propulsive-system  effects  are  conveniently  included  in  the  aerodynamic  components.  Wind- 
tunnel  and  flight-test  aerodynamic  data  usually  include  the  indirect  effects  of  the  jet  flow  or  running  propellers 
upon  the  aerodyuamic  characteristics  of  a  vehicle  configuration.  Methods  for  estimating  these  induced  propulsive- 
system  effects  are  included  in  reference  13. 
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STABILITY  DERIVATIVES 

V 

The  partial  derivatives  that  occur  in  the  expansion  of  aerodynamic  force  and  moment  components  are  commonly  referred 
to  as  “Stability  Derivatives.”  These  quantities  are  useful  in  general  ind  most  applicable  in  the  analysis  of  small- 
disturbance  motions  from  a  steady  reference  flight  condition.  In  this  Section  the  derivatives  and  notation  are  defined. 


Several  systems  of  notation  and  definition  for  stability  derivatives  and/or  parameters  have  been  developed  and  are 
found  in  the  literature.  The  reader  should  be  forewarned  and  reminded  to  be  thorough  and  alert  when  using  published 
works  —  including  the  present  report  —  to  check  the  notation  and  definitions  used  (see  page  129  of  reference  4  or  page 
IV-2  of  reference  11).  The  notation  in  this  report  is  selected  to  be  consistent,  insofar  as  possible,  with  that  used 
by  NASA  and  in  references  3  and  4. 


A  summary  of  the  stability  derivatives  and  notation  used  in  vehicle  stability  and  control  analysis  is  presented  in 
this  Section.  Since  symmetry  of  the  vehicle  and  initial  flight  condition  is  assumed  throughout,  the  number  of  de¬ 
rivatives  is  reduced,  as  indicated  in  table  17.  Unusual  configurations  and  special  problems  may  require  deriva¬ 
tives  that  are  not  included  in  this  Section. 


Data  presented  in  reference  13  include  gen'  ul  and  detailed  information  on  stability  derivatives  and  methods  for 
estimating  values  of  these  derivatives.  Also,  general  and  special  methods  and  analyses  for  evaluation  of  these 
derivatives  are  found  in  many  places  in  the  literature., References  4,  7,  and  15  are  typical  general  references  and 
reference  17  is  a  typical  special  investigation  - 


Three  types  of  stability  derivatives  are  used  in  airplane  stability  and  control  analysis.  The  following  paragraphs 
discuss  each  type.  Notation  for  these  derivatives  and  other  items  used  in  connection  with  stability  analysis  are 
given  in  table  18.  Table  19  contains  the  relationships  used  to  define  the  nondimensional  derivatives  along  body 
axes  and  expresses  them  in  terms  of  stability  axis  derivatives.  Nondimensional  derivatives  along  stability  axes 
are  presented  and  identified  in  table  20.  This  is  the  most  familiar  form  of  the  stability  derivatives. 


DIMENSIONAL  DERIVATIVES -BODY  AXES 


The  derivatives  used  in  the  development  of  small-disturbance  expansions  for  aerodynamic  force  and  moment  rela¬ 
tions  (see  page  49)  were  defined  as  dimensional  derivatives.  These  partial  derivatives  of  the  force  and  moment 
components  are  taken  with  respect  to  perturbations  of  the  significant  velocity,  acceleration,  and  control  variables*. 


*  Coatidn  the  following  for  the  cnee  of  emnli  di*t whence  from  n  steady  Sight  coedition: 

Let  U  =  U0  +  u  sad  P  =  P0  -f-  p 


since 


au,  .  ap, 
au  “d  ap 


ax 

ax 

au. 

t  ax 

au 

ax 

au  “ 

au.. 

au 

+  17 

au  “ 

au 

dN 

aP  ~ 

aN 

ap„ 

apa 

ap 

+  aN 
+  *p 

Je_  _ 

ap  ~ 

aN 

dp 

au 


tre  each  sero  and  and  ~p~~  each  equal  1 
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The  genera!  notation  uses  upper  rate  symbols  with  a  subscript  denoting  the  variable  of  differentiation,  and  body  axes 
are  specified.  In  this  case,  however,  body  axes  may  refer  to  any  axis  system  fixed  to  the  vehicle  and  thus  include  the 
special  cases  of  stability  and  principal  axes.  This,  at  least  in  part,  is  the  origin  of  much  of  the  confusion  in  stability 
derivative  notation.  Also,  it  is  the  reason  for  using  the  prime  to  denote  the  difference  between  body  and  stability  reference 
axes  in  the  transformation  of  stability  derivatives  (tables  8  and  9). 


Dimensional  derivatives  are  used  as  elements  of  the  matrices  shown  in  equations  24  through  29  and  in  various 
equations  throughout  subsequent  Sections.  Some  examples  of  these  derivatives  are 


v  dX  u  M  m  dN  dZ 

X„=-;  M.=  — ;  N„=  — ;  ^ 

i 

The  dimensional  derivatives  are  listed  in  table  17,  with  the  notation  given  in  table  13. 


NONDIMENSIONAL  DERIVATIVES  -  BODY  AXES 


The  use  of  nondimensional  equations  is  usually  convenient  when  aerodynamic  forces  and  moments  are  involved  in  a 
problem.  Hence  it  is  useful  to  define  nondimensional  stability  derivatives  along  body  axes.  As  in  the  preceding  case,  these 
axes  may  be  considered  as  general  body  axes  that  include  the  stability  and  principal  axes  as  special  cases. 


Lower  case  basic  symbols  are  used  to  designate  the  nondimensional  body-axis  derivatives.  This  introduces  some  additional 
possibilities  for  confusion  and  ambiguity  in  the  notation.  For  example,  m„'  is  a  stability  derivative  and  m  without  any 
subscript  denotes  the  vehicle  mass.  Although  this  is  not  a  desirable  situation,  ultimately  it  is  less  confusing  to  maintain  the 
system  of  notation  and  be  wary  of  the  pitfalls  of  ambiguity  in  the  notation  than  to  revise  the  familiar  and  established 
symbols.  In  this  instance  the  mass  m  is  frequently  incorporated  in  the  parameter  r  ana  thus  the  confusion  is  prevented. 


V 


Several  examples  of  the  nondimensional  derivatives  along  body  axes  are  given  below.  A  more  pearly  complete 
lifting  of  these  derivatives  and  the  notation  used  is  given  in  table  19.  The  equations  relating  stability-axis  de¬ 
rivatives  to  the  derivatives  along  body  axes  are  given  in  this  table.  The  relations  given  in  table  19  also  serve  to 
transform  the  derivatives  based  upon  wind-tunuel  axes  to  stability  axes  (a0  =  0).  ' 


yp 


_  dX  1  _  y  dw  1 

"  da  q„S  “  Aw  da'  q,S 

_  JN_  _1 _ „  _dv  1 

~  dp  q«Sb  “  dp  q,Sb 

* 

_  _dY_ _ 1 _ _J _ r  _J_* 

“dp  q*S  qoS  “  ^  2V0 

dM  1 _ 1 

“  98.  q<>Sc  ~ qose 


Note  that  in  the  above  examples  the  divisor  changes  and  the  linear  velocity  disturbance  variables  w  and  v  are 
converted  to  nondimensional  variables  a  and  p,  respectively. 


NONDIMENSIONAL  DERIVATIVES  -  STABILITY  AXES 

To  many  individuals  the  term  “stability  derivatives”  means  the  nondimensional  derivatives  of  aerodynamic  coefficients 
with  reference  to  stability  axes.  These  are  the  familiar  parameters  CL#,  C,^,  C^,  etc.  that  are  used  in  aircraft  stability 
and  control  analysis.  Lift  and  drag  are  the  Z-  and  X-force  components.  Wind-tunnel  data  are  usually  reduced  to  stability 
axes  and  provide  experimental  values  for  many  of  these  stability  patameters. 


*  The  nondimentional  rotary  derivative*  retain  the  distention  of  time  in  the  caae  of  body  axe*,  irhile  in  the  cate' of  (lability  axe*  the 
nondimensional  rotary  velocities  are  naed,  i.e„  . 
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There  is  much  literature  containing  analytical  and  experimental  investigations  of  stability  derivatives  and  parame¬ 
ters.  Reference  13  includes  methods  and  reference  material  for  evaluating  these  quantities.  Chapter  5  of  reference 
4  contains  a  comprehensive  discussion  of  airplane  stability  derivatives. 

The  nondimensional  stability  derivatives  referred  to  stability  axes  are  listed  in  table  20.  The  general  notation 
used  is  given  in  table  18.  These  derivatives  are  grouped  into  the  longitudinal  stability  derivatives  or  parameters 
and  the  lateral  derivatives. 

Included  in  table  20  are  sketches  of  typical  variations  of  the  stability  parameters  with  Mach  number.  This  infor¬ 
mation  was  adapted  from  reference  16.  Also  included  in  the  tabulation  of  nondimensional  stability  derivatives  a>-e 
some  specialized  parameters  such  as  Cn^  and  C Tin  effects  and  importance  of  these  two  derivatives  are  dis¬ 
cussed  in  reference  18. 


TABLE  18 

SYMBOLS  AND  NOTATION 

STABILITY  DERIVATIVES  AND  RELATED  PARAMETERS 


SYMBOL 

DEFINITION 

b 

wing  span 

C 

_ 

(i)  basic  symbol  for  aerodynamic  force  and  moment  coefficients 

(ii)  aerodynamic  cross-wind  force 

Cl,  Cd 

lift  and  drag  force  coefficients,  respectively,  (stability  axes; 

__  L  p  D 

'-i  —  ' 

qS  qb 

Ct,  Ct,  g* 

longitudinal,  side-force,  and  normal- force  coefficients,  respectively, 
(body  axes) 

r  _  X  r  Y  r  __Z_ 

*  qS  ’  T  qS  qS 

C|,  C»,  C. 

rolling-,  pitching-,  and  yawing-moment  coefficients,  respectively 

„  L  ^  N  * 

Cl  —  cl  ‘  —  -  •  C.  -  1 

qSb  qSc  qSb 

Coa,  Ci.jj,  Cm/f, 

G  ,,  etc. 

nondimensional  stability  derivatives  with  reference  to  stability  axea 
(see  table  20) 

Cix,  CiT 

moment-of-iaertia  coefficients  and  product -of-inerti*  coefficient 

c.*,  CiXI 

clx  =  -I§-  c,t=JL 

q-Sb  q.Sc 

c,*-'rfb’  c,“  ~  iSb 

Note:  1.  The  divisor  of  CtT  contains  F  instead  of  b.  ^ 

— 

2.  The  inertia  parameters  must  correspond  to  the  axis  system 
used  in  a  particular  analysis,  i.e.,  body  stability,  or  principal 

_ ssa - - 
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TABLE  18  Continued 


SYMBOL 

DEFINITION 

c 

wing  mean  aerodynamic  chord 

D 

aerodynamic  drag  —  the  aerodynamic  force  in  the  plane  of  symmetry 
along  the  projection  of  the  relative  wind  on  the  plane  of  symmetry. 
Drag  is  positive  in  the  negative  X  (downstream)  direction. 

g 

gravitational  acceleration  constant 

H. 

engine  momentum  ,  oooatarolookwtaa  viewed  from  c  ar 

lx,  It,  Is 

moments  of  inertia  about  X-,  Y-,  an  ’  Z-axes,  respectively 

Ik 

product  of  inertia  with  respect  to  X-  and  Z-axes 

Note:  Moment-of-inertia  and  product-of  inertia  terms  must  correspond 
to  the  axis  system  being  used. 

It 

incidence  of  thrust  line  with  respect  to  XZ-plane  of  body  reference 
system.  Thrust  incidence  is  positive  for  Tain  it  acting  in  the  negative 

Z  (lift)  direction. 

L 

aerodynamic  lift  —  the  aerodynamic  force  in  the  plane  of  symmetry 
perpendicular  to  the  projection  of  the  relative  wind  on  the  plane  of 
symmetry.  Lift  is  positive  in  the  negative  Z  (upward)  sense. 

L.M.N 

aerodynamic  rolling-,  pitching-,  and  yawing-moments  about  X-,  Y-,  and 
Z-axes,  respectively 

Note:  Lift  and  rolling  moment  use  the  same  symbol,  L 

L(  >,  Me  i«  N<  j 

basic  symbols  for  dimensional  moment  derivatives:  subscript  denotes 
variable  of  differentiation  (sec  table  1?) 

aL,  am,  an 

incremental  changes  in  aerodynamic  moments  used  in  small-disturbance 
an,’ 

AMt 

,'itching-mnment  component  of  direct  thrust  force 

Note:  When  direct  thrust  is  included  in  the  aerodynamic  or  total 
moment,  AMt  should  be  deleted. 

C  i,  W{  >,  n<  i 

basic  symbols  (or  nond  intensions!  moment  derivatives  about  body 
sues;  subscript  denotes  variable  of  differentiation  (ace  table  19} 

!•>  bh.  o§ 

direction  cosines  between  body  axes  ana  the  gravity  vector 

w 

Mach  number 

Note:  M  is  alto  used  as  pitching  moment 

Note:  as  M  also  u#ed  as  tbs  basic  symbol  for  the  nondmenakmal 

pitching -moment  derivatives  sboci  bod.f  sxwi. 

t,Q,  R 

robing,  pitching,  a ad  yawing  velocity  component*  (angular)  about 

^  .falxiriJiMiksamatoii - 
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TABLE  18  Continued 


SYMBOL 

DEFINITION 

p>  q.r 

small-disturbance  angular  velocity  components  about  X-,  Y-.  and  Z- 
sxes,  respectively 

q 

dynamic  pressure 

*  2 

Note:  <!  is  also  used  as  the  small-disturbance  pitching  velocity. 

S 

wing  area  or  reference  area  for  aerodynamic  coefficients 

T 

net  direct  thrust  force 

Tc 

T 

thrust  coefficient,  Tj  =  —  — 
qb 

TV 

3T 

d\ 

t 

time 

U,  V,  w 

linear  velocity  components  aloug  X-,  Y-,  and  Z-axea,  respect)  ely 

u,  V,  W 

small  disturbance  linear  velocity  components  along  X-,  Y-,  and  Zaxes, 
respectively 

u' 

small-disturbance  noodimetisionai  longitudinal  velocity  variable,  *■' 

*  *  ( 

t 

V 

total  linear  velocity  of  vehicle  c.g, 

X.  Y,  Z 

aerodynamic  force  components  along  X-,  Y-,  and  Z-axeSy  respectively 

X,  ,.Y,  „Z<  , 

basic  symbols  for  dimensional  force  derivative;  subscript  denotes  vari¬ 
able  of  differcntiaiioo  (see  table  17),  For  example. 

Y  ax  Y  ax  v  ax 

du  at.  3W 

AX,  AY,  AZ 

incremental  changes  in  aerodynamic  force  components  used  in  small 
disturbance  analyua 

Xt.  Zt 

components  of  direct  thrust  force  along  X-  and  Z-uea,  reaped ively 

Note:  When  direct  thrust  ts  included  in  the  aerodynamic  force  ooro 

ponmts,  these  terms  should  be  deleted. 

*i  i.  Ti  i.*<  > 

basic  symbols  for  nostdunesuiooal  force  derivatives  along  body  axea; 
subscript  denotes  the  variable  of  differentiation  (see  table  19) 

»r 

distance  parallel  to  Z-axia  from  vehicle  c.g.  to  tbe  projection  of  the 
thrust  line  in  the  plane  of  symmetry  (positive  for  c~g.  above  throat 
line) 

a 

angle  of  attack  (see  figure  22) 
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TABLE  18  Continued 


SYMBOL 

DEFINITION 

a' 

small-disturbance  velocity  variable,  a'  — 

▼  O 

ar 

angle  of  attack  of  thrust  line,  at  =  a  -|-  St 

p 

sideslip  angle  (see  figure  22) 

P 

small-disturbance  velocity  variable,  /9'  =  -2L 

Vrt 

7 

<  flight-path  angle,  the  angle  between  the  velocity  vector  and  the  plane 

of  the  horizon 

i*,  &r 

change  in  deflection  of  ailerons,  elevator,  and  rudder,  respectively 

p 

air  density 

T 

m  V  m 

time  parameter,  r  =  — — — - ZZ — 2-S- 
pSVo  2  q„S 

*,  e,  * 

orientation  angles  of  vehicle  body  axes  in  yaw,  pitch,  roll  sequence 
(see  page  11) 

Note  :  In  some  special  cases  using  stability  axes  the  flight-path  angle 
■y  is  used  in  place  of  the  pitch  orientation  angle  0.  When 
this  is  done,  the  yaw  and  pitch  angles,  +  and  <J>,  should  be 
used  as  referring  to  stability  axes  also  (see  page  77). 

0.  9,  i 

perturbations  of  vehicle  axes  orientation  angles  *,  0,  4',  respectively. 

In  the  small-disturbance  approximation  0  =  /pdt,  8  —  / qdt,  f  = 
/rdt,  respectively. 

0X,  Oy,  Ojr 

angular  velocities  of  wind  axes;  prime  is  used  to  denote  small-angle 
apr.ro*:'*-.ations  ( see  page  77). 

General  Notes: 

1.  All  angles  and  angular  velocities  are  in  radian  measure. 

2.  Fundamental  units  are  used  throughout,  i.e.,  slugs,  feet,  seconds. 

3.  Throughout  this  Uble  the  symbol  q  denotes  dynamic  pressure  when  multiplied  by 

the  wing  area  (qS). 

4.  The  subscript  o  denotes  steady-state  reference  condition  for  small-disturbance 
analyses. 
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>  TABLE  19 


NONDIMENSIONAL  STABILITY  DERIVATIVES 
’  (BODY  AXES) 

Direct  thrust  terms  included 


Symbol 

Derivative 

Expanded  Form  in  Terms  of  Stability 

Ax  is  Derivatives 

Force 

Velocity 

X„' 

JX  1 
du'  q„S 

Cx„  cos*  a„  4  Cza  sin*  a,, 

4"  ( — Cxa  —  C*n)  cos  a*  sin  a» 

2t 

4 - T  y  c  os  ir  cos  Oo 

m 

x,' 

ax  i 

da'  q0S 

Cxa  cos*  ac  —  Cz,j  sin*  a* 

+  (Cr„  —  Cic)  cos  a*  sin  a* 

2r 

4-  _  Tv  sin  iT  sin  t. 
m 

y>’ 

CtP 

w 

C*„  os*  a.  —  Cxa  sin*  a. 

4-  i  Cxa  ~  C*a)  cos  a»  sin  «. 

2r 

—  —  Tv  sin  iT  cos  a> 
m 

u< 

_az_  _i_ 

dc’  q«S 

Cxa  cos*  a.  4-  Cxa  sin*  a. 

4-  (Cx,  +  Ca,)  cos  a.  sin  a. 

2r  m 

—  —  Tv  un  train  a. 
m 

Force 

Rotary 

ax  l 

aq  q»S 

( Cx„  cos  «.  -  Cs,  sin  «.) 

y» 

(CTp  cos  a.  —  Ct,  sin  a*)  ~- 

7t 

(Crr  cos  «»  4-  Crp  sin  «•) 

AT  * 

az  l 

dq  q»S 

{C*s  cos  <u  4-  Cx,  sin  *.) 

Force 

Control 

-42L  J_ 

d«e  q*S 

Cxet  cos  a*  —  Cs£(  tin 

_JY_  _L 

d3t  q»S 

CxK 

y*r 

c*«r 

«*. 

Css4  cos  if  4*  sin  *• 

i 

TABLE  19  Continued 
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NONDIMENSIONAL  STABILITY  DERIVATIVES 
(BODY  AXES) 


TABLE  19  Consumed 

NGNDIMENSIONAL  STABILITY  DERIVATIVES 
(BODY  AXES) 


Symbol  I  Derivative 


dM 


da 


Expanded  Form  in  Toms  at  Stability 
Axes  Derivatives 


(C-i)  ~2v~ 


1.  The  symbol  q.  is  the  reference  dynamic  pressure. 

2.  The  symbol  q  (without  subscript)  denotes  pitching  velocity  about  the  Y-axis. 

3.  The  stability  axes  derivatives  are  defined  in  table  20 

4.  The  subscript  zero  denotes  a  steady-state  reference  flight  condition. 

5.  Symbols  and  notation  are  given  in  table  18  _ 


TABLE  2 n 

NONDIMENSIONAL  STABILITY  DERIVATIVES 
STABILITY  AXES 

Direct  thrust  effect;  not  included 


I.  LONGITUDINAL  DERIVATIVES 


Symbol 


Derivative 


Typical  Variation  with  Mach  No.  (ref.  16) 


Remarks 


C«,  =  -  ~  M.-2(Co). 

Speed  damping  derivative 


subject  to  aarodastic  effects 


Cx*  sr  (C*).  — 


TABLE  20  Continued 


Symbol 

Derivative 

Typical  Variation  with  Mach  No. 

Remarks 

c»« 

flCi 

da 

-C* 

i 

T  i 

®  i  Mo  a 

C‘a=“  ["toT  +  (Cb)*] 

approximately  the  negative  of  the 
lift-cnrve  slope 

dCm 

da 

c, 

a 

c 

■ .  f 

»  1  M0  2 

static  stability  derivative  that 
fixes  the  stick-fixed  neutral  point; 
this  derivative  is  a  basic  static- 
stability  parameter 

<*) 

*°-i 

i 

*  c 

- - ■ 

'  ^  Mo  * 

CM£  is  important  in  damping  of 
short-period  mode;  this  parame¬ 
ter  is  subject  to  high-speed  aero- 
elastic  effects 

Cx, 

dCx 

/  / 

usually  not  significant 

c«. 

dCi 

•(*) 

1 

usually  neglected;  however,  aero- 
elastic  effect*  may  become  signi¬ 
ficant  at  high  speeds 

dC_ 

< 

pitch-dampiBg  derivative;  thia 
parameter  is  significant  in  the 
short-period  mode 

<*) 

»  i  m0  a 

j€x 

W. 


usually  negligible 


Symbol 


Deri  ra tire 


TABLE  20  Continued 
Typical  Variation  with  Mach  No. 


Remark* 


TABLE  20  Continued 


Symbol 

Ddntbt 

Typical  Variation  with  Mach  No. 

Remarks 

ICr 

usually  negligible 

•(#) 

€•# 

ac, 

;  V 

Cifi 

( 

dihedral  effect  — —  C,p  is  im¬ 
portant  in  lateral-dynamic-stabil¬ 
ity  analysis  ^ 

Mo  1  2 

* 

J&. 

9fi 

°v 

1 

weathercock  static  stability  par¬ 
ameter;  important  effects  upon 
lateral  dynamics 

»  Mo  1  a 

c»i 

9£r  ■ 

9fi 

•c» 

fi 

( 

side-force  damping  derivative  — 

CT0  contributes  to  damping  of 
Dutch  roll 

'  M  ‘  * 

mo 

c*# 

ac, 

special  derivative  that  is  signifi¬ 
cant  at  high  angles  of  attack  on 
higl  ly  swept  and  delta-type  wings 
(reference  10) 

•(*) 

S 

ac. 

see  note  above  for  Cij 

w 

J£l 

aa. 

-C„. 

< 

— ^ 

aileron  effectiveness  —  important 
factor  in  establishing  maximum 
rate  of  roll 

f  T 

>  1  2 
*• 
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TABLE  20  Cow  finned 


Symbol 

Derivative 

Typical  Variation  with  Mach  No. 

Remarks 

c*K 

as. 

-C**a 

i 

adverse-yaw  derivative  —  an  im¬ 
portant  item  in  lateral-directional 
control 

m0  »  i 

Cr«. 

3Cr 

as. 

almost  always  negligible 

C,«, 

ac, 

as. 

C,«e 

« 

usually  has  a  small  bat  signifi¬ 
cant  effect  upon  control  and 
'vnamic-atability  analyses 

M  1  i 

"o 

C 

_d£a. 

as. 

C*v 

1 

rudder  effectiveness  —  important 
to  lateral-directional  control 

ig  i  a 

mo 

Ct«r 

acT 

as. 

Ct«r 

i 

usually  negligible  in  dynamic 
analyses 

w  1  s 

"o 

Notea: 

1.  The  subscript  o  denotes  a  steady-state  reference  Sight  condition. 

2.  The  typical  variation  with  Mach  number  is  adapted  from  reference  16. 

3.  Methods  of  evaluating  stability  derivatives  are  given  in  reference  13. 

4.  Symbols  and  notation  an  gives  in  table  18. 

5.  See  reference  18  for  additional  discussion  of  typical  variations  of  stability  derivatives  with 

Mach  numbers. 
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SECTION  5.  SIMPLIFICATION  OF  THE  EQUATIONS  OF  MOTION 


The  general  equations  of  motion  may  usually  be  simplified  for  many  cases  of  practical  interest.  Certain  terms  become  | 
negligibly  small  or  reduce  to  zero  as  a  result  of  practical  considerations  and  selection  of  appropriate  reference  axes,  j 
The  methods  used  to  simplify  these  equations  are  outlined  in  the  paragraphs  that  follow.  ? 

r 

Equations  describing  the  motion  of  a  rigid  body  are  given  in  Section  2.  In  the  following  pages  the  real  force  and 
moment  expressions  from  Section  4  are  used  to  form  equations  of  motion  for  an  aircraft  operating  in  the  atmos¬ 
phere.  Vehicle  symmetry  and  the  assumption  of  small  disturbances  from  reference  flight  conditions  are  then  used 
to  reduce  the  equations  to  simpler  forms.  The  restricted  equations  of  motion  are  given  in  nondimensional 
form  also. 

Steady  equilibrium  flight  and  linearization  based  upon  steady  and  maneuvering  initial  flight  conditions  are  discussed  as 
special  cases.  - 

GENERAL  SIMPLIFICATION  OF  RIGID-BODY  EQUATIONS 

Complete  equations  of  motion  for  a  body  moving  in  the  atmosphere  are  quite  complex.  Consequently,  it  is  of  practical 
interest  to  simplify  them  in  order  to  facilitate  analysis  of  the  motion  of  a  body.  Vehicle  symmetry  and  restriction  of  the 
motion  to  small  disturbances  from  a  reference  flight  condition  are  used  to  reduce  certain  terms  to  zero  and  to  linearize 
the  equations. 


Initial  flight  conditions  are  referred  to  frequently  in  subsequent  paragraphs.  Terms  pertaining  to  these  conditions  are 
defined  below,  as  given  in  reference  19. 

STEADY  FLIGHT  —  Motion  with  zero  rates  of  change  of  the  linear  and  angular  velocity  components,  i.e., 

U  =  V  =  W  =  P  =  Q^R  =  0. 

Steady  sideslips,  level  turns,  and  helical  turns  are  possible  steady  flight  conditions.  Steady  pitching  flight  is  a 

•  * 

“quasi-steady”  condition  because  U  and  W  cannot  both  be  zero  for  an  appreciable  time  if  Q  is  not 
equal  to  zero. 

STRAIGHT  FLIGHT  —  Motion  with  zero  angular  velocity  components,  P,  Q,  and  R  =  0. 

Steady  sideslips  and  dives  or  climbs  without  longitudinal  acceleration  are  straight-flight  conditions. 

SYMMETRIC  FLIGHT— •  Motion  in  which  the  vehicle  plane  of  symmetry  remains  fixed  in  space  throughout  the 
maneuver. 

The  asymmetric  variables  P,  R,  V,  ♦,  and  ♦  are  all  zero  in  symmetric  flight.  Some  symmetric  flight  con¬ 
ditions  are  wings-level  dives,  climbs,  and  pullups  with  no  sideslip. 

ASYMMETRIC  FLIGHT  —  Motion  in  which  any  or  all  of  the  above  asymmetric  variables  may  have  non-zero  values. 
Sideslips,  rolls,  and  turns  are  typical  asymmetric  flight  conditions.  [ 

The  full  set  of  equations  for  the  motion  of  a  rigid  body  is  given  below.  These  equations  are  “Eulerian”  in  that  they  are  ' 
referred  to  axes  fixed  on  the  body.  Because  the  coordinate  axes  rotate  with  the  body,  the  gravity-vector  components  ; 
depend  upon  the  orientation  of  the  body  with  respect  to  a  fixed  inertial  reference  (Earth  Axes).  Relations  are  thus  j 
included  to  express  the  kinematic  angular-velocity  component  relations  in  terms  of  orientation-angle  rates  of  change. 

| 
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The  equations  of  motion  with  reference  to  general*  body  axes  ((8)  and  (9))  are  combined  with  tbe  real  force  and 
moment  components  (equations  (32)  and  (33)).  The  gravity -force  components  along  body  axes  are  then  obtained 
from  the  set  of  equations  (34).  Angular  velocity  relati.  is  are  found  in  table  16.  Collecting  and  combining  these 
relations  give  the  following  sets  of  aynamic  and  kinematic  equations  for  the  motion  of  an  arbitrary  rigid  body. 


X  -j-  XT  —  mg  sin  0  =  m  (U  —  RV  -f-  QW) 

Y  +  Yt  4*  “g  cos  0  sin  *  —  m  (V  —  PW)  +  RU 
Z+  Zr  +  mg  cos  0  cos  #  =  m  (W  —  QU  +  PV)  J 

L  +  Lr  =  PIx  —  QIxy  —  RIzx  —  QR  (It  —  I*)  —  PQIzx 

—  (Q2  —  R2)  I«  +  RPIrx 

M  +  Mt  =  +  QIy  -  RIyz  -PIxt  ~  RP  ( Is  -  lx )  -  QRIxy 

—  (R*  —  P2)  Iz,-fPQI2T 

N  -f-  W*  =4"  ftl*  — PIzx  —  QIy* —  PQ  (lx  —  Iy)  — RPIys 

—  (P2  —  Q*)  I*  y  4"  QRlxe 

4>  —  P  4-  Q  sin  tsn  8  -}-  R  cos  ♦  tan  0 

0  =  Q  cos  4>  —  R  sin  * 

• 

*  =  Qsin4>sec0-f-R  cos  *  sec  0 


(38) 


(39) 


(40) 


Hie  moment  equations  (39)  become  significantly  simpler  when  consideration  is  limited  to  bodies  having  symmetry 
about  the  XZ-plane.  As  a  consequence  of  this  symmetry,  the  prodact-of-  inertia  terms  IXYand  Its  ore  zero.  The 
thrust  components  YT,  Lx,  and  NT  are  zero  except  for  special  asymmetric-thrust  conditions.  The  dynamic  equa¬ 
tions  for  a  symmetric  body  are  then  the  following  (references  3  and  4): 


X  -f  XT  —  mg  sin  0  =  m  (U  —  RV  -j-  QW) 

Y  -f-  mg  cos  0  sin  <t>  =  m  (V  —  PW  4“  RU  ) 

Z  4-  Zr  4"  ®ig  <*>•  9  *  —  m  (W  —  QU  4~  PV) 


(41) 


L  =  PIX  -  (R  4-  PQ)  Ixi  -  QR  (Iy  -  Is) 

M  -f-  Mt  =  QIt-  (R*-P*)  Irs-RP(Is-Ix)  ’  (42) 

N  =  RIZ  -  (P  -  QR)  I**  -  PQ  (lx  -  Iy) 

(The  angular-velocity  relations  are  unchanged  from  (40 ) ). 

Equations  of  motion  in  the  form  above  may  be  modified  to  use  direction  cosines  of  the  gravity  vector  instead  of  the 
orientation  angles.  The  functions  of  the  angles  <t>,  0,  and  *  are  replaced  by  the  direction  cosine*.  {,,  m»,  n,  of  the 
gravity  vector  from  the  vehicle  body  axes.  Relations  between  the  direction  cosines  and  tbe  orientation  «ngW  are  listed 
below  (see  page  24). 

•Tbe  present  dktnmlon  U  chres  is  terns  ef  general  body  axes.  However,  general  body  exee  may  be  interpreted  cs  any  of  the  uii  ayatenw 
fixed  to  tbe  vehicle,  Le.,  body  axaa,  aubility  axea,  or  principal  axee. 
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fg  =  cos  (mg,  X)  =  —  sin  © 
in g  —  cos  (mg,  Y)  —  sin  4>  cos  0 
n8  =  cos  (mg,  Z)  =  cos  <&  cos  ® 


(43) 


Substitution  of  the  direction  cosines  into  equations  (40)  and  (41)  results  in  an  alternate  set  of  equations  of  motion 
for  a  symmetric  body.  These  equations,  which  can  be  found  in  reference  3,  are  especially  convenient  for  use  with 
an  analog  computer.  Equations  of  motion  using  the  direction  cosine  form  are  as  follows: 


X4-  Xt  4-  n'g =  m  (U  —  RV  4-  QW) 

Y  4-  mg  m8  —  m  (V  —  PW4-  RU  )  • 

Z+  Zr  +  mgn8  =m.(W-QU4-PV) 

L  =  PI*  -  (R  +  PQ)  IM  -  QR  (IT  -Iz) 

M  4-  Mt  =  QIt  —  (Ra  —  P3)  I«  —  RP(Ib  —  Ix) 
N  =  RIZ  —  (P  —  QR)  IM  — PQ(IX— It) 

• 

f,  =  m*R  —  n8Q 
m*  =  n8P  —  /8  R  ► 
n,  =  fcQ  —  mgP 


(44) 


(45) 


(46) 


Nondimensional  forms  of  the  foregoing  equations  of  motion  may  be  obtained  by  simply  dividing  through  by  an  ap¬ 
propriate  divisor.  For  the  force  equations  (41)  and  (44)  the  divisor  is  the  reference  dynamic  pressure  times  the 


reference  area 


V2 

O 


*  V  2  V  2 

The  moment  equations  (42)  and  (45)  are  divided  by  !Ll2-§b  or  .  The  latter  value 

2  2 


is  used  only  with  the  equation  for  pitching  motion.  Further  use  of  non-dimensional  equations  occurs  in  subsequent 
Sections,  after  small-disturbance  approximations  are  introduced. 


Many  problems  of  aircraft  motion  involve  only  small  disturbances  from  a  steady  reference  flight  condition.  Thus  the 
approximations  compatible  with  restriction  of  the  motion  to  small  disturbances  allow  further  simplification  of  the 
symmetric-body  equations  (40),  (41),  and  (42). 


General  notation  for  small-disturbance  analysis  is  as  follows.  Perturbations  of  velocity  and  orientation  variables  are 
designated  by  the  lower  case  symbols  for  these  quantities,  i.e.,  u,  v,  w,  p,  q,  r,  <f>,  4,  and  \p.  Upper  case  symbols  are  used 
with  a  subscript  sero  to  denote  the  reference  values  of  these  variables.  Thus  Uw  W<*  Q*  0M  etc.  are  reference  or  initial 
values  for  velocity  components  and  orientation  angler.  Incremental  changes  in  aerodynamic  force  and  moment  com- 
ponenta  are  denoted  by  the  pertinent  symbol  with  a  prefix  A,  e.g.,  aX,  aZ,  AM,  etc. 


Expansion  of  the  aerodynamic  force  and  moment  components  for  the  small-disturbance  approximation  is  discussed  in 
Section  4  and  summarized  in  table  17. 

In  addition  to  the  perturbation  quantities,  the  approximations  noted  below  are  used  in  the  trigonometric  relations  used 
with  the  condition  of  small  disturbances. 


•in  (0.  4-  0)  —  sin  0,  coa  8  -f  coa  0O  sin  9  ' 

5S  sin  0O  -f-  3  cos  0* 

► 

cos  (0,  4*  C)  =  cos  ©o  coa  9  —  sin  0,  sin  $ 

~  coa  0.  —  0  «n  0. 

Note:  These  relations  are  typical  and  are  applicable  to  small -dkturbanoe  approximations  of  any  regular  variable. 
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If  small-disturbance  notation  and  the  above  approximation  for  trigonometric  functions  in  equations  (40)  and  (42)  are 
used,  the  dynamic  equations  for  small-disturbance  motion  expand  to  the  set  of  equations  below.  Products  of  per¬ 
turbation  quantities  are  neglected. 

Xo  -f-  AX  -f-  XT  —  mg  (sin  @0  -j-  $  cos  0O)  —  m  (U„  4  u  —  RoV,  —  R«v  —  V„r  \ 

4-  Q»W»  +  Qqw  +  W„ql  I 

•  •  j 

Y0  4  AY  4  mg  (cos  €>„  —  0  sin  (sin  <!>„  -j-  <j>  cos  <t>„ )  =m(Vu-f  v  +  R,U„  1 

4  R„u  U„r  —  P0W0  —  P„w  —  W0pl  /  (48) 

Zo  4  4  4  mg  (cos  —  6  sin  @0)  (cos  4>0  —  <£  sin  <t>„)  =  m  ( Wn  4  w  1 

-  QoUo  -  QoU  -  U0q  4  PoV.  4P0v  4  V*)  / 

Lo  4  4L  =  (P0  4  p(  I.x  —  (Ro  4  r  4  PoQo  4  P<>q  4  Qop)  I.xz  ^ 

—  (QA  4  Qof  4  R<>q)  (It  —  W  J 

Mu  4  4M  4  aMt  —  (Q0  4  q)  ly  —  (Ro2  4  2  R<>r  —  P02  —  2  P«p)  lx*  ( 

-  (P^Ro  4  PoT  4  RoP)  (Is  -  lx)  I  (49) 

N0  4  AN  —  (R,4r)  Iz  —  (P»4p  —  QA  —  QoX  —  Roq)  Ixs  I 

—  (PoQo  4  Poq  4  Q«p)  (lx  —  It)  / 


Many  of  the  terms  in  the  above  equations  are  zero  for  initial  conditions  of  steady,  straight,  and/or  symmetric 
flight  Linearization  of  these  equations  for  straight,  symmetric  flight  and  maneuvering  flight  is  presented  in  the 
following  paragraphs. 


SIMPLIFICATION  OF  EQUATIONS  OF  MOTION  FOR  STEADY-FLIGHT  CONDITIONS 

Steady-flight  conditions  provide  the  reference  values  for  many  analyses  of  vehicle  motion.  Hie  terms  used  to  de¬ 
scribe  various  flight  conditions  are  defined  on  page  64.  Hie  equations  of  motion  are  reduced  for  several  steady- 
flight  conditions  in  the  following  paragraphs.  Hiese  relations  for  steady  flight  are  used  subsequently  to  eliminate 
initial  forces  and  moments  from  the  equations  of  motion. 


STEADY,  STRAIGHT  FLIGHT 


This  is  the  simplest  case  of  steady  flight.  All  time  derivatives  are  tero  and  there  is  no  angular  velocity  of  the  body  about 
its  center  of  gravity.  Thus,  setting  all  of  the  time  derivatives  and  the  angular-velocity  components  P,  Q,  and  H  equal  to 
zero  in  equations  (41)  and  (42)  results  in  the  following  equations  for  stead}’,  straight  flight: 

X  4  XT  —  mg  sin  8  =  0 

v  i  Q  .  _  n  (50) 

Y  4  mg  cos  9  sin  *  =  0 

Z  4  Zr  4  mg  cos  8  cos  ♦  =  0 


L  =  0 

M  4  Mt  =  0 
N  —  0 


Note  that  these  equations  are  .applicable  to  the  rteady  sideslip.  The  velocity  components  Y  and  W  and  the  bank  angle 
are  not  necessarily  aero.  When  the  motion  it  restricted  to  symmetric  flight,  the  bank  angle  is  aero.  The  force  equations 
for  steady,  straighMyinmetric  flight  ate  then 
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X4-  XT  —  mg  rin  €r  =  0 

Y  =  0 

Z  +  Z?  -j-  mg  cos  0  =  0 


(52) 


The  momects  are  again  a!!  zero  (51).  I 

STEADY  TURNS  } 

* 

f 

In  the  cose  of  steady  taming  flight  the  dotted  quantities  in  equations  (41)  and  (42)  are^  zero,  as  in  the  preceding 
case. 

Also  the  orientation  aogle  rates  of  change  €  and  *  are  zero  and  the  rate  of  turn  *  is  constant  With  these 
conditions  applied  to  the  dynamic  and  kinematic  equations  for  rigid-body  motions,  the  relations  for  steady  turning  fiignt 
may  be  developed.  However,  in  ma.  cases  it  is  convenient  and  practical  to  consider  only  small  elevation  ang.es  or 
shallow  climbing  and  diving  turns.  : 

Applying  the  above  conditions  to  the  anguiar-velocity  relations  in  table  12  results  in  the  following  angular- velocity 
components  for  a  steady  turning  maneuver.  The  approximation  for  small  elevation  angle  (  &  )  is  indicated. 

P  =  —  4  sin  0  SS  —  *  0 

•  m 

Q  =  ♦  sin  6  cos  0  ~  ♦  sin  $ 

a  • 

R  =  ♦  cos  $  cos  0  ~  ¥  cos  $ 

For  most  cases  of  interest  ¥  may  be  considered  as  a  small  quantity,  so  that  the  products  of  the  angular  velocity 
components  P,  Q,  and  may  be  neglected.  In  addition,  for  coordinated  shallow  turns,  the  side  force  Y  is  zero 
and  the  velocity  corr  onents  V  and  W  are  small.  Hie  equations  for  a  steady,  coordinated,  shallow  turn  become 
(see  reference  4) 

X  -f-  Xt  —  m| 
mg  sin  4>  =  m  i 
Z+  Zr  -f*  mg 

L  =  0 

M  -  Mt  =  0 
N  =  0 

Solution  of  the  second  relation  of  (54)  for  the  rate  of 

* 


STEADY  PITCHING  FLIGHT 

Symmetric  flight  of  an  aircraft  along  a  curved  flight  path  with  a  constant  pitching  velocity  Q  results  in  a  quasi-steady  1 
flight  condition.  The  linear  veloci  /  components  U  and  W  must  necessarily  vary  with  time  in  this  crae.  Thus  with 
the  asymmetric  velocity  components  V,  P,  and  R  and  the  bank  and  yaw  angles  ♦  and  *  all  equal  to  aero,  the  equations  j 
of  motion  for  a  symmetric  body  (41)  and  (42)  reduce  to  the  following:  f 


t 
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(57) 


X-j-  XT  —  mg  sin©  =  m  (U -f  QW) 
Y  —  0 

Z  -j-  Zt  -f  mg  cos  8  =  m  (W  —  QU) 

L  =  M  -f-  Mt  =  N  =  0 


The  above  relations  may  be  used  to  evaluate  initial  conditions  for  a  small-disturbance  analysis.  The  values  of  U.  and  W0 
may  be  Hcen  as  instantaneous  values  end  the  variation  with  time  as  disturbance  quantities  u  and  w,  respectively.  For 
reasonable  values  of  pitching  velocity  the  linear  accelerations  u  and  w  may  be  neglected,  so  that  the  X  and  Z  relations 
above  become  initial  conditions 


(X  -f  XT) o  —  mg  sin  0O  —  m  Q0Wo 
(Z  -f  Zt)o  -f  mg  cos  @0  —  — m  Q„U0 


(58) 


Solution  of  the  s^coi-d  equation  above  provides  a  relation  between  the  initial  pitching  velocity  Q0  and  the  initial  load 
factor  %o  aio.ig  the  reference  Z-axis: 


Qo  = 


_S_  l£_±-7sj.<>- 

U0  \  mg 


(59) 


STEADY  ROLLING  OR  SPINNING  FLIGHT 

In  the  preceding  examples  the  steady-flight  equations  readily  reduce  to  simple  forms  of  the  equations  of  motion.  However, 
the  equations  for  steady  rolls  cr  spins  cannot  be  simplified  without  considerable  oversimplification  of  the  physical  relation 
describing  the  motion. 

The  procedure  outlined  in  references  20  and  21  utilises  only  the  mon.snt  equations  to  evaluate  the  perturbed  motion 
from  steady  roll  and  spin,  respectively.  In  such  car »  the  steady  condition  then  becomes  that  of  moment  equilibrium 

(L=(M+  MT)=N  =  0). 


LINEARIZATION  FOR  STEADY,  STRAIGHT,  SYMMETRIC  INITIAL  FLIGHT 


Simplification  of  the  equations  of  motion  for  small  disturbances  from  a  reference  steady-flight  condition  results  in 
the  sets  of  equations  (48)  and  (49).  These  equations,  when  combined  with  the  expansion  of  aerodynamic  force  and 
moment  components  from  table  17,  form  the  linearized  dynamic  equations  of  rigid-body  motion. 

Iu  the  present  cose  many  terms  of  the  equations  of  motion  are  zero,  and  the  steady,  straight,  symmetric  flight 
equab'ms  (52)  are  used  to  simplify  the  equations  further. 

s  •  _  « 

For  the  steady,  straight,  symmetric  initial  flight  condition  the  quantities  V„  P„  Q„  R*  ♦„  ♦„  U„  V„  sad  W, 
are  all  zero**.  Also,  from  the  steady-flight  equations  (54),  the  initial  momenta  (L*  M,  +  MT#,  N,)  and  side  force 
(Y0)  are  zero.  The  initial  equilibrium  in  the  X  and  Z  directions  is  expressed  by  the  relations 

X*  -f  XT#  —  mg  on  «,  =  0 
and 

Z»-f  Zt0  -f  rug  cos  »,  =  0 


—  {I  q-  7it)  • 

mg 


»*  The  eubecript  *.vo  denote*  initial  oocaitioo. 
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The  ir  itial  velocity  components  U0  and  W„  are  related  to  the  initial  velocity  Y,  and  angle  of  attack  a*  by  the 
equations  below 

U0  =  Vo  COS  do 
Wo  =  Vo  sin  <*0 


Application  of  the  foregoing  conditions  a^d  relations  to  the  small-disturbance  equations,  (48)  and  (49),  reduces 
them  to  the  equations  below.  These  are  the  rigid-body  dynamic  equations  of  motion  for  small  disturbances  from 
steady,  straight,  symmetric  flight  Body  axes  are  used  with  this  form  of  the  equations. 

.  > 

AX  —  mg  cos  ®0  =  m  (u -f- q  V„  sin  ae) 


AY  +  mg  if*  sin  80  -f-  mg  cos  ®o  =  m  (v  r  Vc  cos  oo  —  p  V0  sin  oo) 

AZ  —  mgtfsin0o  =  m(w  —  qVoCosao) 


(60) 


AL  =  plx  —  r  Ixz 

AM  =  q  I* 

*  • 

AN  =  r  Iz  —  p  lx* 


(61) 


The  equations  containing  AX,  AZ,  and  AM  are  commonly  referred  to  as  the  symmetric  or  longitudinal  equations  of 
motion.  The  equations  for  AY,  aL,  and  AM  are  then  the  asymmetric  or  lateral  equations  of  motion.  In  the  above  equa¬ 
tions  the  thrust  contributions  to  the  force  and  moment  component  increments  should  be  included  in  AX,  AZ,  AM,  etc. 

The  expansions  of  the  force  and  moment  component  increments  for  small  disturbances  are  summarized  in  table  17. 
The  expansion  outlined  in  the  table  applies  to  any  orthogonal  reference  axes  fixed  to  the  vehicle;  oowever,  once  ; 
axes  are  established,  the  components  and  derivatives  may  not  be  interchanged  from  one  axis  system  to  another.  ! 
The  transforaiation  relations  necessary  to  change  the  derivatives  to  different  reference  axes  are  given  in  tables 
8  and  9. 

Within  the  restriction  of  small  disturbances,  the  perturbation  angular  velocities  are  given  by  the  following  relations: 


P  =  q  = 


(62) 


Linearised  equations  referred  to  stability  axes  are  readily  obtained  from  the  foregoing  set  of  equations.  Stability  axes 
are  oriented  with  the  velocity  V8  at  the  initial  flight  condition  (see  page  13).  Hence  =  V0  and  W^  =  0. 
The  initial  elevation  angle  of  the  stability  axes  is  the  initial  flight  path  angle  y«.  These  conditions  are  equivalent  to  j 
replacing  0*  by  y,  and  a,  by  0  in  the  Unearned  body-axis  equations.  It  is  important  to  note,  however,  that  all  of  the  ! 
velocity,  force,  and  moment  components  in  the  new  set  of  equations  are  referred  to  stability  axes.  The  linearised  equations  ! 
of  motion  referred  to  stability  axes  are  given  below.  These  equations  are  restricted  to  small  disturbances  from  steady,  | 
straight,  symmetric  (light  \ 


* 

AX  —  mg  6  cos  y,-iaii 

AY  -f  mg  y.  sin  y0  4-  mg  4  cos  y*  =  m  (v  4-  r  Ve) 

AZ  —  mg  9  sin  y«  =  m  (w  —  q  V,) 

•  • 

AL  =  p  lx,  —  r  Ixt, 

AM  =  q  It, 

AN  ~  r  ^  —  p  I  xx. 


(63) 
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(64) 


(AX,  aY, . . AN,  4>,  *,  u,  v,  w,  p,  q,  and  r  are  referred  to  stability  axes  in  these  equation*.) 


Aerodynamic  forces  and  momenta  are  usually  reduced  to  noadimenuonal  coefficient  form.  Hence  it  is  convenient  to 
express  the  foregoing  equations  of  motion  in  nondintensional  form.  Fust  the  farce  and  moment  component  increments, 
AX,  AM,  etc.,  are  expanded  a*  indicated  in  table  17.  The  resalting  equations  are  thob  reduced  to  aoodimenoional 
form  and  tabolated  in  tables  21  and  22.  The  equations  are  regrouped  into  the  longitudinal  and  lateral  equations. 

The  tabular  presentation  of  these  equations  is  taken  from  reference  3.  Coefficients  of  die  small  -disturbance  variable*  are 
arrayed  so  that  the  desired  equation  is  obtained  by  setting  the  sum  of  the  products  of  die  coefficients  and  appropriate 
variables  horixontally  across  each  row  equal  to  aero.  The  appropriate  variable  is  given  at  die  head  of  each  column.  For 
example,  the  force  equation  along  the  body  X-axis  is  the  following: 


d 

Note  that  the  operator  is  included  in  the  coefficient  of  the  variable.  Notation  used  in  the  nondimenslonai  equations  is 
summarized  in  table  18. 


TABLE  21 

NONDIMENSIONAL  EQUATIONS  OF  MOTION  FOR  SMALL  DISTURBANCES 
FROM  STEADY,  STRAIGHT,  SYMMETRIC  FLIGHT 
BODY  AXES 


Longitudinal  Equations 

Disturbance  Variable  Coafficknt* 

n' 

a 

$ 

1 

a 

Force  Equation 

d 

X.' 

(  x.  .  \d  n 

(  2r  **  fit 

\  * 

Along  X-axis 

2r  dt 

2r 

—  fcose, 

\  Tr 

\ 

Force  Aquation 

u’  _  d 

(£+— *}s- 

1 

/  * 

Along  Z-axia 

2r 

2r  dt 

2r 

Moment  Equation 

About  Y-axis 

BW' 

— 

dt 

dt  dt* 

- 

Lateral  Equations 

Disturbance  Variable  CsSjalmta 

ft 

# 

* 

t 

Force  Equation 

Along  Y-axis 

2r  dt 

+  -£-ca*«. 

IBaM 

ESI 

St 

tr 

St 

tr 

Moment  Equation 

About  X-axis 

V+-*sv 

r  d  c  e 

QEOI 

m 

Hi 

Moment  Equation 

About  Zads 

V+4r*r 

Ur-i-Csg  J 

•s 

■A  | 

Symbols  md  notation  am  ghun  k  labia  18  .  4 
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TABLE  22 

NONDIMENSIONAL  EQUATIONS  OF  MOTION  FOR  SMALL  DISTURBANCES 
FROM  STEADY,  STRAIGHT,  SYMMETRIC  FLIGHT 


STABILITY  AXES 

(UMPUTIBD  FORM) 


Longitudinal  Equation* 

Disturbance  Variable  Coefficients 

o' 

tl 

$ 

a 

Force  Equation 

Along  X-axis 

2  (Cs). 

+  *TA 

—  (ClI.  +  Co, 

cvs:t+,w- 

c»* 

Force  Equation 

Along  Z-axis 

2  (Cl). 

Cl#  +  (Co). 

-*>  i 

+  (Cl),  tan  >. 

Moment  Equation 

About  Y-axis 

0 

C-. 

+  C  —  — 

+  lj-A2V.  dl 

• 

C  d  <P 

^  TV.  dt  ~C,T  dt* 

c-a 

Lateral  Equations 

Disturbance  Variable  Coefficients 

* 

f 

a. 

a. 

Force  Equation 

Along  Y-axis 

m 

r  b  d 
»  XV.  dt 

+  «*). 

-4-(G.).tan  y. 

Ct*. 

Cr*r 

Moment  Equation 

About  X-axis 

c,t 

c,»w.  * 

c,  -Li 

r  2V.  dt 

+  G*.  £ 

Cm. 

C-r 

—  lleeaest  Equation 

About  Zautio 

£ 

^  TV.  dt 

-On  i 
dt* 

c*. 

Symbols  an i  netatiau  me  gfeon  in  l able  18  ■ 
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LINEARIZATION  FOR  STEADY  MANEUVERING  FUGHT 


There  are  certain  typea  of  problems  in  analysis  of  vehicle  motion  in  which  the  assumption  of  small  disturbances  from  a 
steady  maneuvering  flight  condition  is  the  most  efficient  method  of  approach.  The  limitation  to  a  steady,  straight 
symmetric  initial  flight  condition  is  unnecessarily  restrictive. 

General  equations  for  amall-diaturbance  motions  from  steady  flight  are  derived  and  summarized  concisely  in  ref¬ 
erence  22.  These  equations  are  equations  of  motion  for  small  disturbances  from  steady  turning,  pitching,  rolling, 
or  longitudinally  accelerating  flight.  Table  23  gives  the  general  dimensional  equations  for  smalbdistnrbance  mo¬ 
tion.  The  next  table  ( 24 )  summarizes  die  conditions  to  be  used  in  the  general  equations  of  table  23  for  die  several 
types  of  steady  initial  flight  conditions. 

These  equations  readily  reduce  to  those  used  in  special  cases  treated  in  the  literature,  e.g.,  references  20  and 
23.  Also,  the  small-disturbance  equations  of  the  preceding  Section  for  steady,  straight,  symmetric  flight  may  be 
obtained  from  table  23.  Stability  axes  are  the  reference  axes  for  this  case  and  the  notation  is  defined  in  table  18. 

TABLE  23 

GENERAL  EQUATIONS  OF  MOTION  FOR  SMALL  DISTURBANCES  FROM  STEADY  FUGHT 

STABILITY  AXES 


Codfidmt  of  Ditfarhuce  Variable 


a 

B 

B 

p 

4 

r 

*+(*: 

).= 

n 

B 

Q. 

0 

BD 

— g0aa7.de  A.  S  <  )  A 

AY 

■ 

n 

H 

V.  —  [  gds T»  + 

m 

H 

u 

■ 

1  OMR  7a  m  M»  )  l  1  RE 

0 

(£).]><  >* 

M 

m 

s 

-Q. 

B 

9 

i 

fcaav.  dsAw  J  (  )  A 

i 

* 

* 

0 

AL 

E 

B 

B 

fl 

h,  -  —  Q*  In, 

R.  (Is,  —  It,)—  P«  In, 

-In.^+  Q.  (L,  It.) 

AM 

s 

0 

0 

0 

R*  (I*,  —  U.)  +1P.  In, 

iT  AL  .1 

P.(l*.-k,)-IR.!„. 

ss 

B 

B 

fl 

_In,iLJ  +  Q.(I,,-IB,) 

P.  (It.  -  Ik.)  +  R.  In. 

Note:  0  ss  J  f  A 


#  =  f«* 

♦  as  |t  A 

Symbols  and  aomtiaa  Mad  aa  taMe  18 
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$*•. '  W  ■  y '  -i V  *  -t;  tv*.;  V*  r1?.  »> 


TABLE  24 


IAL  CONDITIONS  FOB  GENERAL  SMALL-MSTURBANCE  EQUATIONS  OF  MOTION 


untin^i 


OHMy  oCrtlQM 
Spaasstric  Flight 


Stand?  Tareirg  Flight 


STABILITY  AXES 


(Ceadkiaas  for  ate  witk  table  23  » 


Initial  Yalae  of 
Velocity  Components  aod 
Orientation  Ancles 


«.  =  V,;V.  =  W.  =  (f)  =0 
P.  =  Q.  3S  R.  s  0 


Urn  =  ■».;  ♦.  =  0 


U.  =  ».;V.  =  W.=  (f.)  =. 

P.  =  -+.  tin  if. 

Q.  =  ♦.  eaa  y.  m  ♦. 

R.  s  t.cMy.eaa*. 

■y.  =  7.;  ♦.  = 


u.  =  v.;v.  =  w.  =  (tjL)  =. 

P.=  P.;Q.  =  R.  =  0 


if«  —  7»» 


V.  =  ,.  =  (^)  =. 

P.  =  R.  =  Os  Q.  =  (a?  —  oaa  t.) 

7.  =  >.;♦.  =  0 


U.  =  V.;  V.  =  W.  =  0 
P.  s  Q.  =  L  =  0 


(*).- 


Same  aa  presented  on  page  70 
(see  equations  (63)  and  (64)). 


This  Mgla  aaadMaa  is  alw 
stand?,  dans  A  i*  0. 


Ia  addition  to  die  equation*  of  motion  given  in  table  23  them  are  other  specialized  forms  of  these  equations.  The 
study  of  small  disturbances  from  steady  pitching  flight  may  be  conveniently  made  with  vehicle  body  axes  for  ref¬ 
erence.  Below  are  die  equations  for  analysis  of  snail-disturbance  motion  from  ‘steady*  pitching  flight*  as  given 
in  reference  3.  Body  axes  are  used  and  notation  is  defined  in  table  18. 


AX  —  9  mg  co#  80  —  m  (u  -f  q  V0*in  a,  -f  w  Q0) 

AY+  f  ingiin#*.-}- *mgcose,=  m  (v  +  rV.cos«e  — pV.sinaJ 

AZ  —  9  mg  sin  «*  =  m  (w  —  q  V.  coe  «*  —  u  Q,) 

AL  =  I*p  —  Ix*r  +  Q<>r  (I*  —  It)  —  Qj>  lx* 

AM  —  IyO 

AN  =  I*r  —  Imp  -f-  Q«p  (It  —  lx)  +  Q«f  lx* 

Note:  These  equations  may  be  developed  from  (48)  and  (49). 


) 


(65) 


U,  =  V.cos  V.  =  0;  W.  =  V.sin  «, 
sIsoP.IU,#.  =  0 


(66) 


The  steady  pitching  velocity  is  given  by 


In  the  above  relation  other  the  initial  pitching  velocity  Q,  or  the  initial  load  factor  ns  may  be  pacified. 


A  case  of  — disturbance  from  a  steady  fight  condition  oecars  ia  the  analysis  of  dm  spinning  motion  of  an 

aircraft  An  analysis  of  the  dynamics  and  stability  of  flat  spins  is  contained  ia  reference  21.  Bsc  ansa  of  the  spacialiand 
nature  of  the  analysis  and  the  likelihood  of  confution  is  the  symbols  sad  notation,  the  equations  are  not  repeated  hen. 


ADDITIONAL  SKCKAUZCD  POUKS  OP  THE  EQUATIONS  OP  MOTION 

Rigid-body  equations  of  motion  have  been  dsvalopsd  and  specialissH  for  many  pacific  and  special  problems  No  attesnpt 
has  been  mads  to  collect  them  aB  hare.  However,  some  gswaraBy  nsefnl  farms  of  these  equations  aw  summarised  in 
titisSaotfaiL 


•4. 
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LARGE  DISTURBANCES  FROM  STEADY,  STRAIGHT,  SYMMETRIC  FLIGHT.  (References  3  and  19) 


In  the  equations  of  motion  presented  below,  the  approximation  of  small  disturbances  has  been  limited  to  the  linear 
velocity  components  U,  Y,  and  W.  The  remaining  orientation  angles  and  angular-velocity  components  are  not  restricted 
to  small  values  as  in  the  case  of  general  small-disturbance  motion.  Many  practical  problems  in  aircraft  motion  may  be 
analysed  under  these  conditions. 

The  initial  condition  used  is  steady,  straight,  symmetric  flight;  hence  V0(/?0),  POT  Q<»  and  are  zero.  The  small, 
nondimensional  linear  velocity  disturbances  are  denoted  by  u',  /?\  and  a,  and  the  reference  axes  are  the  vehicle  body 
axes.  These  equations  are  developed  from  the  symmetric  rigid-body  equations  (41)  and  (42)  and  put  in  nondimen- 
sional  form.  Special  notation  used  in  these  equations  is  defined  in  table  18. 

(If- -5-)' “'+ + (-t-« ■*)«+ V R - * 


+  3-  «,  —  sin  6  -f  ~S~  sin  0*  =  0 
2r  V* 


76 


oWV ; < ■.  •;  rP-l-.'Ts  S-- 


Note:  1.  Engine  angular  momentum,  H*,  is  included  in  this  set  of  equations. 
Note:  2.  Symbols  and  notation  are  defined  in  table  18 . 

Note  3.  Reference  axes  are  vehicle  body  axes. 


EQUATIONS  OF  MOTION  ALONG  WIND  AXES 


Generally  the  dynamic  equations  of  motion  along  wind  axes  are  too  cumbersome  for  use  in  vehicle-motion  analysis.  The 
variation  of  inertia  parameters  with  orientation  angle*  precludes  tny  extensive  exploitation  of  the  simplified  aerodynamic 
terms  along  wind  axes.  Hie  dynamic  force  equations  and  kinematic  relations  are  sometimes  useful,  however,  and  are 
therefore  given  here  (equations  (68)  and  (69)).  These  are  taken  from  reference  3,  hut  they  may  be  developed  di¬ 
rectly  from  equations  (40)  and  (41).  (Symbols  and  notation  are  defined  in  table  18.) 


T  —  D  —  mg  sin  y  =  m  V 

— T/?  -f  C  -+■  mg  sin  *  cos  y  —  m  V  (/9 -f*  R  —  P*1 

=  m  V  (O'*) 

— T(«+i,)—  L  -f-  mg  cos  ♦  cos  y  —  m  V  («*  —  Q  -f  P/J) 

=  mV(-0'T) 


(68) 


♦  =  O'*  -f  (n't s*n •  -f*  O'* cos#)  tan  y  —  O'*  +  *  sin  y" 
y  —  (fj  cos  ♦  —  O'*  sin  # 

+  =  (O' t  sin  ♦  -f  cos  #)  sec  y 


(69) 


From  the  second  and  third  equations  above,  the  following  relations  for  the  rates  of  change  of  angle  of  attack  and  of 
sideslip  angle  are  obtained. 


«  =  Q-P0-n't 
&  —  P« —  R  +o'i 


(70) 


A  practical  application  of  the  above  equations  occurs  in  the  simplified  analysis  of  inertial  coupling,  e^.,  that  of 
reference  24.  In  this  case  it  was  desired  to  develop  a  simplified  analysis  that  provided  a  quick  and  ample  method  for 
surveying  the  dynamics  of  a  rolling  airrrsft  Problem  areas  could  subsequently  Le  more  thoroughly  and  rigorously 
investigated-  The  simplified  analysis  was  then  made  by  using  the  shove  force  equations  along  wind  axes  and  the  moment 
equations  along  principal  axes  (equation  (13)). 
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SECTION  6.  SOLUTION  OF  THE  EQUATIONS  OF  MOTION 

The  equations  developed  and  presented  in  the  preceding  Sections  describe  the  motion  of  a  particle  mass  and  that  of  a 
rigid  body.  Solutions  of  the  complete  equations  are  not  always  possible  or  may  be  impractical  for  the  problem  under 
consideration.  Several  methods  of  simplifying  these  equations  are  given  in  Section  5. 

Methods  for  solving  differential  equations  found  in  many  standard  mathematics  texts  may  he  applied  to  find  solutions 
of  the  equations  of  motion.  Some  general  methods  for  solving  the  equations  of  motion  are  outlined  in  the  paragraphs 
that  follow.  Included  in  this  Section  is  a  brief  discussion  of  computer  methods  and  of  some  approximate  solutions.  The 
approximation  formulas  are  useful  for  preliminary  estimates  of  dynamic  stability  characteristics. 


ANALYTICAL  METHODS 

Solution  of  the  simplified  equations  of  motion  by  analytical  methods  is  possible  in  many  cases.  The  simplified  equations 
are  generally  a  system  of  ordinary  linear  differential  equations  having  constant  coefficients. 

* 

Use  of  the  direct  method  of  solution  is  outlined  for  the  linearized  small-disturbance  equations.  The  Laplace  transform 
method  is  also  outlined  and  a  matrix  method  noted. 

Analytical  methods  for  solution  of  nonlinear  systems  of  differential  equations  are  not  included.  References  25  and  26 
present  analytical  methods  for  obtaining  solutions  of  the  motion  in  nonlinear  dynamic  systems. 


DIRECT  METHOD  OF  SOLUTION 


The  direct  method  of  solution  for  a  system  of  ordinary  linear  differential  equations,  such  as  the  small-disturbance  equa¬ 
tions  of  motion,  is  described  and  illustrated  in  Chapters  6  and  7  of  reference  A  This  procedure  is  reviewed  below  with 
the  longitudinal  equations  of  motion  as  an  example. 

Equatioua  of  notion  for  small  disturbances  are  separated  into  a  set  of  longitudinal  equations  (symmetric)  and  a 
set  of  lateral  equations.  These  equations  are  given  in  nondimensional  form  in  tablea  21  and  22.  The  longitudinal 
equations  from  table  22  are  used  below  to  illustrate  the  direct  method  of  solution. 


If  the  control  remains  fixed,  I  is  aero  and  the  longitudinal  equations  reduce  to  s  system  of  simultaneous  ordinary  homo¬ 
geneous  differential  equations.  These  equations  then  have  the  dependent  variables  u\  a',  and  9  as  functions  of  time,  the 

independent  variable. 

The  solution  for  the  dependent  variables  is  then  assumed  to  be 


u'  —  u'„e*'  j 
a'  «—  o,„  e*’  \ 

9  ) 


(71) 


Substituting  these  relations  into  the  longitudinal  equations  of  table  22  results  in 


the  set  of  equations  below. 


[2(Cd).  -f  2tA]  +  (CD.  -  (CtM  .;e*‘  +  (C -  0 
2(Cl).  u'e*‘  -j-  [Cl,  4-  (CdL  +  2tA ]  «;e*‘  -  2tA  *e*’  -  0 


\ 

I 
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Note:  CD,  and  Ct^  are  assumed  negligible.  Steady  level  symmetric  flight  is  assumed  (y0  —  0).  Notation  is  given  in  table 
18. 

The  factor  e*‘  is  common  to  all  terms  in  the  above  equations  and  may  be  divided  out*  The  result  is  then  a  set  of  linear 
simultaneous  algebraic  equations  in  the  variables  u'M  oV,  and  0O  with  a  parameter  A  to  be  determined.  The  condition 
upon  A  required  for  nonzero  values  of  the  dependent  variables  is  that  the  determinant  of  coefficients  of  equations 
(72)  be  zero.**  Thus, 

[2{Ca)o  +  2rA]  [CD<S-(CL)0]  (CL). 

2(Cl)0  [CLa+  (CD)o  +  2rAi  —  2tA 

0  +  Cni  -$y-  A  J  [Ca,  2\r  A  -  C«T  *‘] 


Expansion  of  this  determinant  results  ir«  the  characteristic  equation  for  the  solution.  This  is  a  fourth-degree  polynominal 
in  the  parameter  A. 

AA4  +  Ba>  +  Ca*  +  DA  +  E  -  0  (74) 

where: 


A  =*  -4raC, 


8  .  -2  rC,T  (3CDj>  4  C L.)  4  4  C„  4  ^  C.4j 

C  --2C,r(CDoCLti-  CLCDJ  -  2CD/  C,r-2C^>  4  4r 


[Cd^v0 


C-.  t  7C  C- 


4  2tCdo  C%)  4  2tClJ^  <;„)  t  4t,c^ 

d  ■  +  w.^)*  + rjC^(hc-)  -X^[  W0C-,)*  «tC».C.. 


E  .  2C,  *  C., 


The  roots  of  this  equation  are  the  values  of  A  corresponding  to  the  modes  of  motion. 

The  roots  of  (74)  may  be  real  or  complex.  The  complex  roots  necessarily  occur  in  conjugate  pairs  and  denote  an 
oscillating  mode  of  motion.  Each  real  root  corresponds  to  a  pore  convergence  or  divergence  without  any  oscilla¬ 
tion***  Convergence  or  divergence  of  each  mode  of  motion  is  established  by  the  sign  of  the  real  root  or  the  real 
part  of  the  comnlex  conjugate  roots.  Four  types  of  motion  are  possible,  as  illustrated  in  figure  24,  for  the  func¬ 
tion  0  ~  e** 


•The  anlulio*  e*1  =  0  i*  trivial, 

••Not*  lKat  the  eolation  (or  »ny  o(  the  'triable*  in  equalise  (72)  would  knee  the  determinant  in  the  numerator  equal  to  aero.  Thu*  for 
a  variable  to  have  a  nonzero  value  the  denominator  determinant  nu  u  be  aero.  The  resulting  indeterminate  form  may  then  bare  nonzero 
value*.  Alto  it  should  be  noted  that  the  general  solution  o f  the  longitudinal  equations  involve*  a  forcing  function  such  as  a  control  pulse, 

so  that  the  right  side  of  equation  (72)  is  not  ail  aero. 

•••Whenever  the  characteristic  equation  it  of  as  odd  degree,  it  must  have  at  hurt  one  real  root. 
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The  behavior  of  the  dependent  vari abler  u',  a,  and  6  may  be  determined,  once  the  roots  of  the  characteristic 
equation  (values  of  x)  are  known.  The  variation  of  those  quantities  with  time  is  then  given  by  equations  (71) 
for  each  mode  of  motion  of  the  system,  with  additional  constant  multiplying  factors  depending  upon  the  input. 

Is  addition  to  the  time  history  of  each  dependent  variable,  several  quantitative  parameters  that  describe  the  mo¬ 
tion  may  be  determined  from  the  roots  of  the  characteristic  equation.  These  items  are  the  period  and  the  time  to 
halve  (or  doable)  the  initial  amplitude.  The  cycles  to  halve  (or  double)  the  amplitude  of  oscillatory  motion  are 
also  of  interest.  Table  25  lists  these  items  for  both  real  and  complex  conjugate  roots  of  the  characteristic 
equation.  f 


(*)  PURE  DIVERGENCE 

\  *8  BEAL  AND  POSITIVE 


(b)  PURE  CONVERGENCE 

K  18  BEAL  AND  NEGATIVE 


(o)  DIVERGENT  OSCILLATION 
A.  IS  IMAGINARY  WITH 
POSITIVE  REAL  PART*/?  (A) 


HQURE  34  TYPES  OP  MOTION  FOR  DIFFERENT  BOO  TS  OF  THE  CHARACTERISTIC  EQUATION 


SO 


TABLE  25 

QUANTITATIVE  CHARACTERISTICS  OF  MODES  OF  MOTION 


Characteristic  of  the  Motion 

Real  a 

Complex  A 

Time  to  halve  or  double  amplitude 

0.693 

0.693 

w 

wm 

Period 

Not  Applicable 

2w 

-  W51 

Cycles  to  halve  or  double  amplitude 

Net  Applicable 

-Olio  Siili 

0110  K(A)| 

LAPLACE  TRANSFORM  METHOD  OF  SOLUTION 

Use  of  the  Laplace  transformation  in  the  solution  of  linear  differential  equations  has  several  advantages  when  compared 
with  the  direct  method  of  solution.  This  method  of  solution  of  the  small -d isturbance  equations  of  motion  is  explained 
and  illustrated  in  references  4  and  11.  There  are  numerous  texts  that  contain  the  mathematical  development  of  the 
Laplace  transform,  such  as  references  27  and  28.  These  references  also  provide  additional  examples  and  tables  of 
Laplace  transforms. 

The  primary  advantages  of  using  Laplace  transformations  to  obtain  solutions  of  systems  of  linear  differential  equations 

are: 

1.  Initial  conditions  are  introduced  directly  into  the  solution  in  order  to  avoid  the  evaluation  of  the  constants  of 
integration  required  by  direct  methods. 

2.  In  problems  involving  several  dependent  variables  the  solution  for  one  variable  may  be  obtained  independently. 
The  literature  provides  detailed  treatment  of  the  Laplace  transform  method.  Examples  sad  disease  ion  of  this 
method  applied  to  vehicle  motion  analysis  are  found  in  references  4  and  11.  Included  in  these  references  are 
methods  of  presentation  and  interpretation  of  results  in  terms  of  both  the  basic  variables  and  the  transformed 
variables. 

To  illustrate  the  correspondence  between  the  Laplace  transform  method  and  the  direct  method,  the  longitudinal 
small-disturbance  equations  from  table  22  are  transformed  below.  A  bar  is  used  over  die  symbol  to  denote  the 
transformed  variable.  For  example, 

U'f9)-./[u'ft)] 

Applying  the  Laplace  transform  to  the  longitudinal  equations  from  table  22  results  in  the  following  transfonned 
equations.  The  derivatives  Cq,  and  Cl,  are  again  assumed  negligible  and  initial  level  flight  ie  used  (y«  —  0). 
The  control  is  considered  to  be  fixed  so  that  $  *  0. 

[2{C0)o  +  2ra]  u'  +  [CD.  -  (CL)0J  S'  +  2r  o'(0) 

2(CL)0  u'  +  {Cl,  4-  (CD)0  4- 2rs]  S'  -2«?-2r  [*'(0)  -#(0)] 

+ Cmiw; 8  ] « 4-  [°-.  w:*  ~  *  ]f- 

4~  C*,  2^  '  ^(0)  —  Qy  [1(0)  4-  fd(01J 

(Notation  in  the  above  equations  is  defined  in  table  18.) 

u  the  real  part  of  the  complex  number  x  tod  it  related  to  the  modal  at  or  amplitude  of  the  voder  ropreoeatxtioa  of  the  amahar  la 
the  complex  plane. 

0(M  i»  the  imaginary  part  of  the  complex  number  X  end  it  related  to  the  utfular  velocity  of  dm  vecear  reprmenlrlteei  of  the  mmnhw  in 

the  complex  piene. 


•  1 


Note  the  similarity  between  these  equations  and  those  in  the  preceding  Section  (equation  (72)).  The  coefficients 
of  the  transformed  variables  are  die  same  as  those  obtained  with  the  direct  method  when  A. is  replaced  by  s. 

For  an  nndistnrbed  steady-state  condition  the  terms  on  the  right  side  of  equations  (75)  are  all  zero  and  the  trivial 
case  exists.  However,  if  a  disturbance  such  as  a  control  pulse  or  a  gnat  io  int  need,  a  dynamic  motion  problem 
is  generated. 

Suppose  an  aircraft  in  steady  level  flight  encounters  a  gust.  The  tenn  a!  (0*)*  is  then  different  from  zero,  and 
the  solution  of  equations  (75)  for  the  pitch  angle  may  be  expressed  in  determinant  forms.  Thus, 


[2(Cd)o  ~f-  2ts] 

[2(Cl)0] 


[CDa-  (CL)0] 
[CL  -f-  (CD)  o  +  2rs] 


Cra  -f-  Cm 


C 

EvT8J 


[2(Cd)o  +  2ts] 
^[2(Cl)0] 


[CB.-(Ct).] 
[CL  *+■  (CD)o  +  2ts] 


[Cl)0] 

[  2rs] 


Cm  -f-  Cm 


•  2V0  8 


S  —  Ciy  8* 


Expansion  of  the  above  determinants  results  in  a  quotient  of  two  polynomials  in  the  transformed  independent  variables. 
The  denominator  determinant  expands  to  the  characteristic  equation  of  the  direct  method  (see  equation  (74)). 

The  solution  for  the  pitch  angle  6  as  a  function  of  time  requires  application  of  the  inverse  transformation  of  6(a).  Thus, 

*(t)  -  X-'  [*(*)] 

In  order  to  simplify  the  inverse  transformation,  it  is  usually  expedient  to  separate  the  expansion  of  equation  (72) 
into  partial  fractions.  This  procedure  then  requires  finding  the  roots  of  the  denominator  or  characteristic  equation. 

The  zeros  (roots)  of  the  denominator  of  equation  (76)  have  the  same  significance  as  the  roots  of  the  characteristic 
equation  (74)  in  determining  the  modes  of  motion  as  shown  in  figure  24. 

MATRIX  METHOD  OF  SOLUTION 

A  method  of  solving  the  equations  of  motion  using  matrices  is  presented  in  reference  29.  This  is  a  procedure  more 
readily  adapted  to  machine  computation  methods  than  to  analytical  methods. 

Briefly,  the  procedure  consuls  of  a  stepwise  integration  of  the  differential  equations  with  a  Maciaurin  series  expansion 
used  in  each  computation  step  to  achieve  any  desired  degree  of  accuracy.  This  method  is  a  rather  specialised  technique, 
and  the  reader  is  referred  to  the  cited  reference  !cr  the  detailed  explanation  of  the  method. 

COMPUTER  METHODS  OF  SOLUTION 

The  development  of  modern  machine-computing  equipment  has  opened  Use  way  for  many  new  and  varied  analyses  to 
be  undertaken.  Problems  that  are  impractical  to  solve  by  lengthy  hacdcottoputarion  methods  are  readily  computed  by 
high-speed  digital  computers.  Problems  that  involve  nonlinear  equations  may  be  solved  quickly  on  an  analog  computer. 
The  use  of  machine  computation  methods  also  permits  more  variables  (degrees  of  freedom)  to  be  considered  and 
reduces  the  number  of  approximations  or  coemption*  that  must  be  made  in  order  to  facilitate  solution  of  a  problem. 
Machine  computation  method*  thus  provide  a  large  increase  in  thr  amount,  scope,  and  accuracy  of  analysis  possible  in 
many  prob! om¬ 
it  f»  beyoad  the  acope  vf  this  repart  to  present  a  complete  discourse  on  machine-comparing  methods.  The  para¬ 
graphs  that  follow  provide  some  general  background  information  and  references  for  detailed  treatment  of  tbe  subject. 


•All  of  the  taitia!  oouditiMw  «adt  *»«.«»,  a' (01,  ocoKfia*  os  the  right  tide  ef  equations  (T»)  jbsald  be 

at  +  *>  m  e*. 


itstpre—d  •>  she  value 


DIGITAL  COMPUTER 


Electronic  digital  computers  have  been  developed  to  a  very  high  degree  as  fast,  automatic  computing  systems.  The 
digital  computer  is  a  device  that  automatically  performs  the  basic  operations  of  arithmetic.  It  performs  these  operation., 
in  a  sequence  prescribed  by  the  program  for  a  given  problem.  Since  the  digital  computer  functions  as  a  mechanical 
desk-type  calculator,  it  is  capable  of  very  precise  computation  (many  significant  figures).  Any  problem  that  can  be  set 
up  for  hand  computation  can,  in  principle,  be  programed  for  an  automatic  digital  computer. 

Reference  30  offers  a  thorough  presentation  of  the  principles  and  features  of  digital  computers  and  data  proc¬ 
essing.  The  matrix  method  and  the  solution  of  the  characteristic  equation  in  the  preceding  Sections  are  examples 
of  calculations  that  may  readily  be  programed  for  a  digital  computer.  Iterative  processes  and  decision-making 
routines  may  be  incorporated  into  a  digital  program. 

The  digital  computer  can  perform  very  complex  calculation  routines  that  involve  comparison  with  previously  computed 
or  reference  data  and  can  then  choose  a  procedure  according  to  one  of  several  alternate  subroutines.  This  is  accomplished 
very  rapidly  and  very  precisely.  The  digital  computer  is  used  most  advantageously  in  stability  and  control  calculations 
for  making  large  numbers  of  calculations  of  a  given  type,  such  as  the  response  to  arbitrary  control  functions  or  the 
dynamic  behavior  of  a  flight  vehicle  with  a  well-defined  automatic  control  system. 


ANALOG  COMPUTER 

Analog  computers  are  a  combination  of  electrical  and  mechanical  components.  Hiese  components  are  selected  and 
arranged  so  that  the  differential  equations  of  the  analog  system  are  dual  to*  the  differential  equations  for  the  problem 
being  studied.  Electrical  components  are  used  in  most  analog  computers.  Reference  4  gives  a  resume  of  analog  compo¬ 
nents  and  their  basic  function.  References  30,  31,  and  32  are  comprehensive  texts  or  handbooks  covering  the  design  and 
application  of  analog  computer  systems. 

Certain  features  of  analog  computers  are  quite  different  from  those  of  digital  computers.  The  electrical  analog  system 
operates  with  either  the  current  or  the  potential  in  a  component  circuit  representing  a  variable  of  the  problem  being 
studied.  Thus  data  are  continuous  and  all  operations  are  simultaneous,  while  the  digital  computer  must  follow  a  pre¬ 
scribed  sequence  of  operations  (program)  on  distinct  pieces  of  data.  The  accuracy  of  an  analog  computer  depends  upon 
the  precision  and  quality  of  its  components.  Analog  computers  are  generally  less  accurate  than  digital  computers. 

The  analog  computer  has  been  used  extensively  in  airplane  stability  and  control  analysis.  It  is  readily  adapted  to  solving 
the  equations  of  motion  when  nonlinear  characteristics  must  be  included. 

A  very  useful  application  of  analog  computers  is  the  flight  simulator,  since  it  can  calculate  in  real  time.  This  device 
extends  the  analog  simulation  to  include  duplication  of  the  cockpit,  controls,  and  flight  instruments.  Simulators  have 
been  built  tor  many  different  types  of  aircraft  and  used  for  flight  research  and  for  familiarisation  and  training  of  pilots 
and  flight  crews. 


APPROXIMATE  SOLUTIONS 

Frequently  approximate  solutions  of  the  equations  of  motion  are  useful.  Preliminary  estimates  and  quick  evaluation 
of  flying  qualities  often  require  drastic  simplifications  of  an  analysis.  The  paragraphs  that  follow  present  some  useful 
approximate  solutions  of  the  equations  of  vehicle  motion. 


•Of  the  time  form  •». 
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APPROXIMATE  FORMULAS  FOR  SMALL-DISTURBANCE  MOTION 


Relations  resulting  from  approximate  solutions  of  the  equations  of  motion  are  listed  in  table  26.  These  formulas 
are  developed  from  the  small-disturbance  equations  of  motion  along  stability  axes  in  table  22.  Approximate  quan¬ 
titative  values  for  the  characteristics  of  the  normal  modes  of  motion  are  provided  by  these  relations.  The  informa¬ 
tion  in  table  26  is  adopted  from  a  similar  tabulation  in  reference  3  and  utilizes  the  notation  of  table  18. 

TABLE  26 

APPROXIMATE  FORMULAS  FOR  SMALL-DISTURBANCE  MOTION 


I  Low  Frequency 
Root  Pair 
j  or  I  ( Phugoid  i 


Symmetric 

Modes 


Lateral 


or 


Asymmetric 

Modes 


High  Frequency 
Root  Pair 
(Short-Period 
Model 


-  Small  Real 
Root 

(Spiral  Model 


Large  Real 
Root 
(Roll 

Subsidence  I 


Root  Pair 
( Dutch  Roll  I 


Period  (sec) 


0.138  V„  (fps) 
or 

0.234  V„  (knots) 


Damping  Ratio 
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CiXz  C,,  Ct„ 
C,v 


(tbla  ia  a  relatively  poor  approxiBatloD) 


c 

'»  2V 


Notation  is  defined  in  table  18  . 

Stability  axes  are  the  axes  of  reference. 

[lamping  ratio  is  the  ratio  of  damping  to  critical  damping. 


APPROXIMATE  FORMULAS  FOR  RESPONSE  TO  CONTROL  INPUT 

Approximate  solutions  for  the  response  to  control  input  and  for  maximum  accelerations  are  useful  for  preliminary  esti¬ 
mation  and  checking  of  vehicle  motion.  Several  items  are  included  below  that  provide  estimates  of  response  to  control 
deflection,  maximum  acceleration,  and  roll  rates.  Notation  for  the  relations  in  this  Section  is  given  in  table  18 

1.  LOAD  FACTOR  DUE  TO  CONTROL  DEFLECTION 

In  horizontal  symmetric  steady  flight  the  derivative  of  normal  acceleration  (load  factor)  with  -respect  to  control  deflec¬ 
tion  is  given  approximately  by 

dn  V-  C-«  C**  _  (7 

dS  g  C-,„  C„.  r  —  CIU  (Ci  c  4-  2t) 

2V„  2V„ 
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2.  MAXIMUM  ROLL  VELOCITY 


The  maximum-rolling-velocity  approximation  is  obtained  from  the  linearized  equation  for  moments  about  the  X-axis, 
which  is  solved  for  the  steady-state,  single-degree-of-freedom  caae.  The  approximate  equation  for  maximum  steady 
rolling  velocity  is  then 


(78) 


3.  MAXIMUM  ANGULAR  ACCELERATION 


Maximum  angular  accelerations  resulting  from  control  actuation  are  sometimes  needed  in  aircraft  design  work.  A  simple 
relation  for  maximum  angular  acceleration  is  given  below. 

First  it  is  assumed  that  the  applied  moment  coefficient  is  represented  by  the  function  illustrated  in  the  sketch  below. 


A  single-degree-of-freedom  approximation  may  be  used;  however,  static  stability  is  neglected  in  the  case  of  yawing  and 
pitching  motion.  Under  these  conditions  the  maximum  angular  acceleration  occurs  at  t  —  t,  and  is  given  for  the  rolling 
case  by  the  equation 


Pm  ax 


c‘-  r 

C,,«.  L« 


(79) 


Similarly,  the  equations  for  maximum  pitching  and  yawing  accelerations  are 


C-. 

Cm,  t« 
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C'Y 

C.r  t, 

c,s 


(80) 

(81) 
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SECTION  7.  SPECIAL  PROBLEMS 


INSTRUMENT  READINGS 

In  the  analysis  and  the  automatic  control  of  vehicle  motion  it  is  frequently  desirable  —  or  even  necessary  —  to  utilise 
several  types  of  instrumentation.  Instruments  may  be  used  to  indicate  the  attitude  of  a  vehicle  or  to  measure  velocity 
and  acceleration  components. 

The  following  Sections  present  relations  and  equations  that  are  useful  in  the  resolution  and  interpretation  of  instrument 
readings.  These  relations  are  adapted  from  references  9  and  19 . 


ATTITUDE-MEASURING  INSTRUMENT  READINGS 

Vehicle  attitude  is  usually  determined  from  a  system  of  gyro  instruments.  The  rotation  angles  of  the  instrument  about 
the  inner  and  outer  gimbal  axes  are  related  to  the  vehicle  orientation  angles  *P,  0,  and  4>.  The  conventional  free 
vertical  and  directional  gyro  system  is  shown  in  figure  25. 


VERTICAL  GYRO 
SPIN  AXIS 


VERTICAL  GYRO 


°D„ 


DIRECTIONAL  GYRO 


FIGURE  35  CONVENTIONAL  FREE  VERTICAL  AND  DIRECTIONAL  GYRO  SYSTEMS 

From  the  above  figure  and  the  relations  for  the  orientation  of  vehicle  body  axes,  figure  19  and  table  11,  the  gyro 
equations  below  are  developed. 
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Vertical  Gyro 

Gv,  *=  sin (cos  t|  sin  0  -f-  *>n  «i  cos  *  co»  ®) 
sin  4*  cos  0 


tan 


- s 

V  COS  Cf  cos  <t> 


COS 


0  —  si 


sin  «i  sin 


hT5") 


(82) 


Directional  Gyro 

C(),  _  sin'1  [cos  «,  (cos  T  cos  4>  sin  0  +  sin  V  sin  4>)  +  sin  *,  cos  V cos  0] 

sin  'P  cos  4>  —  sin  4>  sin  0  cos  4* _  ~1 

s  T  —  sin  <i  (sin  *P  sin  4>  -f  cos  »P  cos  4>  sin  0  J 


Gd  =  tan 


'[■ 


COS  <1  cos 


0  cos 


(83) 


where 

GV(,  G0,  are  rotation  angles  about  the  inner  gimbal  axes  of  the  vertical  and  directional  gyros,  respectively. 

GV(),  G„o  are  the  rotation  angles  about  the  outer  gimbal  axes  of  the  vertical  and  directional  gyros,  respectively. 

„  is  the  angle  between  the  outer  gimbal  axis  of  the  vertical  gyro  and  the  vehicle  X-axis  (this  is  also  the 

angle  between  the  vehicle  Z-axis  and  the  directional  gyro  outer  gimbal  axis).  The  subscript  1  denotes 

reference  to  instrument  axes. 

T',  0,  $  are  the  orientation  angles  of  the  vehicle  body  axes  as  defined  on  page  11. 

The  cyro  systems  shown  in  figure  25  and  analyzed  in  equations  (82)  and  (83)  are  for  the  most  simple  form  of  free 
gyro.  More  complex  attitude-  and  direction-sensing  instrumentation  is  used  in  many  advanced  vehicles.  The  output 
indications  of  gyro- instrumented  stable  platforms  in  terms  of.  vehicle  attitude  are  derived  in  reference  33.  The 
attitude  output  signals  of  other  fire-control  and  navigational  devices  are  discussed  in  references  34,  35,  and  36. 


VELOCITY-MEASURING  INSTRUMENT  READINGS 


Velocity  components  are  generally  measured  by  instruments  that  are  not  located  at  the  vehicle  center  of  gravity.  In 
addition,  the  orientation  of  these  instruments  may  not  coincide  with  the  vehicle-orientation  reference  axes  (body  axes). 
Thus,  even  after  instrument  e.rors  and  position  errors  (sidewash,  upwash,  etc.)  are  accounted  for,  the  velocity  com¬ 
ponents  of  the  vehicle  center  of  gravity  are  not  given  directly  by  these  instruments. 

A  general  set  of  instrument  axes  may  be  used  having  its  origin  located  by  a  vector  r  from  the  vehicle  center  of  gravity. 
These  axes  may  be  oriented  with  respect  to  the  vehicle  body  axes  system  by  the  angles  0„  and  4>t  as  indicated  in 
figure  26. 


The  velocity  vector  on  instrument  axes  is  given  by  the  equation 

V,  ~  V-HttX  r 

In  this  equation  the  separate  vectors  are  given  by  the  relations 

Vi  =—Ui  i|  +  Vi  ji  -f-  W[  k| 

V  —  Ui  +  Vj  -f-  Wk 
«>  -  Pi  +  Qj  +  Rk 

r  «=*  xi  -j-  y j  “I- 


( Instrument  velocity! 
(Vehicle  linear  velocity) 

( Vehicle  angular  velocity) 

( Instrument  location  vector) 


(84) 


(85) 
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Hie  general  notation  in  these  equations  is  defined  on  page  9,  and  the  subscript  1  denotes  instrument  axes. 

Expressing  equation  (84)  in  Cartesinn  form  and  applying  the  transformation  matrix  for  Case  I  of  table  2  to  de¬ 
scribe  the  orientation  of  the  instrument  axes  results  in  the  equations  below  for  the  velocity  components  along 
iastrssmct  axes. 


Z 


OBDEE  Or  SOTATION8 


FI  CUBE  M  QEMBBA1,  DfSTBUMBlfT  AXES 


u,™  (U  —  Ry  -f  Qt)  cos 4i  cos  in  +  ( V  —  Px  4-  R*)  cos  0i  sin  —  (W  —  Qx  4-  Py)  »»» 

V!  —  <U  —  Ry  4~  Q*)  (cos^i  sin  <h  sin  $x  —  sin  cos  <fc) 

-|-  ( V  —  Pt  -f  Rx)  (sin  sin  sin  $i  +  cos  cos 
(W  —  Qx  -|-  Py)  sin*,  cos  0j 

Wj  —  (U  —  Ry  -4-  Q»)  (cos  *i  coo  *i  sin  0i  4~  •»  sin  *,) 

4-  (V  —  Pi  4-  Rx)  (>»n  *i  cos  ’  sirs  0,  —  cos  *i  sin  *j) 

4-  (W  —  Qx  4-  Py)  cos  0|  cos  *i 
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The  foregoing  equations  can  be  greatly  simplified  for  moot  cases  of  interest  Usually  and  9i  will  be  sero  and 
certain  of  the  x,  y,  and  s  may  be  negligible.  Also  the  limitation  to  small  disturbances  permits  the  following 

simplification  of  the  expressions  for  the  angles  of  attack  and  sideslip,  respectively. 


M 


a.-uur’^i)  « 

*  “ tan''  (tit)  ~ 

ACCELERATION-MEASURING  INSTRUMENT  READINGS 

The  acceleration-measuring  instruments  are  located  and  oriented  in  much  the  same  way  as  the  velocity  instruments, 

i.e.,  displaced  from  the  vehicle  center  of  gravity.  Therefore  the  instrument  accelerations  must  be  related  to  the  vehicle 
center-of-gravity  acceleration. 

Linear  accelerations  along  instrument  axes  are  given  by  the  equation 

Vi-V  +  «,XV  +  iXt  +  -»X(aXr)  (88) 

The  vectors,  in  addition  to  those  given  in  equation  (85),  are  expressed  below. 

•  •  •  w 

Vi  —  Ui  ij-f  Vji+Wk,  (Instrument  acceleration) 

V  —  UI  -)-  Vj  -f  Wk  (Vehicle  linear  acceleration) 

„  •  *  • 

«*  —  PS  -(-  Qj  -j-  Rk  (Vehicle  angular  acceleration) 

Acceleration  components  along  instrument  axes  arc  expressed  below.  These  equations  are  obtained  from  equation 
(88)  and  Case  1  of  table  2.  The  general  notation  is  defined  on  page  9,  and  the  subscript  I  denotes  quantities 
referred  to  instrument  axes. 


U,~  [U- 
+  [V- 
-[W 

Vi-[U- 
+  [V- 
+  [W 

W,-[U 
+  [V- 
+  [W 

The  above  relations  can  be  simplified  in  moat  instances,  aa  was  the  cam  with  the  previous  velocity  component  equations. 

There  are,  however,  additional  factors  to  be  considered  in  acceleration-measuring  devices  Time  factor*  are 

1.  Effect  of  gravity  force  on  suspended  or  pivoted  mass  (seismic  element). 

2.  The  dynamics  of  the  device  itself;  the  instrument  has  spring  and  damping  forces  acting  an  the  muf ended  mass. 

3.  The  effect  of  acceleration  of  the  instrument  mounting,  when  a  pivoted  mam  ie  wed.’ 

Detailed  analyte*  of  the  above  items  are  beyond  the  scope  of  the  present  report.  The  reader  is  referred  to  die 

presentations  in  references  19  and  37  for  these  details. 


•  RV  -f-  QW  —  x  (R*  -f-  Q*)  +  Y  (PQ  ~  R)  -M(RP  +  Q)]cos*,cos*  | 

PW  -f  RU  -f  *  (PQ  4-  R)  -y  P*4  R*)  4*  *  (QR  —  P)  ]  cos  I,  am  pi 

-  QU  -f  PV  +  x  (PR  —  Q)  +y  (QR  +  P)  -  *  (Q*  4-  P*)  J 

.  RV  +  QW  —  x  (R*  +  Q*)  4-y  (PQ-R)  +*(RP  +  Q)]  icosfi  sin*,  sin#,  -  siafccoe*,) 
PW-f-RU  -f  x  (PQ  +  R)  -  y  (P*  4-  R*)  +*  (QR  -  P)]  (sin  *,  sin  *,  sin  #,  4- eos*,cos*i) 

-  QU 4-  PV  4-  x  (PR  —  Q)  4-  y  (QR  4-  P5  -*  (Q*4-P*)]  •»*«»*  I 

-  RV  4-  QW  -  x  (R*  4-  Q*)  4-  y  (PQ  —  R)  4-  «  (RP  +  Q)  J  («*#!  cos  sin  ^  4-  fi) 

-  PW  4- RU  4- x  (PQ  4-  R)  -  y  (P*  4-  R*)  4-  *  (QR  -  ft]  (sm  *  cos  sin  -  cos  *,  sin  *,) 

-  QU  4-  PV 4-  x  (PR  -  $)  4-  y  (QR  4-  P)  -  ■  (Q*  +  P*)]  c« ^ 


•9 


FUEL  SLOSH 


Fuel  slosh  within  partially  filled  tanks  is  known  to  affect  the  dynamics  of  manned  aircraft  and  missiles.  Fuel 
slosh  introduces  additional  degrees  of  freedom,  owing  to  the  relative  motion  of  the  fuel  mass  and  the  airframe. 
This  Section  of  the  report  is  concerned  with  the  effects  of  fuel  slosh  on  the  rigid-body  modes  of  vehicle  motion. 
The  literature  on  fuel  slosh  contains  numerous  references  to  the  effect  of  foe!  slosh  on  the  flutter  problem.  The 
latter  is  not  treated  here.  Furthermore,  in  some  applications,  such  as  flexible  boosters  for  large  ballistic  missiles 
or  space  vehicles,  fuel-slosh,  rigid-body,  and  body-bending  effects  may  all  be  inseparably  coupled.  The  approach 
described  in  this  Section  would  obviously  require  extensions  to  include  the  effects  of  structural  flexibility  in 
such  cases. 


CONDITIONS  UNDER  WHICH  FUEL  SLOSH  HAS  BEEN  FOUND  TO  BE  SIGNIFICANT 


The  addition  of  fuel-slosh  degrees  of  freedom  to  the  equations  of  vehicle  motion  complicates  the  analysis,  as  may 
be  seen  subsequently.  As  a  guide  to  the  need  for  this  complication,  a  brief  summary  is  presented  in  table  27  of 
conditions  under  which  fuel-slosh  effects  have  been  found  to  be  significant. 


TABLE  27 

SOME  CONDITIONS  UNDER  WHICH  FUEL  SLOSH  HAS  SIGNIFICANT 
EFFECTS  ON  VEHICLE  RIGID- BODY  MODES  OF  MOTION 


Fuel  Muss 
Total  Mms 

Fuel -Tank 
Location 

Fuel-Tank 

Shape 

Fuel-Slosh 

Natural 

Frequency 

Vehicle 

Rigid-Body  ;Vlode* 
Affected 

Remarks 

Ref 

>  0.25 

Forward  of 
Vehicle  c.g. 

Any 

Approximately 
equal  to 

rigid-body  mode 

lateral -directional 
(Dutch- Roll) 
oscillation 

The  mode  damping  is 
reduced  in  this  case. 
Unstable  roots  can 
appear,  leading  to 
a  limit  cycle  (snaking) 

38 

>0.10 

Any 

Large  spanwise 
dimension 

Not  applicable 

Spiral  mode,  in 
horiaontal  Sight 

Spiral  divergence 
occura,  as  if  the  vehicle 
had  negative  dihedral. 

The  long-term  response  to 
directional  control  is 
reversed 

38 

>0.10 

Any 

Large 

longitudinal 

dimension 

Approximately 
equal  to 
rigid -body  mode 

Long-period 
(phugoid) 
oscillation,  in 
horiaontal  flight 

The  mode  damping  is 
reduced  in  this  case 

39 

>  0.25 

■  . . . .  J 

Forward  of 
vehicle  c.g. 

Any 

Approximately 
equal  to 

rigid-body  mode 

Yaw  or  pitch 
oscillation  in 
low-speed  vertical 
flight,  as  for  rocket 
take-off 

The  mode  damping 
is  reduced 

•8s«  p«f«  84. 

As  is  implied  in  table  27,  fuel-sloah  coupling  with  the  short-period  longitudinal  mode  of  rigid-body  motion  is 
usually  negligible  in  horizontal  flight.  Fuel-sloah  effects  may  be  negk  cted  in  studies  of  this  mode.  Coupling 
exists,  of  course,  for  this  mode  of  motion  in  vertical  flight,  as  in  the  take-off  phase  of  many  liquid-fueled, 
rocket-powered  vehicles. 
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The  B'ieady-sUte  relationship  between  fuel  mass  center  shift  and  fuel-tank  acceleration  for  closed-top  rectangular  fuel 
tanks  is  given  reference  38.|These  results  show  that  only  slight  tank  accelerations  can  produce  near-maximum  fuel  mass 
center  shifts,  in  many  practical  cases.  As  an  example,  for  a  typical  height-to-length  ratio  of  0.08  for  a  wing  fuel  tank, 
80  percent  of  the  maximum  possible  fuel  mass  shift  is  attained  for  a  lateral  acceleration  of  0.1  g,  with  the  tank  half  full. 

It  is  concluded  that  the  coupled  motions  of  the  airframe  and  sloshing  fuel  masses  are  generally  significant  for  small 
disturbances.  For  large  vehicle  disturbances,  involving  large  values  of  fuel-tank  acceleration  in  a  horizontal  plane,  slosh¬ 
ing  fuel  tends  to  act  as  an  off-center  fixed  mass,  without  dynamic  coupling  to  the  airframe,  or  with  discontinuous  coupling. 

MOTION  OF  SLOSHING  FUEL 

In  the  analysis  of  fuel-slosh  effects  on  vehicle  flight  dynamics,  sloshing  fuel  masses  are  generally  represented  as  mass- 
spring-damper  single-degree-of-freedom  dynamic  elements.  The  natural  frequencifc  of  the  analog  elements  correspond 
to  the  lowest  or  fundamental  modes  of  fuel  slosh. 

More  complex  analogs  could  be  constructed  to  represent  higher  frequency  fuel-slosh  modes  in  addition  to  the  fur, da- 
mental  mode.  The  forces  applied  to  the  airframe  by  the  higher  frequency  modes  of  fuel  slosh  ere  relatively  small.  As  a 
general  rule,  only  the  fundamental  mode  is  represented  in  practice. 


Available  data  on  the  fundamental-mode  natural  frequencies  for  several  tank-shapes  are  summarized  in  figure  28. 
The  data  in  this  figure  are  in  dimensionless  form.  For  convenience,  an  auxiliary  chart  is  presented  as  figure  29, 
for  the  fundamental-mode  natural  frequencies  of  open-top  rectangular  tanks,  in  terms  of  physical  dimensions. 


Available  data  on  the  forces  applied  to  the  airframe  by  sloshing  fuel  is  more  limited  at  present  than  the  corre¬ 
sponding  data  for  natural  frequency.  For  the  purj>ose  of  this  report,  the  applied-force  data  are  presented  in  terms 
of  the  ‘effective*  fuel  mass.  Hie  effective  fuel  mass  mf  iu  defined  in  relation  to  the  single-degree-of-freedom 
analog  of  figure  27-  If  the  actual  fuel  tank  and  the  analog  are  given  the  same  horizontal  and  rotational  motions, 
tne  effective  fuel  mass  mf,  equal  to  the  concentrated  mass  in  the  analog,  provides  the  same  reactions  on  the 
container  (acting  through  the  spring  and  damper)  as  the  reactions  on  the  fuel  tank  applied  by  the  sloshing  fuel. 

For  reference,  open-top  rectangular  data  of  figure  30  were  developed  from  the  equivalent  pendulum  concept  of 
reference  40.  Evaluation  of  the  transfer  function  relating  applied  force  to  linear  and  rotational  input  tank  motions 
showed  that  the  effective  fuel  mass  mf  is  equal  to  the  equivalent  pendulum  mass  of  reference  40. 


SLOSHING  FUEL 


FIGURE  2T  SINGLE-DEQBES-OF*  FREEDOM  ANALOG  TO  SLOSHING  FUEL 
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aOKIEONT  VL  CIBCUL  \B  CYLINDRICAL.  TAJIK  (r«f.  4 


UPRIGHT  CIRCULAR  C  TLINDRICAL  TANK  <»»?.  4 1  ) 


FIGURE  3  8  (CON  TO.) 
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AS  AN  APPROXIMATION.  USE  THE 
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FIGURE  30  EFFECTIVE  MASS  OF  FUEL  SLOSH  FUNDAMENTAL  MODE  rOR  VARIOUS  FUEL-TANK  SHAPES 
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